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[57] ABSTRACT 

A pulse-echo, immersion method for ultrasonic evaluation 
of a material which accounts for and eliminates nonlevelness 
in the equipment set-up and sample thickness variation 
effects employs a single transducer and automatic scanning 
and digital imaging to obtain an image of a property of the 
material, such as pare fraction. The nonlevelness and thick- 
ness variation effects are accounted for by pre-scan adjust- 
ments of the time window to insure that the echoes received 
at each scan point are gated in the center of the window. This 
information is input into the scan file so that, during the 
automatic scanning for the material evaluation, each 
received echo is centered in its time window. A cross- 
correlation function calculates the velocity at each scan 
point, which is then proportionalized to a color or grey scale 
and displayed on a video screen. 

15 Claims, 5 Drawing Sheets 
(1 of 5 Drawings in Color) 
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PULSE-ECHO ULTRASONIC IMAGING 
METHOD FOR ELIMINATING SAMPLE 
THICKNESS VARIATION EFFECTS 

ORIGIN OF THE INVENTION 

The invention described herein was made by an employee 
of the United States Government and may be manufactured 
and used by or for the Government for governmental pur- 
poses without payment of any royalties thereon or therefor. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The invention relates to ultrasonic evaluation of material 
properties. More particularly, the invention relates to non- 
destructive ultrasonic evaluation of materials by measuring 
velocity using a single transducer pulse-echo immersion 
system, automatic scanning and digital imaging, which 
provides a video image of the sample in color or grey scale 
which is a map of a material property such as porosity 
fraction. 

2. Background Of the Disclosure 

Nondestructive evaluation applicable to evaluating prop- 
erties of materials such as ceramics, metals, plastics and 
various composites are known to those skilled in the art and 
include x-radiography, ultrasound or ultrasonic evaluation, 
and thermal methods. These methods provide an efficient, 
quasi-quantitative measure of material homogeneity, but 
often lack the precision necessary for microstracture evalu- 
ation of high-performance materials, such as high tempera- 
ture oxidation resistant ceramics and the like. The develop- 
ment and use of materials for high-performance applications 
requires detailed, quantitative knowledge of microstructural 
and compositional variability for defining acceptable levels 
of variability and for rejecting those materials and processes 
that yield sample-to- sample and within-s ample variations 
likely to result in unacceptable property (e.g., strength, 
thermal conductivity, oxidation resistance, resistance to 
spalling, etc.) variations. Such variability must be precisely 
characterized either directly in terms of property measure- 
ment or indirectly through microstructural characterization 
where microstructure-property relations have been previ- 
ously established. 

Repeated, uniformly spaced ultrasonic contact measure- 
ments have been successful for quantifying and mapping 
inhomogeneity in various ceramics (e.g., SiC, A1 2 0 3 , 
YBa 2 Cu 3 0 7 and Si 3 N 4 ) and metals in terms of ultrasonic 
material properties such as reflection coefficient, velocity 
and attenuation coefficient as mentioned, for example, by 
Roth, et. al. in Quantitative Mapping of Pore Fraction 
Variations in Silicon Nitride Using an Ultrasonic Contact 
Scan Technique , NASA TP 3377 (1993). This publication 
describes quantitatively characterizing material (e.g., Si 3 N 4 ) 
microstructure in terms of actual ultrasonic wave param- 
eters. The wave parameters include reflection coefficient, 
attenuation coefficient and velocity. A post-scan interactive 
data display system is used for comparing ultrasonic prop- 
erties at different locations within samples and viewing the 
resultant ultrasonic images. Further refinement of this pro- 
cess is disclosed by Roth, et al. in PSIDD: A Post-Scan 
Interactive Data Display System for Ultrasonic Scans , 
NASA TM-4545 (1993). This process relates to contact 
scans and does not disclose how to account for thickness 
variations in the sample being measured. Piche discloses a 
single transducer immersion method for evaluating plastic 
using a technique in which 16 scan points are pulsed for the 
sample and the results evaluated using regression analysis 
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[L. Piche, Ultrasonic Velocity Measurement for the Deter- 
mination of Density in Polyethylene , Polymer Eng. &. Sci,,v. 
24, n.17, p. 1358-58 (Dec. 1984]. This method does not 
relate to forming an image of the sample property, nor does 
5 it provide an experimental technique that automatically 
accounts for nonlevelness and thickness variation during a 
scan procedure required to form an image. Consequently, a 
need still exists for a method which will permit ultrasonic 
material evaluation that will account for nonlevelness and 
10 thickness variations in the material, require only a single 
transducer, eliminate problems associated with physical 
contact between the transducer and sample or buffer rod, and 
display, on a video screen in gray scale or color, an image of 
the scanned material which is a map of an internal structural 
15 property of the material, such as porosity fraction. 

SUMMARY OF THE INVENTION 

The invention relates to a method for nondestructive 
ultrasonic evaluation of materials by measuring velocity 
20 using a pulse-echo immersion system with automatic 
scanning, echo cross-correlation and digital imaging to 
obtain a grey scale or color image of the sample. The 
velocity values obtained for each scan point are scaled on a 
grey or color scale and displayed on a video screen which 
25 shows a material property, such as porosity fraction. Prior to 
the automatic scanning, nonlevelness in the set-up and 
sample thickness variation effects are accounted for and 
eliminated by insuring that the echoes at each scan point are 
first gated and input into a scan parameter file in a computer, 
30 so that during the subsequent automatic scanning each 
received echo is centered in the time window set for it 
While it is possible, but not practical to do a manual prescan 
at each and every scan point needed for a two dimensional 
video image of the material property being evaluated, many 
35 sample thickness variations are in the form of a uniform 
thickness variation from one edge to another. In this case of 
uniform thickness variations from one edge to another, 
preliminary scans are performed along a single line in both 
the x- and y-directions of the sample to provide slant 
40 correction factors. The slant correction factors are input into 
the scan parameter file so that any wedge-shape variations 
are taken into account during the automatic scanning for the 
material evaluation, to insure that each echo received during 
the automatic scanning is centered within the time window. 
45 A single transducer is used in a preferred embodiment of the 
invention. 

In the immersion method of the invention, the material to 
be evaluated is surrounded by a liquid and positioned over 
an acoustic reflector which is also immersed in the liquid. An 
50 ultrasonic wave of a known frequency is transmitted through 
the liquid and four separate echoes are recorded and evalu- 
ated at each scan point. Each echo is received as an analog 
waveform which is digitized and stored in a computer. The 
echoes received, digitized and stored during the sample 
55 evaluation scans are the first two succesive echoes reflected 
off the back surface of the sample, the first echo reflected off 
the front surface of the accoustic reflector in which the 
received wave has passed through the sample, and the forth 
is the first echo reflected off the front surface of the reflector 
60 with the sample not present, so that the received wave does 
not pass through the sample. This means that at least two 
separate scans must be made, with and without the sample 
present between the transducer and reflector. However, as a 
practical matter it is difficult from both a hardware and 
65 software perspective to accomplish this in just two scans and 
obtain maximum time resolution and thus maximum accu- 
racy. Consequently three or four separate scans are 
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performed, with three being faster and four being more 
accurate. The choice is left to the discretion of the practi- 
tioner. In the embodiment in which four separate scans are 
performed during the sample evaluation, the ultrasonic wave 
goes through both the liquid and the material during the first 
three scans. For the forth scan the material sample is 
removed so that the transmitted wave is reflected off the 
front surface of the reflector without going through the 
material. Although the order is not important, it is conve- 
nient to receive the first echo reflected off the back surface 
of the material during the first scan and the second succes- 
sive echo reflected off the back surface of the sample during 
the second scan. During the third scan in which the trans- 
mitted wave goes through both the immersion liquid and the 
material sample, the first echo reflected off the front surface 
of the reflector is received. The first echo reflected off the 
front reflector surface is received during the fourth scan 
when the sample is not present This process is repeated at 
a plurality of scan points sufficient to produce a video image 
of a microstructural property, such as porosity, of the mate- 
rial. After the scanning is completed, the digitized wave- 
forms are retrieved from the computer and the time delay 
between the first two successive echoes received from the 
back surface of the material at each scan point is determined. 
The time delay between the two different reflections or 
echoes received off the reflector (with and without the 
transmitted wave going through the sample) is also deter- 
mined for each scan point The wave velocity at each scan 
point is then calculated from the time delays and the speed 
of the transmitted wave in the liquid. The velocity values for 
all of the scan points are scaled to corresponding propor- 
tional color or grey scale values which are then displayed on 
a video screen or cathode ray tube (CRT). Thus, in this 
embodiment of the invention, four separate scans are made 
at each scan point to separately receive, as analog 
waveforms, the first two successive ultrasonic echoes off the 
back surface of the sample and the first echo off the front 
surface of the reflector both with and without going through 
said sample; digitizing and storing the waveforms; retrieving 
the digitized waveforms; determining the received time 
delay between the first two successive sample back surface 
echoes and between the two reflector front surface echoes; 
calculating the wave velocity at each scan point; scaling the 
calculated velocities to corresponding proportional color or 
gray scale values, and displaying the resulting image. The 
wave velocity at each scan point is calculated from 

v=c (-#- +l ) 

wherein 2t is the received time delay between the two 
successive sample back surface echoes, wherein At is the 
time delay between the two different echoes received from 
the reflector with and without going through the sample, and 
wherein c is the speed of the transmitted ultrasound wave in 
the liquid. The embodiment in which three separate scans 
are made is similar to that in which four separate scans are 
made, with the difference being that during the first scan the 
first two succesive reflections or echoes off the back surface 
of the sample are received, digitized and stored. It will be 
noted that above equation does not include the sample 
thickness value. This means that the thickness of the sample 
need not be measured or known. 

As set forth above, prior to the two, three or four scans 
during which the sample is evaluated, nonlevelness and 
sample thickness variations are accounted for and eliminated 
by pre-scans to insure that the received reflections or echoes 
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are within their set time windows to provide a complete 
waveform for evaluation and cross-correlation to accurately 
obtain the time delay data used in calculating the velocity 
values. In the case of a sample having a thickness variation 
5 in the form of a uniform thickness variation from one edge 
to another, preliminary scans are performed along a single 
line in both the the x- and y-directions of the sample to 
provide slant correction factors. The slant correction factors 
are input into the computer scan parameter file so these 
io variations are taken into account during the automatic scan- 
ning for the material evaluation. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The file of this patent contains at least one drawing 
15 executed in color. Copies of this patent with color 
drawing(s) will be provided by the Patent and Trademark 
Office upon request and payment of the necessary fee. 

FIG. 1 schematically illustrates the spatial relationship 
between the transducer, liquid, material sample, reflector 
20 plate and the transmitted and reflected ultrasonic echo waves 
in the practice of the invention. 

FIGS. 2(a) and 2(b) graphically illustrate the amplitude 
and time delay of the received analog ultrasonic echoes 
25 reflected off the material sample and reflector. 

FIGS. 3(a) and 3(b) graphically illustrate respective first 
and second sample back surface echoes in which the second 
is inverted with respect to the first. 

FIG. 4 is a block diagram schematically illustrating the 
30 instrumentation used in the invention. 

FIGS. 5(a) and 5(b) are photographs of video grey scale 
displays of a thickness based ultrasonic velocity image of a 
ceramic according to the prior art method and the method of 
the invention, respectively. 

55 DETAILED DESCRIPTION 

Referring to FIG. 1, transducer 10 is schematically shown 
as partially immersed in a liquid 12 which is the immersion 
fluid. Material sample 14 is shown positioned in the fluid 12 
40 between the transducer 10 and reflector plate 16 in container 
or tank 18 located on top of x- and y-direction motorized 
stages 19 and 20. Supports 15 and 15' maintain sample 14 
above reflector 16. The transducer 10 and motorized x- and 
y-direction stages 19 and 20 are electrically connected by 
45 means not shown to a pulser-receiver (not shown) and to 
means (not shown) for moving the x- and y- stages in their 
respective directions. Similarly, liquid 12 is connected to 
means not shown for maintaining the temperature of the 
liquid preferably within at least about ± 1° F. of the tem- 
50 perature at which the scans are to be run. It is possible to 
maintain the liquid temperature within ±0.1° F. The better 
the temperature control, the more accurate the results will 
be. For example, if the temperature of the immersion liquid 
is ±1° F. and the liquid is water, a 1.5% error in velocity is 
55 possible. If the porosity fraction or other property of the 
material at a particular point is such as to result in a velocity 
value difference in the sample of 2%, only a 0.5% micro- 
structural velocity difference might be detected if a ±1° F. 
temperature variation is present during the scan. At each 
60 scan point an ultrasonic wave 22 of a known frequency is 
transmitted from transducer 10 through liquid 12 and into 
material sample 14. Entering material 14 causes part of wave 
22 to be reflected (not shown) off the top surface of the 
sample, with the rest of the wave passing through the 
65 material as 23. Part of wave 23 continues through the 
material and to the top surface of accoustic reflector 16 as 
24, is reflected back off the top surface of reflector 16 as 25, 
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passes back through the sample 14 and returns to the 
transducer 10 as wave 26. A portion of wave 23 is reflected 
off the back surface of the sample and returns to the 
transducer as 27. Part of the wave 23 reflected off the back 
surface of the material is reflected off the top surface, returns 5 
to the back surface, is again reflected back to the top surface 
and exits as wave 28. Waves 27 and 28 are the first two 
successive back surface reflected waves used in the method 
of the invention at each scan point. Not shown is the wave 
transmitted through the liquid and reflected back to the 10 
transducer without going through the material. This wave 
which is not shown and wave 26 are the two reflector front 
surface echoes used in the method of the invention. Motor- 
ized stages 19 and 20 form part of an automated scanning 
system which incrementally moves in both the x- and 15 
y-directions to obtain an ordered array of points across the 
entire surface of the material sample. A 20 MHz, broadband 
transducer was used in the practice of the invention. Broad- 
band transducers emit a broadband frequency content domi- 
nated by a center frequency. That is, they are made to emit 20 
at a nominal frequency proximate that of the design fre- 
quency (e.g., 20 MHz), with a Gaussian fall-off on either 
side of the nominal center frequency. Thus, a 20 MHz 
broadband transducer will also emit frequencies slightly 
above and below the nominal center frequency of 20 MHz. 25 
In the Piche article referred to above, although the two 
different reflector front surface echoes are captured and 
recorded, the first front surface echo and the first back 
surface echo are captured and recorded. This is different 
from the method of the invention which captures and records 30 
the first two successive sample back surface echoes and not 
the first front sample surface echo. Further, Piche does not 
use automatic sample scanning or digital imaging. FIGS. 

2 (a) and 2 (b) graphically illustrate the reflected waveforms 
received and displayed on the CRT of an oscilloscope as 35 
time domain analog waveforms. Turning to FIGS. 2(a) and 
2(b ) , the intensity or strength of the received waveform is 
displayed as voltage amplitude, which is the ordinate of the 
graph, and the received time delay as the abscissa. In this 
representation, Bl, B2 and M' refer to waves 27.28 and 26 40 
of FIG. 1, respectively, with M" representing the wave 
transmitted through the liquid and reflected off the front 
surface of the reflector without passing through the material 
sample. The time delay between the first two successive 
echoes reflected from the back or bottom surface of the 45 
material back to the transducer, Bl and B2, is readily 
obtained, as is the time delay between the two reflections 
received from the front surface of the reflector, M’ and M". 
Since the velocity of the ultrasonic wave is faster in denser 
media than in less dense media, voids, delaminations, poros- 50 
ity and other density variables within the material are 
obtained as a function of the speed of the wave, which is 
determined by the time delay between the first two succes- 
sive echoes received which have been reflected off the back 
of the material, and the time delay between the two different 55 
reflections from the front surface of the reflector. As set forth 
above under SUMMARY, the speed or velocity of the 
transmitted wave traveling through the material sample is 
determined according to the simple equation: 

60 

v=c (-i- +1 ) 

wherein 2 t is the received time delay between the two 
successive material sample back surface echoes, wherein at 65 
is the time delay between the two different echoes received 
from the reflector with and without going through said 
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sample, and wherein c is the speed of the transmitted wave 
in the liquid. This equation is accurate for a single point 
measurement Prior art ultrasonic velocity scan techniques 
such as that of Roth et. al. in the quantitative mapping 
publication referred to above, assume that the material 
sample is of uniform thickness and do not take into account 
nonlevelness and material thickness variations as does the 
method of the invention. 

In the practice of the method of the invention, tank 18 may 
be made of any suitable material. Gear plastic such as 
polymethylmethacrylate (e,g., Lucite or Plexoglass) has 
been found useful. The sample tank contains a suitable 
elastic liquid, such as water, as the immersion fluid to 
provide an accoustic coupling between the transducer, mate- 
rial and reflector plate. Since the x-, y-direction scans made 
across the sample surface in the method of the invention can 
take a significant amount of time compared to that for a 
single point measurement and since the speed of sound in a 
liquid is also a function of temperature, the water is main- 
tained at a constant temperature during the scanning. This is 
readily accomplished simply by using a constant tempera- 
ture regulating means, such as a constant temperature water 
circulator, for maintaining the desired temperature constant 
during the ultrasonic scanning. It is convenient to keep the 
temperature of the water at about ambient or 68° F. ±1° F. 
during the scan, although other temperatures may be used if 
desired, as long as the temperature is maintained within no 
more than ±1° F. In the case of distilled, deionized water, the 
wave velocity may be obtained from published tables. 
However, tap water may be used as long as the velocity in 
the water is actually measured. The reflector is placed on the 
bottom of the tank. Other immersion liquids may be used, if 
desired, such as Dow Corning 704 vacuum pump oil. 

The reflector is a solid plate of material having an 
accoustic impedance significantly different from that of the 
liquid or water. A flat plate of tungsten (e.g.. Vie" -Vi" thick) 
is preferably used, because tungsten has an accoustic imped- 
ance almost two orders of magnitude higher than water in 
units of g/cm 2 -sec. The use of a tungsten plate results in the 
highest possible reflection amplitude of any solid material 
for the echoes reflecting off the front surface of the reflector 
plate. This large difference in accoustic impedance is impor- 
tant when attempting to obtain ultrasonic echoes that have to 
travel into and through immersion liquid and the sample, 
bounce off the reflector plate, and travel back through the 
liquid and sample to the the transducer for reception. High 
frequency ultrasound provides greater time resolution than 
lower frequencies and is therefore more desireable for 
greater accuracy of the velocity of the ultrasound through 
the sample and corresponding velocity image. The higher 
the frequency, the greater the velocity accuracy. By having 
the highest reflection amplitude possible, it is possible to use 
the method of the invention (a) at higher frequencies where 
attenuation through the sample is greater than if using lower 
frequencies and (b) with materials that significantly attenu- 
ate ultrasound, such as composite materials. By high fre- 
quency ultrasound is meant from 1-100 MHz, typically 
3-50 MHz and more typically 10-30 MHz. 

The material sample is easily positioned over the reflector 
plate by using spacers on top of the plate and placing the 
material on top of the spacers. It is important that the spacers 
have the same height or thickness so that the material is as 
level as possible. Lucite is available as sheets which are very 
uniformly thick and it is convenient to use 0.5" cubes of this 
plastic as spacers. The material sample, such as a plate of 
silicon nitride ceramic, is placed on the plastic spacers over 
the tungsten reflector plate prior to scanning. 

FIG. 4 schematically illustrates, in block diagram fashion, 
the basic system and instrumentation used for the scanning 
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and ultrasonic imaging according to an embodiment used in 
the practice of the invention. 

Thus, referring to FIG. 4 , the basic instrumentation 
includes a transducer 50 , a pulser-receiver 52 , and a pro- 
grammable waveform digitizer 54 having associated with it 
a vertical voltage amplifier 56 , programmable time base 58 , 
and anolog and digital monitors 60 and 62 , respectively. 
Also included are a time delay or synthesizer 64 , an image 
processor 66, an X, Y, Z controller 68, computer or central 
processing unit (CPU) 70 , and video display 72 . In the 
embodiment shown, the computer 70 is a CPU with terminal 
74 and associated video display 76 also forming part of the 
system. Monitors 60 and 62 , along with digitizer 54 , voltage 
amplifier 56 and time base 58 also serve as respective analog 
and digital oscilloscopes. The time synthesizer, time base, 
voltage amplifier and waveform digitizer are all general 
purpose interface [IEEE- 488]bus (hereinafter “GPIB”) pro- 
grammable and interconnected via GPIB cables. The com- 
puter 70 is programmed in Fortran and contains an image 
processor system which is connected to the video color and 
gray scale display 72 . The computer controls the GPIB 
instrumentation and acquiring of the desired waveforms via 
the GPIB. The process includes data acquisition, analysis/ 
calculation, image processing and display. The Fortran 
software, with callable routines in in IEX-VMS interface 
software to communicate with the GPIB instruments is 
written for instrument control and waveform acquisition 
following the method of Generazio, et. al. in Interfacing 
Laboratory Instruments to Multiuser Virtual Memory 
Computers , NASA Technical Memorandum 4106 (1985), 
the disclosure of which is incorporated herein by reference. 
The Fortran programs used in the practice of the invention 
including the scanning program, the analysis and cross- 
correlation program, the grey scale imagemaker program, 
and the display program are contained in the attached 
Appendix. The wave form digitizer is a Tektronics 7912 AD 
Programmable Digitizer along with a Tektronics 7A16 P 
Programmable Vertical Voltage Amplifier (voltage base) and 
a Tektronics 7B90 P Programmable Time Base. The time 
delay (time synthesizer) is a Hewlett Packard Model 5359A 
and the X, Y, Z programmable stepper motor controller and 
associated tables is a Klinger Scientific C-1.22. The ultra- 
sonic pulser-receiver is a Panametrics Model 5601 and the 
transducer used in the scanning of the silicon nitride ceramic 
disk in the example below is a Panametrics 20 MHz, 
longitudinal, unfocused, broad band transducer. The com- 
puter is a Digital Equipment model Microvax IL 

In the practice of the invention, the image processer is a 
Grinnell Systems Grinnell 274 Image Processing System 
and the Grinnell Systems GMR Series Software Package, 
Release 2.2, Jun. 19,1981, available from McLoud 
Associates, 165-F Croftich Lane, Campbell, Calif. 95008, 
the disclosure of which is also incorporated herein by 
reference. Two video displays are used, one of which is a 
DEC VT340 terminal, which is the user terminal attached to 
the computer, and a Mitsubishi 20LP is the video display 
monitor attached to the image processing system. High level 
VAX Fortran software used for driving the system is 
included in the Appendix as set forth above. The Grinnell 
library of Fortran subroutines is called from this high level 
software. The video display shows ultrasonic images. 

The pulser-receiver applies the voltage pulse to the trans- 
ducer to generate the ultrasonic waves into the sample and 
to the reflector plate and also receives the raw ultrasonic 
echo waveforms from the transducer. The approximate times 
where the echo waveforms are expected to occur are deter- 
mined a priori (prior to the automatic scanning for the 
material evaluation) using the time systhesizer to find and 
position the echo waveforms on the oscilloscope. The time 
base and voltage amplifier are used to modify the time and 
voltage scales to view the waveforms on the oscilloscope. 
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Both the time base and the time synthesizer are externally 
triggered by the pulser-receiver (a+2 volt synchronizing 
pulse). Triggering occurs on the positive slope of the pulse. 
The time base is adjustable over the range of from 1 
5 psec-500 msec/div on the oscilloscope, with the optimum 
setting for each waveform determined a priori and input to 
a data file in the computer. The pulser-receiver output is 
connected to the voltage amplifier. The voltage amplifier, 
selectable over the range 50 mV-lV/div, is automatically 
adjusted by the digitizer so that the entire received analog 
io waveform is digitized with maximum amplitude fit onto the 
digital waveform monitor. The digitizer digitizes each wave- 
form received into 512 point arrays (at a sampling rate 
ranging from 0.512 -1.024 GHz depending on the time base 
time/division setting). Each waveform is acquired 64 times 
15 and averaged to obtain a smoother waveform with averaged 
noise levels using a Fortran algorithm included in the 
scanning program in the Appendix and which is also found 
in the NASA Technical Memorandum 4106 referred to 
above. The X and Y positional and Z intensity outputs from 
the waveform digitizer are attached to the analog and digital 
monitors 60 and 62 . The analog monitor is used for the 
prescans and the digital for the automatic scanning. 

As set forth under the SUMMARY, prior to the two, three 
or four scans during which the sample is evaluated, nonlev- 
elness and sample thickness variations are accounted for and 
25 eliminated by pre-scans to insure that the received reflec- 
tions or echoes are within their set time windows to provide 
a complete waveform for evaluation and cross-correlation to 
accurately obtain the time delay data used in calculating the 
velocity values. That is, during the nonlevelness and mate- 
30 rial thickness variation scans, the operator notes if the time 
delay of each echo received at each scan point is such that 
it is no longer centered within the oscilloscope time window. 
If a received echo is not centered within the time window on 
the scope, this is noted and the time window changed for 
35 each such echo received until the received echo time domain 
waveform is completely within the new time window set for 
it to insure that the complete time domain waveform is 
captured or gated completely within the new window. This 
time delay information at each scan point is inputted into the 
scan parameter file and recalled during the actual scanning 
40 during the material evaluation, to automatically adjust the 
time delay for the received echoes at each scan point so that 
each echo received during the scanning is centered within 
the time window set for it. This is very time consuming to 
do for each scan point. However, in the case of a sample 
45 having a thickness variation in the form of a uniform 
thickness variation from one edge to another, preliminary 
scans are performed along a single line in both the the x- and 
y-directions of the sample to provide slant correction factors. 
The slant correction factors are input into the computer scan 
50 parameter file so these variations are taken into account 
during the automatic scanning for the material evaluation. It 
is important that the echo at each scan point is centered in 
its time window, because the whole pulse or echo time 
domain waveform is needed to give the precise time delay 
55 between echoes for the cross-correlation which provides the 
velocity value. In doing this for a wedge shaped sample, 
prior to the automatic scan, the transducer scans along two 
straight lines over the sample, once in the the x-direction and 
once in the y-direction, during which an operator notes the 
echo received from the first and last scan points, starting 
60 from the first scan point which is generally at one comer of 
the area defined for scanning. The time difference from the 
sample end-to-end in each of the x- and y-directions of the 
first and last scan point is noted by the operator who then 
adjusts the time base for each echo if needed to insure that 
65 it is centered within the time frame set for it This is done so 
that so that each echo received during the scans in which the 
material is being evaluated is centered (gated) within the 
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oscilloscope time window set for it so that the received 
waveform is displayed with the maximum possible ampli- 
tude on the CRT and still have the complete waveform. This 
permits the maximum time resolution of individual echoes 
without losing any part of the time domain waveform which 
appears on the CRT screen as a function of voltage 
(amplitude) and time, wherein time is the x- axis and voltage 
or amplitude is the y- axis. It is important and forms an 
aspect of the invention that the complete waveform or pulse 
echo be captured or “gated” on the CRT screen in order to 
perform an accurate cross-correlation later on in the proce- 
dure of the process of the invention. The cross-correlation of 
echoes provides the precise time delay between received 
echoes or pulses which is required to calculate the velocity 
or speed of the ultrasound in the material evaluated which, 
in turn, provides the information to gray or color scale the 
velocity data into a digitized map of the material density. 
This slant correction procedure also allows an accurate 
evaluation to be made without the need for specialized 
leveling equipment. These x- and y-direction time window 
corrections are called slant correction factors and they are 
inputted into the scan parameter file in units of “nsec/pm” 
where (a) the number of nsec is the time extent from sample 
end-to-end that is required to keep the specific echo centered 
and is determined using the (a) time synthesizer to reposition 
echoes in time and (b) the number of pm is the distance 
traveled by the transducer for which this slant factor is 
determined. By way of an illustrative, but nonlimiting 
example, the first scan point ( 0 . 0 ) along the x- direction may 
have a Bl echo centered at a time=6.77 psec, while the last 
scan point ( 40 , 0 ), may have Bl centered at time=7.14 psec. 
If the x- direction scan line length is 40 mm, the x-direction 
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wherein W^is the correct delay time window at a particular 
scan location, T 7 is the time delay at the the initial scan 
5 location, X sc and ?sc are the x- and y-direction slant 
correction factors, X SN and Y SN are the scan point numbers 
in the x- and y-directions, and X SI and Y SI are the x- and 
y-direction scan increments. With many samples it has been 
found that the slant correction factors turn out to be the same 
10 for the Bl ,B2 echoes and the slant correction factors for the 
M', M" echoes are the same. However, for some samples 
(e.g., thick samples), they may not be the same. In such cases 
a first x- and y-direction scan is made for the Bl echoes and 
1S a second x- and y-direction scan made for the B2 echoes. 
The same holds for the M' and M" echoes for which two 
separate scans are made in the x- and y-directions. 

A scan parameter file is input into a computer which 
contains all of the information necessary to automatically 
20 scan the material sample being evaluated. This information 
includes a predefined and ordered array of scan points over 
which to run the scan. By way of an illustrative, but 
nonlimiting example, in an example of the method of the 
invention in which the material being evaluated was a 
25 monolithic ceramic wedge, the scan consisted of a 41 
(X-direction) by 81 (y-direction) grid of measurements for a 
total of 3,200 scan points, with each measurement or scan 
point separated by 1 mm (x-) and y- scan increment). 
Information input into the scan parameter file (NOTHICIC_ 
ALLSHAPE1.DAT) includes the following: 


C **TTILE** NOTHICK ALLSHAPE1DAT 
C ** SCAN INCREMENT (uM) IN X-DIRECTION IS: 

1000. 

C ** SCAN INCREMENT (uM) IN Y-DIRECTION IS: 

1000. 

C ** SCAN LENGTH (uM) IN X-DIRECTION IS: 

40000. 

C ** SCAN LENGTH (uM) IN Y-DIRECTION IS: 

80000. 

C ** X-DIRECTION SLANT CORRECTION FACTOR (nsec/uM) FOR Bl & B2 ECHOES IS: 

-0.0055 

C ** Y-DIRECTION SLANT CORRECTION FACTOR (nesec/uM) FOR REFLECTOR ECHOES IS: 
-0.0055 

C ** Y-DIRECTION SLANT CORRECTION FACTOR (nsec/uM) FOR Bl & B2 ECHOES IS: 

-0.00175 

C ** Y-DIRECTION SLANT CORRECTION FACTOR (nsec/uM) FOR REFLECTOR ECHOES IS 

0.0 

C ** TIME LOCATION (uSEC) OF Bl ECHO AT SCAN ORIGIN IS: 

52.83 

C ** TIME LOCATION (uSEC) OF B2 ECHO AT SCAN ORIGIN IS: 

52.31 

C ** TIME LOCATION (uSEC) OF REFLECTOR ECHO W/SAMPLE PRESENT AT SCAN LOCATION 
IS: 

69.46 

C ** TIME LOCATION (uSEC) OF REFLECTOR ECHO W/O SAMPLE PRESENT AT SCAN LOCATION 
IS: 

72.48 

C ** IMMERSION FLUID VELOCITY (cm/uSEC) IS: 

0.148 

C ** B2 PHASE-INVERTED WRT Bl (Y/N)?: 

N 

C ** M" PHASE INVERTED WRT M' (Y/N)?: 

N 


slant correction factor is obtained from (7.14-6.77)/40 psec/ 
mm. It should be noted that slant correction factors can be 
negative as well as positive numbers. The location of the 
time window during scanning for the material evaluation is 
automatically adjusted via computer control by using the 
formula: 


As set forth above, in the method of the invention, the 
transducer is activated so that the first front surface echo off 
the sample, the first two successive back surface echoes, and 
the first echo off the reflector plate are all seen in the 
oscilloscope display at the same time by adjusting the the 
time base to the appropriate time per division setting and 
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adjusting the time synthesizer delay time. Viewing the first 
front surface echo off the sample enables the operator to 
know if the back surface echoes are also on the CRT screen. 
The unfocused transducer is positioned above the sample at 
a distance determined initially by the natural focal distance. 5 
When using an unfocused transducer a good initial starting 
height is approximately one to two inches above the sample. 
The reflector plate front surface echo may be low in ampli- 
tude compared to the sample back surface echoes, so that the 
the pulser-receiver gain/attenuation or vertical amplifier gain 10 
settings may have to be increased to see this echo. It is 
important not to confuse the echoes off the front surface of 
the reflector plate with the second set of echoes originating 
from the front and back surfaces of the sample. The second 
set of echoes originating from the front and back surface of 15 
the sample will always occur at twice the delay time where 
the first set of these echoes appears. For example, if the first 
set of echoes begins at 50 msec on the digital oscilloscope, 
the second set will begin at 100 msec. If using, for example, 
a three milimeter thick sample placed on 0.5" thick plastic 20 
supports on the reflector plate, the first reflector echo will 
occur at about 20 msec after the time where the first set of 
echoes originates and thus the reflector echo will be seen at 
about 70 msec in this illustration. Another way to note the 
reflector plate echo is to raise and lower the sample while 25 
noting the location of the stationary echo corresponding to 
the stationary reflector plate. It is essential to have reflector 
plate echoes that will not interfere with the second set of 
echoes originating from the front and back sample surfaces. 
Attention is next focused on the first back surface echo from 30 
the sample, Bl. The echo is centered in the oscilloscope time 
window to obtaining maximum time resolution by adjusting 
the time base time per division and the time synthesizer 
delay. The synthesizer time is recorded and inputted into the 
scan parameter computer file. This procedure is repeated for 35 
the second back surface sample echo, B2, the first front 
surface echo off the reflector plate with the sample present, 

M* and the first echo off the reflector plate with the sample 
removed, M". The next step is to account for and eliminate 
any nonlevelness in the set-up and also sample thickness 40 
variations. This done at each scan point, except that in the 
case of uniform thickness variations in the sample, the slant 
correction factors outlined above are determined and input 
into the scan parameter file in the computer. 

The scanning is then automatically performed through the 45 
remainder of the scanning points previously inputted into the 
scan parameter file in the computer using a program written 
in Fortran and IEX GPIB to perform the scanning and also 
to obtain maximum vertical voltage resolution of the 
received ultrasonic waveforms. Scanning is accomplished 50 
through the use of computer controlled x-, and y- micros- 
canning tables used to reposition the sample in the x- or 
y-direction in a 1 mm increment (other increments may be 
used at the convenience and discretion of the practitioner) 
for the next measurement. The ultrasonic waveforms 55 
received are then digitized (512x512 pixel resolution) at 
each scan location and stored successively in the scan data 
file in the computer. Four separate ultrasonic scans are 
performed at each scan location. As set forth above, the 
echoes are Bl (first echo off sample back surface, obtained 60 
in first scan), B2 (second echo off sample back surface, 
obtained in second scan), M* (first echo off front surface of 
reflector plate with the sample present, obtained in third 
scan), and M" (first echo off front surface of reflector plate 
without the sample present, obtained in fourth scan). Each of 65 
the four echoes is obtained in a separate scan to obtain the 
maximum time resolution for each echo by setting, before 


12 

each scan during the nonlevelness and material thickness 
variation procedure, the optimum time per division setting 
on the oscilloscope time base that allows maximum time 
resolution. The minimum number of scans for this thickness- 
elimination procedure is two, but the time per division 
setting for only two scans cannot be obtained in this case as 
the time per division setting would be fixed for all three 
echoes obtained in the first of the scans using this scan 
procedure. 

The following is an algorithm of a scanning program 
which accounts for and eliminates the nonlevelness of the 
set-up and uniform thickness variation effects of the sample 
in the resulting ultrasonic image displayed on the video, the 
code for which is included in the Appendix. 

1) Determine the scan lengths and scan increments in the 
x- and y-directions, time positions of echoes at scan 
origin, slant correction factors, and immersion fluid 
velocity. 

2) Edit NOTHICK _J\LLSHAPE1.DAT FILE, which is 
the scan parameter file, and input information from 1) 
above. 

3) Start scanner fortran program on computer which 
automatically does the following: 

A) Initialize all GPIB instrumentation, which includes the 
time synthesizer, digitizer, time base, voltage amplifier, 
Klinger X, Y stages. 

B) Perform scan to digitize Bl echoes and store in file 

I) Digitize Bl at scan origin 

Adjust voltage base for echo with maximum 
amplitude in video/oscilloscope window 
Move Klinger tables under transducer in x- direc- 
tion specified x- direction increment 
Time synthesizer moves to delay time position 
determined by Bl, B2, slant correction factors. 
This results in echo in video being centered in 
the Tektronics analog video/oscilloscope dis- 
play and subsequently digitized and stored. 

Time 

where T 0 =correct delay time window at a 
particular scan location S^=x- direction slant 
correction factor (nsec/pm) 

N^=scan point number in x- direction 
Lf=x- direction scan increment (pm) 

Sy=y- direction slant correction factor (nsec/pm) 
Ny=scan point number in y- direction 
Iy=y- direction scan increment (pm) 

II) Repeat I) until one scan line in x- direction is 
completed. 

HI) Increment transducer in y- direction specified y- 
direction increment and repeat I)— H). 

IV) Repeat I)— HI) until y- scan length is traversed 
and scan is completed. 

V) Return Klinger tables to scan origin. 

C) Perform scan to digitize B2 echoes and store in file 
by repeating steps B(I-V) 

D) Perform scan to digitize reflector echoes with 
sample present and store in file by repeating steps 
B(I--V), but using reflector echo slant correction 
factor 

E) Remove sample. Perform scan to digitize reflector 
echoes without sample present and store in file by 
repeating steps B(I-V), but using reflector echo slant 
correction factor 

4) Start velocity calculation Fortran program on computer 
to produce a file of velocities at each scan location by 
performing the cross-correlation algorithm. 
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5) Start image formation Fortran program on computer 
which results in a file of values between 0 and 255 
which scale directly with the velocity values. 

6) Start image display program which brings grey scale 
level image up on video. 

Before initiating the scanning procedure, the temperature 
of the water or other immersion fluid is measured. If the fluid 
is water, published tables or graphs of temperature and 
velocity can be used to determine the velocity of the 
ultrasound in the constant temperature water bath. If the 
immersion liquid is a liquid other than water, or if a more 
precise temperature than that available in published graphs 
and tables is desired, the velocity of the ultrasound in the 
liquid is determined by recording the times (T P ) where 
ultrasonic peaks occur for two different vertical positions 
(Z1 and Z2) of the transducer above the reflector plate. The 
velocity, V, is then determined from 

V={Zl-Ziy(T P \-T P 2) 

The phase relationships of (a) B1 compared to B2 and also 
(2) the reflector front surface echo with the sample present 
(M') compared to that without the sample present (M") are 
examined. These phase relationships are important for the 
computation of the velocity image of the scanned sample. 
The quantity 2 t is obtained by cross-correlating echoes B1 
and B2 which is defined as the precise time delay between 
the B1 and B2 echoes. If B1 and B2 are phase inverted with 
respect to each other, the time occurrence of the minimum 
in the cross-correlation function is used to obtain 2 t. If M* 
and M" are phase inverted with respect to each other, the 
time occurrence of the minimum in the cross-correlation 
function is used to obtain At. Otherwise, at the time 
occurence of the maximum in the cross-correlation function 
is used. 

FIGS. 3 (a) and 3(b) graphically illustrate the case in 
which B2 is phase inverted with respect to Bl. The same 
holds for M' and M". Phase relationships generally remain 
the same throughout the scan, unless significant discrete 
micro structural defects are encountered by the ultrasonic 
wave. 

After the scan has been completed and all the received 
echoes have been digitized and stored in the scan data file in 
the computer, they are recalled from the data file to perform 
the velocity image calculation for each scan location. In 
performing this cross-correlation, an overlap method is used 
by the computer based on a cross-correlation program using 
Fast Fourier transforms published in pages 415 and 416 
(Correlation and Autocorrolation Using the FFT) in the book 
Numerical Recipes — The Art of Scientific Computing , by 
Press, et al., 1988 Edition, Cambridge Univ. Press.. The 
Fortran program used is in the Appendix. Echoes M’ and M" 
are also cross-correlated to obtain At where where M' is the 
echo reflected off the reflector plate front surface with the 
sample present, M" is the echo reflected off the reflector 
plate front surface without the sample present, and At is the 
time delay between them. If M' and M" are phase inverted 
with respect to each other, the time occurrence of the 
minimum in the cross-correlation function is used to obtain 
At. Otherwise the time occurence of the maximum in the 
cross-correlation function is used. The velocity, V, at each 
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scan location is then calculated from the equation referred to 
above. The velocity value for each scan location is sequen- 
tially stored in the computer. After the scan is completed the 
velocity values are scaled on a gray or color scale with a 
5 value directly proportional to the velocity values, with the 
highest and lowest scale values corresponding to the highest 
and lowest velocity values. 

The invention will be further understood with reference to 
the example below. 

10 

EXAMPLE 

In this example, a sample of silicon nitride ceramic was 
evaluated using the thickness based velocity image method 
disclosed in the NASA Technical Memorandum TP 3377 
15 referred to under Background. This method is based on a 
velocity, cross-correlation ultrasonic imaging method with- 
out the pre-scan to account for and eliminate nonlevelness in 
the set-up and sample thickness variations. In this method, 
only the first two sample back surface echoes are captured 
20 and evaluated. The silicon nitride ceramic was 3.5 mm thick 
with a uniform 300 micron thickness gradient. Very coarse 
time scaling was used so that the Bl and B2 echoes stayed 
in the time window while the sample thickness changed as 
the scan proceeded. 

25 The same silicon nitride ceramic sample was also scanned 

and velocity imaged on a grey scale according to the method 
of the invention which included the prescans to eliminate 
set-up and sample thickness variations and which also 
captured and cross-coirolated both the first two successive 
30 sample back surface echoes and the two different reflector 
front surface echoes. 

In both cases, the 20 MHz broad band transducer was 
used, the immersion liquid was water and the back plate was 
35 tungsten as set forth above. 

FIGS. 5(a) and 5(b) are photographs of video grey scale 
displays of the thickness based ultrasonic velocity image of 
a ceramic according to the prior art method, and an image 
according to the method of the invention which included the 
4 q prescans, respectively. Referring to FIG. 5(a), it is seen that 
the top defect is masked due to that part of the sample being 
thicker than the bottom part Also, the defect near the bottom 
is not too discernable and the lower portion is very light due 
to it being thinner. In marked contrast and as shown in FIG. 
45 5(b), the method of the invention clearly and correctly 
illustrates the defect areas, including resolution of the upper 
defect and an overall porosity gradient in the sample. It is 
believed that this demonstrates the efficacy and improve- 
ment to the art of the invention. 

50 It is understood that various other embodiments and 
modifications in the practice of the invention will be appar- 
ent to, and can be readily made by, those skilled in the art 
without departing from the scope and spirit of the invention 
described above. Accordingly, it is not intended that the 
55 scope of the claims appended hereto be limited to the exact 
description set forth above, but rather that the claims be 
construed as encompassing all of the features of patentable 
novelty which reside in the present invention, including all 
the features and embodiments which would be treated as 
equivalents thereof by those skilled in the art to which the 
invention pertains. 
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nothick SCAN.FOR —4 SCAN METHOD 
Scan and data aquisition for the modified ultrasonic Immersion scanner 
by Don J.Roth 

PLEASE USE CLNEW.COM TO COMPILE & LINK 


Latest Update 20-scp-94 

integers SCANSTEP,INDXSTEP 
integer^ SCANAX1S, INDXAXIS, ZSTEP 
INTEGERS 117 

character VOLANS* l ,TSCHEME* l ,UOLANS* 1 
CHARACTER LOOPANS*! ,SHAPEANS* 1 ,ZiGANS* 1 ,FLAG* l 
character TBUF*8, SCHEME* 1 ,FILENAME*34,SCANMODE* 1 
real DELAY(0:4),VOLTSET(0:4),MAXVALF,MAXVALT,VOL(4) 
real slant Lslant2.slant3,slant4 

integers NSCAN, N1NDX,PLACE,ZIGC0UNT,JXXX,SCANKNT 
integers A(5 1 2).NA VE, CHECKP, WWFLAG 
INTEGER*2 1 15,1 16.I26,I27.MAXFS1P,MAXF$1 PT 
INTEGERS 136,137.146,147 

1NTEGER*2 L00PTHRUJI,L00PTHRU_S,LLL,1T.I,JJJ,II 
BYTE WRK0(80) 
common /1BLK/ WRKQJBIJFFER 
common /SBLK/ A, DELAY, VOLTSET,MAI 
data CHECKP/O/ZIGCOUNT/1/ 

DATA 16/1/1 1 5/1/1 1 6/1 /[ 1 7/1 /PLACE/O/ 
data 126/1/127/1/136/1/137/1/146/1/147/1/ 

10010 fonnat(A) 

10012 format( A2) 

10020 format( I) 

10030 format(F) 

OPEN( unit=S ; file-TXAl:', status=‘NEW’ ) 

CALL ENTERPARAM( SCANSTEPJNDXSTEP, NSCAN,NINDX,NAVE,VOLTAGEL 
1 ,VOLTAGEU,VOLTAGE_MAX,SHAPEANS,LOOPANS,LOOPTHRU_H.LOOPTHRU_S, 

1 ZIG ANS,SCHEME,TSCHEME,UZEROO,FlLEN AME,I 1 5,SCANMODE,SLANTl,SLANT2 

1 , SLANTS, SLANT4) 

4328 SCANAXIS - 1 
INDXAXIS = 2 
ZSTEP - 10 

CALL STRTGP1B 
DO 1 800 1 ADDR= 5, 1,-1 
1 800 CALL INIT1NSTR( 1ADDR ) 

C 

C Setup of 7912: begin with V/D ^.5, then automatically find 
C the best intensity and Digitize Defects 

C 

C NOTErMAY NEED TO CHANGE THIS V/D TO S TART 
C 

2000 CALL TEKGTL 

* *** ****** ****** * **** ********** ********** ***** ************* 


I 
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10 FORMAT (3F) 

READ (1 0,65403) DELAY l INOTHICK - STARTING 1ST BACK SURFACE ECHO 

WINDOW TIME 

READ (10,65403) DELAY2 INOTHICK - STARTING 2ND BACK SURFACE ECHO 

WINDOW TIME 

READ (10,65403) DELAY3 INOTHICK - STARTING (W/SAMPLE) REFLECTOR ECHO 

WINDOW TIME 

READ (1 0,65403) DELAY4 INOTHICK - STARTING (WO/SAMPLE) REFLECTOR 

ECHO WINDOW TIME 

READ (10,65402) VOLANS 
READ (10,65403) VOL(l) 

READ (10,65403) VOL(2) 

READ (1 0,65403) VOL(3) 

READ (10,65402) UOLANS 
READ (10,65403) VOL(4) 

65402 FORMAT (A) 

65403 FORMAT (F) 

51 FORMAT(A32) 

52 FORMAT(A2) 

53 FORMAT(A4) 

54 FORMAT(A2) 

TYPE *, 'DELAY TIMES (B1,B2,RS,RNS)’ 

TYPE VDELAYl “'.DELAY 1 
TYPE *,'DELAY2-.DELAY2 
TYPE *;DELAY3=’,DELAY3 
TYPE *,'DELAY4-, DELAY 4 
TYPE ' 

CLOSE (10) 

DELAY(1)=DELAY1/(1 *10**6.) 

DELAY (2)=DELAY2/( 1 * 1 0 * *6.) 

DELAY (3)=DEL AY3/(1 * 1 0.**6.) 

DELAY (4)=DEL A Y4/( 1 * 1 0 «6.) 

TYPE*,’ 1 
TYPE YDELAY 
TYPE ' 

q*** ********* ****** ********************************************* 

CALL IN1T1NSTR(0) 

CALL SETVOLTDIV(O.S) 

CALL TEKGTL 
CALL GETMAl(MAI) 

CALL TEKGTL 
CALL PUTTIME(DELAY(1)) 

CALL TEKGTL 
CALL GETBESTMAI( MAI ) 

CALL D1GDEF 

! Preliminary Digitizations to get Watch level and Timcsct (Tiinc/div) 

CALL 'TEKGTL 

CALL AUTOSETVOLTS( MAI, VOLTSET(l),NUMINT ) 

C TYPE 'AUTOSET VOLTS’ 

IF( VOLTSET( l).EQ. 999. )THEN 

TYPE * ' DIGITIZER IS SCREWED UP 

GOTO 2000 

ENDIF 


a 
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CALL GETGR]D{ TIMESET, X ) 

CALL GETSA( NAVE,MAI,A ) 

CALL TEKGTL 

TYPE *, , Timeset=V TIMESET,' Delay-', DELAY(l), 

+ ' Voltset=',VOLTSET(l) 

TYPE VI 16= Ml 6,' 1 1 7= ',11 7 

C*********************** 1 ,******,( e ,:*,, )(# **,:**** (! **, t , [ * 3; , t! , )) ,, t **., # 

WRITE(1 6,rec~l 1 6,fmt=53)T!MESET 

116-116+1 

C 

C Get WATCH - standard ground level 
C 

TYPE*,’ 1 
type*; ’ 

TYPE*; ’ 

TYPE*; f 
W = A(1)/NAVE 
WATCH- W- 100. 

WR1TE( 8.12030 )WATCTI 

type*; * 

TYPE*, 1 * 
type*; ’ 

TYPE *; ’ 

12030 format( ' Minimum acceptable ground level =',F ) 


C SCANSTEP- 1000*SCANSTEP 

C INDXSTEP = 1000* IN DX STEP 

CALL SETXYZ( SCANSTEP, INDXSTEP, ZSTEP ) 

WR!TE( 5,12140) 

12140 formate fit ### SCANNING ###'//) 

IDIR^l 

WWFLAG=0 

do 2600 noth l =1,4 ! nothickness mod 

if (nothJ.ge.2)Lhen 

CALL MOVORGXY (ZIGANS,SCANAXIS,NSCANJNDXAXIS,NINDX) 
REWTND(12) ! .DATHS FILE 

if (zigans.eq;y')zigcount=i 

ENDIF 

IF (NOTHJ.EQ.4)THEN 
type *,’ 

type *; 1 

type *,’ * 

type *; ******* Finished w/ Sample Sean ************ » 

type +; *** Remove Sample from tray & Hit <RET> to perform water scan’ 

type *,' 

type ’ 

type *; * 

Read (5,10020) IO 
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ENDIF 

DO 2400 1=1 .NINDX ! Outer loop over index 

DO 2200 J=1,NSCAN ! Inner loop over scan 

IF (NOTH_l.ge.2)THEN 

READ(12)JJ,ILFLAG,SCANKNT !READ WHETHER 'H' OR ’S' 

type ’ 
type f 
type ' 

type *, ’second time around for\jj,ii 
type */FLAG -,FLAG 
type V ’ 
type V * 

ENDIF 


IF ((SCANMODE.EQ.'P'.OR.SC ANMODE.EQ.'L’).AND.I.GT. 1 )GOTO 2400 !MOD FOR LINE / POINT 

SCAN 

IF (NOTH JLEQ. I )TYPE */I= ',1,' AND J=\J 
PLACE=PLACE+1 

WRITE( S, 13000 )J, I, PLACE ! 

************************************+*********** :************£ 

IF (SHAPEANS.EQ.’N')THEN 

CALL TAKEDATA_0(FILENAME > VOLANS,UOLANS,VOLJI 6,1 1 7,126,127, 

I 136,137,146, 147, NAVE, WATCHJ,JaiothJ,SLANTI , 

1 SLANT2,SLANT3,SLANT4, SC ANSTEP,INDX STEP, ZIGANS,N SCAN) 

ELSEIF (SHAPEANS.EQ/YyrHEN 
IF (NOTH_I.EQ.2.AND.FLAG.EQ/H')GOTO 3429 

CALL T AKEDATA(PLACE,SCHEME,TSCHE ME, VOl AN S, VOL, UZEROO, FILENAME, 

I UOLANS,I26,I27,NOTH_I,I36, 137,146, 147, 

1 SL ANTI, SLANTS, SLANT3,SLANT4,SCANSTEP,INDXSTEP,ZIGANS,N SC AN, 

1 RCMAX,MAXFS1P,MAXVALF,RF1,RF2,MAXFSIPT,MAXVALT,RT1,RT2,I16,I17, 

1 N AVE, WATCH, J, 1,11 5,TT) 

ENDIF 

C* *** ******* ******************************* ****************** **£ 


3429 CONTINUE 

IF (NOTH LEQ.l .AND.IT.LT.4)THEN 

TYPE *,716= ',136,' 1 17= ',117 

ELSEIF (NOTH J.EQT. AND. 11\EQ.4)THEN 

TYPE *,'146= ',146, ' 147= ',147 

ENDIF 

IF( J.EQ.NSCAN )GOTO 2200 ! 

IF (LOOPANS.EQ.'N')GOTO 22001 ! DON'T PICK UP XDUCER EACH TIME 

IF (NOTHJ.EQ.l)TYPE *,TT*\IT 
IF (NOTH_I.EQ.l)TYPE ' 


H 
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IF(N0THJ.EQ.1)TYPE V 1 

IF (NOTHJ.EQ.l)TYPE VVOLANS (WRT "ON SAMPLE")= ’,VOLANS 
IF (NOTH I.EQ.1)TYPE VUOLANS (WRT "NOT ON SAMPLE)= '.UOLANS 
IF (NOTH__I.EQJ)TYPE ’ 


c ************* MOD FOR ZIGZAG SCAN ******************************c 
22001 IF (ZIGANS.EQ. I Y'.AND.ZIGCOUNT.EQ.O)THEN IZJGZAG SCAN 

IF (SCANMODE.EQ.'P')GOTO 12306 (FOR POINT MEASURE, DON'T MOVE X-AXIS 
CALL MO VX YZ(SCAN AXI S 1 ) ! MOVE X-AXTS BACKWARD 
ELSEIF (ZIGANS.EQ.'Y r .AND.ZIGCOUNT.EQ4)TH£N 

IF (SCANMODE.EQ.'P^GOTO 12306 IFOR POINT MEASURE, DON’T MOVE X-AXIS 

CALL MO VX YZ(SC AN A X I S , I ) IMOVF. X-AXIS FORWARD 

ENDIF 

12306 CONTINUE 

CALL WAIT2(1000) !ADD DELAY SO THA I DATA IS NOT TAKEN BEFORE 
C KLINGER STAGES STOP MOVING 

CALL WA1T2(1000) 

CALL WA1T2(1000) 

CALL WAIT2(1O0O) 

CALL WAIT2(1000) 

CALL WA1T2(1000) 

CALL WAIT2(1000) 

CALL WA1T2(1 000) 

CALL WAIT2(1000) 

CALL WAIT2(I000) 

CALL WAIT2(1000) 

CALL WAIT2(1000) 

CALL WAIT2(1000) 

CALL WAIT2(1 000) 

IF (ZIGANS.EQ.'Y’)GOTO 12305 
£*******************¥****************£• 

IF (SCANMODE. EQ.'P’)GOTO 2402 IFOR POINT MEASURE, DON'T MOVE X-AXIS 
2402 CONTINUE 

IF (LOOPANS.EQ. , N')GOTO 2200 
12305 CONTINUE INOTHICK 
2200 CONTINUE 


IF( I.EQ.NINDX )GOTO 2400 ! 

C IFOR LINE SCAN OR POINT MEASURE, DON'T MOVE Y-AXIS 

IF (SCANMODE.EQ. , P'.OR,SCANMODE.EQ. , i:)GOTO 2401 
CALL MOVXYZ( INDXAX1S, 1 ) ! Move 

CALL WAIT2(1000) 

CALL WAIT2(1000) 

CALL WAIT2( 1000) 

CALL WAIT2(IOOO) 
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CALL WAIT2( 1000) 
CALL WAIT2(1000) 
CALL WAIT2(1000) 
CALL WA1T2(1000) 
CALL WAIT2(1 000) 
CALL WAIT2(UK)0) 
CALL WAIT2(1000) 
CALL WAIT2( 1000) 
CALL WAIT2( 1000) 


2401 NTOGO - NINDX-I ! 

IF (ZIGANS.EQ. f N')TIIEN 'NORMAL SCAN 

CALL MOVORG(SCANAXIS,NSCAN,VOLTAGEL,VOLTAGEU, 

1 LOOPTHRIJ_H) ’Move back to the X origin 

£******************* MOD FORZIGZAG SCAN *******************£ 

ELSEIf (Z1GANS.EQ/Y’)THEN 
IF (Z1GCOUNT.EQ. 1 )THEN 
ZIGCOUNT=0 
GOTO 15021 

ELSEIF(ZIGCOUNT.EQ.0)THEN 

ZIGCOUNT-1 

END1F 

END IF 

£** * ***** ********** *** ** * ******* **£ 

15021 CONTINUE INOTHICK 

2400 CONTINUE ! — 

2600 continue 

CALL TIMB(TBUF) 

\VR1TE( 8,13100 )TBUF 
CLOSE( 8 ) 

CLOSE( 6 ) 

CLOSE(14) 

JXXX=L1 !!!!!!!!!!!!!!!!!! Y-LOCATION 

IF (SCAN MODE.EQ. , L\OR.SCANMODE.EQ.’P’)JXXX- 1 !MOD FOR LINE SCAN & POINT 
MEASURE 

WRITEC 1 5,R£X>I 1 5,FMT=52)JXXX 
II 5-1 15+1 
C ENDIF 
STOP 

13000 format( ' ’,13,13,15) 

13100 format( ' ’,A8 ) 
end 


SUBROUTINE TAKEDATA_O(FILENAME > V0LANS > U0LANS J VOL,116, 117,126,127, 
1 136, 137,146,147, NAVE, WATCH, I, Lnoth_I,SLANTl, 

L SLANT2,SLANT3,SLANT4,SCANSTEP,INDXSTEP,ZIGANS,NSCAN) 

integer^ A(512), WFLAG,I1 6,126, WWFLAGJ,J.IL0CX,NSCAN,1L0CY 1 136,I46 
INTEGERS 1 1 7SUB,1 1 7,I27,SCANSTEP,INDX STEP 

INTEGERS 127SUB,I37SUB,I47SUB 


C O 
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real DELAY(0:4), VOLTSET(0:4),VOI ,TAGE(4),VOL(4),SLANTl ,SLANT2,SLANT 
REAL SLANT3,SLANT4 

character DAY*9, TIM*9,VOLANS*l ,UOLAN$* 1 ,C1*80,F1LENAME*34,ZIGANS* 1 
common /SBLK/ A, DELAY, VOLTSET, M A I 
C 

DATA 114/0/ 

C 

if(noth_I.eq.l)then 

OPEN( unit=14, file~'[ROTH. MENUjDONSCAN_INTERP.LOG', 

+ stati^'NEW’jACCESS^SEQUENTIAL 1 , 

+ FORM- UNFORMATTED') 

endif 
C 

51 FORMAT(A32) 

52 FORMAT(A2) 

53 FORMAT(A4) 

54 FORMAT(A2) 

WFLAG = 0 

10600 format(",II5 ) 

if (nothj.eq J)then 

IX=l 

IY-l 

ciseif (noth_l.eq.2)THEN 

IX— 2 

IY-2 

elseif (noth l.cq.3)THEN 

IX=3 

IY-3 

elseif (noth_I .eq.4)TI IEN 

IX=4 

IY*=4 

ENDIF 

DO 600 IT — IX,IY 

IF (IT.EQ. 1 ,OR.IT.EQ.2)THEN 

SLANTX-SLANT1 

SLANTY=SLANT3 

ELSEIF (1T.EQ.3.0R.1T.EQ.4)THEN 

SLANTXR=SLANT2 

SLANTYR-SLANT4 

ENDIF 

C ***** MODIFICATION FORZIGZAG SCAN ******************************0 
IF (ZIGANS.EQ.'Y')THEN 

IF (LEQ.2.0R.I.EQ.4.0R.1.EQ.6.0R.1.EQ.8.0R.I.EQ. 10 
I .OR.I.EQ.12.0R.LEQ. I4.0R.I.EQ. I6.0RXEQ. 1 8.0R.I. 

1 EQ.20.0R.I.EQ.22.0R.I.EQ.24.0R.I.EQ.26.0R.I.EQ.28.0R. 

1 IJEQ30X)RXEQ,32.OR.I.EQ.34.ORJ.EQ.36.OR.l.EQ.3S. 

1 OR.I.EQ.40.OR.I.EQ.42.ORXEQ.44,OR,I.EQ.46.OR.LEQ. 

1 48.0R.I.EQ.50.0R.LEQ.52.0R.I.EQ.54.0R.LEQ.56.0R. 

1 I.EQ.58.0R.LEQ.60.0R.LEQ.62.0R.LEQ.64.0R.1.EQ.66 
I .OR.LEQ.68.OR.T.EQ.70.OR.LEQ.72.OR.I.EQ,74.OR.I. 

1 EQ.76.0R.I.EQ.78.0R.I.EQ.80.0R.I.FQ.82.0R4.EQ.84.0R 


7 
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l .LEQ.86.0R.I.EQ.88.0R.I.EQ.90.0RT.EQ.92.0R.L 
1 EQ.94.0R.I.EQ.96.0R.I.EQ.98.0R.I.EQ.100)THEN 
ILOCX~NSCAN-J IFOR ZIGZAG SCAN, REVERSE BACK FOR DISPLAY 
ELSE 
ILOCX==M 
ENDIF 

elsciRzigans.eq.'N^then Igeorge wood addition 
ILOCX=J-J ! George Wood addition 
ENDIF 
ILOCY=I-I 

IF (IT.EQ. I .OR.IT.EQ.2)THEN ! CORRECT DELAY WINDOW FOR NONLEVELNESS 
DELAY_COR R= DELAY (IT)+(SLANTX * I LOCX* SCANSTEP)+{SLANTY*[LOC Y * I NDXSTEP) 
ELSEIK (IT.EQ.3.0R.IT.EQ.4)THEN 

DELAY_CORR-DELAY(IT)+(SLANTXR*ILOCX*SCANSTEP)+(SLANTYR*ILOCY*JNDXSTEP) 

ENDIF 


TYPE*; • 

TYPE*; ' 

TYPE*; ' 

type *;j=*,j; delayc,it;)= \delay(it) 

type *;DELAY_CORR= \DELAY_CORR 
TYPE*; ’ 

TYPE*; ' 

TYPE*; ' 

CALL PUTTIME{ DELAYCORR ) I Set delay 
100 CONTINUE 

1F((1T.NE.4.AND.V0LANS.EQ;U , ).0R.(ITEQ.4.AND.U0LANS.BQ. , U , ))THEN 
c !!!! USER-DEFINED VOLTAGE SETTINGS 

VOLTS=VOL(IT) 

CALL SETVOLTDIV(VOLTS) 

GOTO 1132 
ENDIF 

CALL AUTOSETVOLTS( MAL VOLTS ,NUMFNT ) ! Set V/D 

U 32 IF(( VOLTS.GT. 1 .0).OR.( VOLTS.LT.O.O 1 ))THEN 

WRITE(5,1 1 3 1) VOLTS 

1131 F0RJV1AT(‘+7BAD VOLTAGE SETTING\E1 0.5) 

CALL TEKRESET(MA1) 

GOTO 100 
ENDIF 

CALL GETSA(NAVE,MAI,A) ! Get waveform 

WATCH1 - A( I)/ NAVE ! Make sure its acceptable 

TYPE V ‘ 

TYPEV 1 
TYPE*; ’ 

WRITEC5, 1 1 33)WATCH I , WATCH 

1133 FORMATC+VCURRENT WATCH LEVEL 1S’.F10J;MIN <F10.5) 

TYPE *; • 

TYPE*,’ 1 


% 
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TYPE*, 1 ' 

IF ((1T.NE.4.AND.V0I.ANS.EQ;U , ).0R.(IT.EQAAND.U0LANS.EQ;U , ))G0T0 56732 
C !!!! USER-DEFINED VOLTAGE SETTINGS 

C IF (VOLAN$.EQ/U')GOTO 56732 !! SKIP WATCH 


IF( WATCH1.LT. WATCH )THEN 
WWFLAG - WWFLAG+1 
IF (WWFLAG.EQ.l)THEN 

OPEN( unit=88, fiIe-tROTH.MENUlDONSCAN.LOG’, 

+ $laLus=’NEW’ ) 

ENDIF 

WRITEC 5.10200) 

CALL TIME(TIM) 

CALL DATE(DAY) 

TYPE V ' 

TYPE*, 1 1 
TYPE*; * 

type *, 'watch 1 =\ watch I.’ < watch-, watch 
type*; ■ 
type*; ' 

TYPE*; ’ 

TYPE*; ' 

WRITEC 88,10201 )FILENAME i DAY,TIM,JJ,Cl 
1 0201 formate \A7 *,a; X- M3,’ Y= U3;MAY LEAD TO BAD PROPERTY VALUE’ ) 

1 0200 foimatC WAVEFORM BELOW "WATCH” ; MAY LEAD TO BAD PROPERTY VALUE’) 

C GOTO 100 

ENDIF 

£***:M****************** ************************* ************ 

56732 CONTINUE 

IF (nothj.eq.l)thcn 

500 WRITE(16,REC=I16.FMT=53)DELAY_C0RR 
116=116+1 

WRITE(I6,REC=II6,FMT=53)VOLTS 

116=116+1 

C 

DO 54321 IJI=1 ,512 
1 1 7SUB=((1 1 7- 1 )*5 1 2)+I JI 
54321 WRITE( 1 7,REC=1 17SUB,FMT=54) A(IJI) 

117=117+1 

elseif(nothJ.eq.2)then 

5500 WRJTE(26,REC=l26,FMT=53)DELAY_CORR 
126=126+1 

WR1TE(26,REC=126,FMT=53)V OLTS 
126=126+1 
C 

DO 54329 1JI-1,512 
I27SUB=((I27-1)*5 12)+IJT 
54329 W R1T E{2 7, REC= 12 7 SUB . FMT=54) A(I J I) 





ooo^o^^^oo 
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127=127+1 

elseif (noth_Leq.3)then 

5800 WRITE(36,REC=I36.FMT=53)DELAY_CORR 
I36-I36+1 

WR1TE(36.REC=J36.FMT=53)V0LTS 

136=136^-1 

C 

DO 58329 UI=l,5t2 
I37SUB=((r37-1 )*5 12)+IJI 
5S329 WRJTE(37,REC=137SUB,FMT=54) A(IJT) 

137=137+1 

elseif (noth_l.cq.4)thcn 

5900 W RITE(4 6,REC= 146, FMT=5 3 )DEL A Y_C ORR 
146=146+1 

WRITE(46,REC-i46,FMT=53)VOLTS 

146-146+1 

C 

DO 59329 IJ1=1, 512 
I47SUB=((I47-1 )*5 1 2)+TJT 
59329 W R ITE(47,REC=I4 7 S U B , F MT=54) A(IJ1) 

147=147+1 

endif 

£***************************** ************** ***************£ 

IF (IT .EQ . 1 ) VO LT AG E( l )= VOLT S 
IF (IT.EQ.2)VOLTAGE(2)=VOLTS 
IF (IT. EQ. 3) VOLT AG £(3)= VOLTS 
IF (IT.EQ.4)VOLTAGE(4)=VOLTS 

C* ************************** *************************** ***£ 

600 CONTINUE 

£***************************************************£ 

IF (I.EQ.0.AND.J.EQ.0) GOTO 20000 'SKIP NEXT STEP (NOISE MEASURE) 
CHECK FOR ERROR IN FS2 VOLTSET 

IF ( V OLTAGE( I ).EQ..0 1 ) THEN 
WRITE (14) U,' VOLT FOR FS2= ’,VOLTAGE(l), 

+ '#INT.PO = \NUMINT 
GOTO 100 
ENDIF 

CHECK FOR ERROR IN B1 VOLTSET 

c IF (VOLTA G E( 1 ).LT.VOLTAGE(2)) THEN 

c WRITE (14) U' VOLT FOR Bl = WOLTAGE(l), 

c ^ ' &INT.PO = '.NUMINT 

C GOTO 100 

c ENDIF 

7 FORMAT (A4) 

£****************«******** ********* *************** ** 

20000 CONTINUE 


ID 
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RETURN 

end 

C 
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SUBROUTINE TAKEDATA(PLACE.SCHEME,TSCHBMK,VOLANS,VOL.UZCRCX),FrLENAME, 

1 UOLANS,I26,I27,NOTHJ,136, 137,146,147, 

I SLANT1,SLANT2,SLANT3,SLANT4,SCANSTEP,INDXSTEP.ZIGANS,NSCAN, 

1 RCM AX,M AXF S 1 P,M AXVALF,RF 1 ,RF2, MAXFS 1 PT,MAX VALT,RT1 ? RT2,I 1 6. 

1 117, NAVE, WATCH J,I,i 15,IT) 

intcgcr*2 A 1 1 (5 1 2) > I26,I ) J,CLOCX ; NSCAN, WFLAGJ 1 6,MAXFS 1 PT, PLACE 
INTEGERS WATCII_COUNT I IIITJLOCY,l36,I46 
INTEGER*2 MAXFS IP, NAVE,A(5i2),n\SCANKNT 
INTEGERS n7SUB,n7JSTAT,T27 ; SCANSTEP,lNDXSTEPJ37,147 
INTEGERS I27SUB,I37SUB,147SUB 

REALM ASPECF I( 1 024),PHASE( 1 024),MAXVALF,MAXVALT > VOL(3),SLANT3,SLANT4 

real *4 DELAY(0:4),VOLTSET(0:4).VOLTAGE(4), SLANT, SLANT1,SLANT2,FI(5 12) 

COMPLEX CF1 ( 1 024),CSPEC(1 024) 

character SCHEME* 1 ,FLAG* I ,FILENAME*34,UOLANS* J .ZIGANS* l 
CHARACTER VOLANSMJSCHEMEM 
common /SBLK./ A, DELAY, VOLTSET, MAI 
C 

DATA 1 1 4/O/HJT/O/ 

C 

IF (NOTHJ.EQ.l)THEN 

OPEN( umt=i4, fHeMROTH.MENUjDONSCAN_INTERP.LOG’, 

•I- status =f NEW\ ACCESS-SEQUENTIAL'. 

+ FORM- UNFORMATTED') 

ENDIF 

C 

5 1 FORMAT(A32) 

52 FORMAT(A2) 

53 FORMAT(A4) 

54 FORMAT(A2) 

WFLAG = 0 

C WATCH_COUNT=0 

10600 format( ",I15 ) 

IF (noth l.eq.l ,and.PLACE.EQ.l)THEN 
OPEN( unit=12, ii le=FI LEN AME//TI S'. status^EW', 

+ form- UNFORMATTED’, ORGANIZATION- SEQUENTIAL 1 ) 

SCANKNTH) 

ENDIF 


if (noth_l .eq. I )then 

IX=1 
1 Y— 1 

elseif (noth_I.eq.2)THEN 

IX-2 

1Y=2 

elseif (noth_l.eq.3)THEN 

IX=3 

1Y=3 

elseif {noth_I.eq.4)THEN 

IX=4 

IY=4 

ENDIF 


11 



41 


5,629,865 


42 


10386 DO 600 IT- IX, IY 


C ***** MODIFICATION FOR ZIGZAG SCAN ******************************c 
IF (ZIGANS.EQ.V)THEN 

IF (I.EQ.2.0R.LEQ.4.0R.I.EQ.6.0R.I.EQ 8.0R.I.EQ. 1 0 
1 X3R.I.EQ.12.0R.I.EQ.14.0R.I.EQ.16.0RXEQ.18.0RJ. 

1 EQ.20.OR.I.EQ.22.OR.1.EQ.24.OR.I.EQ.26.OR.I.EQ.28.OR. 

I I.EQ.30.0R.1.EQ.32.0R.I.EQ.34.0R.I.EQ.36.0R.I.EQ.38. 

I OR.I.EQ.40.OR.I.EQ.42.OR.I.EQ.44.OR.I.EQ.46.OR.1.EQ. 

1 48.OR.I.EQ.50X)R.LEQ.52X>R.l.EQ.54.QRJLEQ.56.OR. 

I I.EQ.58.0R.I.EQ.60.0R.I.EQ.62.0R.I.EQ.64.0R.I.EQ.66 
I .ORJ.EQ.68.ORJ.EQ.70.OR.I.EQ.72.ORXEQ.74.OR.I. 

1 EQ.76.OR.I.EQ.78.OR.I.EQ.80.OR.I.EQ,82.OR.I.EQ.84.OR 
1 .I.EQ.86.OR.I.EQ.88.OR.I.EQ.90.OR.I.EQ.92.OR.I. 

1 EQ.94.OR J.EQ.96.0R.I.EQ.9S.ORJ.EQ. I00)THEN 
I LOCX=N SCAN- J !FOR ZIGZAG SCAN, REVERSE BACK FOR DISPLAY 
ELSE 

ILOCX-J-1 

END1F 

eIseif(zigans.eq.'N')then !george wood addition 
I LOCX-J- 1 ! George Wood addition 
ENDIF 
ILOCY=I-I 

IF (IT.EQ. 1 .OR.IT.EQ.2)THEN 

SLANTX=SLANT 1 

SLANTY=SLANT3 

ELSEIF (IT.EQ.3-OR. IT.EQ.4)THEN 

SLANTXR=SLANT2 

SLANTYR=SLANT4 

ENDIF 

IF(NOTH_I.EQ.l)TYPE V ’ 

IF (NOTHJ.EQ.l)TYPE */IT=\IT 
IF (NOTHJ.EQ. 1)TYPE */PLACE=\ PLACE 
IF (NOT! IJ.EQ.l )TYPE *,'SCHEME=’, SCHEME 
IF (NOTH _I.EQ.1)TYPE *' ’ 

IF (IT.EQ.LOR.IT.EQ.2)THEN !CORRECT DELAY WINDOW FOR NONLEVELNESS 
DELAY CORR=DELAY<IT)+(SLANTX*ILOCX*SCANSTEP)+(SLANTY*ILOCY*INDXST£P) 
ELSEIF (IT.EQ.3.0R.IT.EQ.4)THEN 

DELAY_CORR-DELAY(ITMSLANTXR*lLOCX*SCANSTEP)+(SLANTYR*TLOCY*INDXSTEP) 
ENDIF ~ 

TYPE*/ 1 
TYPE*.' • 

TYPE*: ' 

TYPE *:!=■, j; DELAY(',1T;)- \DELAY(IT) 

TYPE */DELAY_CORR= \DELAY_CORR 
TYPE*,’ ' 

TYPE*: ’ 

TYPE*,' ' 


[3 
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IF (PLACE.GE.l.AND.SCHEME.EQ.'B'.AND.NOTH_I.EQ.1)TIlEN 
IIIT=2 'CORRECT DELAY WINDOW FORNONLEVELNESS 

DELAY_CORR-DELAY(2)+(SLANTX*ILOCX*SCANSTEP)+(SLANTY*lLOCY*INDXSTEP) 

CALL PUTTIME(DELAY CORR) ! Set delay FOR B2 

!!! USER-DEFINE VOLTAGE SETTINGS? 

IF (VOLANS.EQ/U')THEN 
VOLTS=VOL(2) 

CALL SETV OLTD1 V (V OLTS) 

ENDIF 

IF (VOLANS.EQ.'A')CALL AUTOSETVOLTS( MAI,VOLT5,NUMINT )!Sct V/D 
GOTO 1132 
ENDIF 

C IKEEP TRACK OF POSITION WHERE DATA IS ACTUALLY TAKEN WITH A COUNTER 
IF (IT.EQ. 1 .AND. PL ACE. GE. 1 . AND. SCHEME. EQ.X)SC AN KNT=SC AN KNT+ 1 
IF (NOT H I . EQ . I )T YPE */ ' 

IF (NOTH I.EQ. 1)TYPE */X=’,J/ Y=\I/ SCANKNT=',SCANKNT 
IF (NOTTI_I.EQ.l)TYPE */ ' 

CALL PUTTIME( DELAY_CORR ) ! Set delay AS USUAL 

100 CONTINUE 

IF ((VOLANS.EQ.'U'-AND.CIT.EQ.l .OR.IT.EQ.2.0R.IT.EQ.3)).OR. 

1 (U0LANS.EQ.TJ'.AND.1T.EQ.4))THEN !!! USER-DEFINE VOLTAGE SETTINGS 
VOLTS=VOL(IT) 

CALL SETVOLTDIV(VOLTS) 

GOTO 1 132 
ENDIF 


CALL AUTOSET VOLTS( MAI,VOLTS,NUMINT ) ! Set V/D 

1132 lF((VOLTS.GT. 1.0).OR. (VOLTS. LT.0.01))THEN 
WRITE(5,I I31)VOLTS 

1 13 1 FORMATC+VBAD VOLTAGE SETTINGS 10.5) 

CALL TEKRESET(MAI) 

GOTO 100 
ENDIF 

CALL GETSA(NAVE,MA1,A) ! Get waveform 
WATCHl = A(l)/ NAVE ! Make sure its acceptable 

TYPE*/ • 

TYPE*/ ' 

TYPE*/ ' 

WRITE(5 ? 1 133) WATCH 1 , WATCH 

1 133 FORMATC+'/CURRENT WATCH LEVEL IS',F10.5,'MIN =‘,F10.5) 
TYPE*/ * 

TYPE*/ ' 

TYPE*/ ' 


IF ( VOLA N S .EQ/U' . OR . UOLAN S. EQ .TJ ' )GOTO 56732 !! SKIP WATCH 


IH 
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[F( WATCH1.LT. WATCH )THEN 
WFLAG = 1 
WRITE( 5,10200) 

FLAG-H’ !H = SAMPLE HOLDER 
TYPE *,’FLAG-,FLAG 
C 

C {COUNTER FOR BAD WATCH LEVEL, IE. ON WHICH WAVEFORM (B1 OR B2) DOES IT OCCUR 
C 

C IF (SCHEME.EQ.'X'.AN D. IT.EQ. 1 ) WATCH_COUNT~l 

C IF (SCI IEME.EQ/X'. AND.IT.EQ.2) WATC1 l_COUNT=2 

C TYPE VWATCH_C0UNT-,WATCH_C0UNT 

GOTO 20000 

1 0201 formate \A/ \A,‘ X= U3,' Y= M3 ) 

1 0200 formatC WAVEFORM BELOW "WATCH", GOING TO NEXT SCAN POINT') 

C GOTO 100 

ENDrF 

Q** ****** *+* + **************#************»********;«***** ****** 

56732 CONTINUE 

£****************************** ***************************** 


IF (NOTH_I.EQ.l)TYPE ' 

IF (NOTH I.EQ. 1 )TYPE * 'PLACE-', PLACE 
IF (NOTH I.EQ.1)TYPE */SCHEME=\SCHEME 
IF (NOTHJ.EQ.l)TYPE */TSCHEME=',TSCHEME 
IF (NOTHJ.EQ.l)TYPE */ZEROO=\ZEROO 
IF (NOTH J.EQ.l)TYPE */IIIT=MllT 
C IF (NOTH I.EQ.1)TYPE */WATCH_COUNT=',WATCH_COUNT 
IF (NOTH I.EQ.l)TYPE */ ' 

C **** EXAMINE B2(t) APPROACH FOR SAMPLE HOLDER ******* C 
IF (PLACE.GE. ) . AND.SCHEME.EQ/B’. AND.IIIT.EQ.2)THEN 
AA “ 0. ! 

DO 499 111= 1 ,5 12 ! Subtract average from wave, 

A1 1 (III)— A(III) 

A1 1(III)=A1 1(1II)/NAVE ! 

499 AA-AA + A1 1(111) ! REALWAVE - Waveform in volts 

ZERO = A A/5 1 2. ! 

IF (NOTHJ.EQ.l)TYPE */ZERO-',ZERO 
FF-0. 

DO 599 111=1,512 

FI (III) = ( REAL(A1 1(111)) - ZERO )*( VOLTS )*10./512. 

599 FF=FF t ABS(F 1 (III)) 

ZEROO-FF/512. 

IF (NOTHI.EQ.l)TYPE ’ 

IF (NOTH I.EQ.1)TYPE */ZEROO=',ZEROO 
IF (NOTHLEQ.l)TYPE 1 

c ******** MOD FOR AUTOMATIC NOISE THRESHOLD ******************* 

C SET UZEROO - 2*(F!RST "NOISE" LEVEL DETECTED ~> ASSUME WE ARE STARTING 
C ON SAMPLE HOLDER) 

I F (TSCHEME.EQ/A'. AND.PLACE.EQ. 1 )THEN 

UZEROO=2*ZEROO 

IF (NOTH_l.EQ.l)TYPE ‘ 

IF (NOTHJ.EQ. 1)TYPE • ' AUTO THRESHOLD VOLTAGE NOISE LEVEL =',UZEROO 
1F(N0THJ.EQ.1)TYPK*/ ' 
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ELSEIF (TSCHEME.EQ.'M'.AND.PLACE.EQ. 1 )THEN 
IF (NOTH_I.EQ. 1 )TYPE ' 

IF (NOTHI.EQ. 1)TYPE MANUALLY-SET THRESHOLD VOLTAGE NOISE LEVEL =',UZEROO 
IF (NOTH _I.EQ. 1)TYPE ' 

ENDIF 

Q *********************************** ***** + ****>|r###****ii(*i|:* 

IF (ZEROO.LT. UZEROO)THEN !!!< EX..002 VOLTS - ON SAMPLE HOLDER 

56120 FLAG-H' !H = SAMPLE HOLDER 

IF (NOTHJ.EQ.l)TYPE *.'FLAG-,FLAG 
GOTO 20000 

ELSEIF (ZEROO . G£ . UZERO 0)T I IEN !!!> EX..002 VOLTS - ON SAMPLE 

56121 FLAG— S' !S - SAMPLE 

IF (NOTHJ.EQ.l)TYPE VFLAG=AFLAG 
SCHEME-'X' 

GOTO 10386 ISTARTTAKEDATA LOOP AGAIN 

ENDIF 

ENDIF 

IF (NOTH I.EQ.l)THEN 

5 00 WRITE( 1 6 ( REC=I 1 6,FMT=53)DELAY_CORR 

116=116+1 

WRITE(16,REC=116 > FM'J -53) VOLTS 
116=116+1 
C 

DO 54321 IJI=1,5I2 
I17SUB=((I 17-1 )*5 L2)+IJI 
54321 WRITE(17,REC=I17SUB,FMT=54) A(IJI) 

117-117+1 

elseif (noth_l.eq.2)then 

5500 WRITE(26,REC=I26,FMT=53)DELAY CORK 
126=126+1 

WRITE(26,REC=I26,FMT=53)VOLTS 

126=126+1 

C 

DO 54329 IJI=1,512 
I27SUB=((I27-1)*512)+IJ1 
54329 WRITE(27,REC=I27SUB,FMT=54) A(IJI) 

127=127+1 


elseif (noth_Leq.3)then 

5800 WRITE(36,REC=I36,FMT=53)DELAY_CORR 
136=136+1 

WRlTE(36 f REC=I36,FMT=53)VOLTS 

136=136+1 

C 

DO 58329 I J I~ 1 , 5 i 2 
I37SUB=((I37- 1 )*5 1 2)-IJ] 

58329 WR1TE(37,REC=137SUB S FMT=54) A(IJI) 
137=137+1 


elseif (nothJ.eq.4)then 


l(t 
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5900 WRITE(46,REC=146.FMT=53)DELAY CORR 
146=146+1 

WRITE(46,REC=I46,FMT* t 53) VOLTS 
146=146+1 
C 

DO 59329 1JI=T.512 
I47SUB=(( 147-1 )*5 12)+IJl 
59329 WRITE(47,REC=I47SUB,FMT=54) A(1JI) 
147=147+1 


cndif 


C 

C WRITE( 8,10600 )A(200) 

£************************************************ ******** **£ 

IF (IT.EQ. 1 )VOLTAGE(l )= VOLTS 
IF ( IT.EQ .2 ) VOLTAGE(2)= V OLTS 
IF (IT.EQ.3)VOLTAGE(3)=V OLTS 
IF (lT.EQ.4)VOLTAGE(4)=VOLTS 
12 **** ********************************** *******************(^ 

600 CONTINUE 

(^* ******** ******************************* ********** 

7 FORMAT (A4) 

£****** ******************************************* ***(2 
20000 IF( WFLAG.EQ.l )THEN 
CLOSE(8) 

OPEN( unit=8, file=TXAl:\ status='OLD' ) 

ENDIF 

IF (SCHEME.EQ.X^SCHEME-B 1 ! RESET BACK TO 'B' FOR NEXT SCAN POINT 
C IF (I.NE.O.AND.J.NE.O)THEN 


C !RESET COUNTERS TO EXCLUDE DATA GATHERED ON POINT 
C WHERE BAD WATCH LEVEL WAS FOUND 


C IF (WATCH_COUNT.EQ. l)THEN !B 1 HAD BAD WATCH LEVEL 

C SCANKNT=SCANKNT-1 

C 116=116-2 

C 117=117-1 

C ELSEIF (WATCH_COUNT.EQ.2)THEN !B2 HAD BAD WATCH LEVE1 , 

C SCANKNT-SCANKNT- 1 

C 116=116-4 

C 117=117-2 

C ENDIF 

if (noth 1 .eq. 1 )WRITE(1 2)JJ,FLAG,SC ANKNT 1NOTE LOCATION WITH FLAG & LEAVE SUBR. 
IF (NOTHJJEQ.l)TYPE 1 
IF (NOTHJ.EQ.l)TYPE *' ' 

IF (NOTH I.EQ. 1)TYPE */ ’ 

IF (NOTH J.EQ. 1 )TYPE *'NOW IN TAKED ATA :','X=', J , ' Y=,!/ FLAG-, FLAG 
IF (NOTH J.EQ. 1 )TYPE * ' ’ 


11 
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IF (NOTH I.EQ. 1 )TYPE 
IF (NOTH_I.EQ. 1 )TYPE * ' 
IF (NOTH_I.EQ. 1 )TYPE 
C ENDIF 

C WATCH_COUNT-0 

RETURN 
end 


C 


it 
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SUBROUTINE PRESSURE! VOLT AGEL,VOLTAGEU) 

C 

C**** PRESSURE=20*VOLTAGE **** (PLOWER.LE. VOLTAGE.GE.PU PPER) 
C 

C Adjust Z axis to get good pressure (PS1) 

C «< 1 1 = DOWN >» 

C 

C PUPPER - .85 ! Upper Pressure =17 

C PLO WER ~ .75 ! Lower Pressure = 1 5 

C PUPPER = .3 ! UPPER PRESSURE - 6 

C PLO WER = .2 ! LOWER PRESSURE = 4 

C PUPPER - .5 ! UPPER PRESSURE = 10 

C PLOWER = ,4 ! LOWER PRESSURE = 8 

C PLOWER = .5 ! LOWER PRESSURE = 10 

C PUPPER = .6 ! UPPER PRESSURE- 12 

C PLOWER = .7 ! LOWER PRESSURE - 14 

C PUPPER - .8 ! UPPER PRESSURE - 16 

C PLOWER - .6 ! LOWER PRESSURE- 12 

C PUPPER - .61 ! UPPER PRESSURE = 12.2 

PLOWER-VOLTAGEL 
PUPPER-VOLTAGEU 
200 CALL GETFLUKE( P ) 

IF( P.GE.PLOWER .AND. P.LE.PUPPER )GOTO 900 
IF( P.LT.PLOWER )CALL MOVXYZ( 3, 1 ) 

IF( P.GT.PUPPER )CALL MOVXYZ{ 3, -I ) 

GOTO 200 
900 RETURN 
end 
C 


n 
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SUBROUTINE PRESSURE MAX(VOITAGE_MAX,LLL,FILENAME t J f U15,CHECKP) 
C 

CHARACTER FILENAME*34 
INTEGERS CHECKP,LLL,LINCR,U,JXXX,115,mT 
52 FORMAT(A2) 

P_MAX=VOLTAGE MAX 
200 CALL GETFLUKE( P ) 

IF( P.GTJP_MAX)THEN 
CHECKP=CHECKP-*-l 
IF (CHECKP.EQ. I )LINCR-LLL/3 
IF (CHECKP. EQ .2 ) LINCR=LLL/2 
IF (CHECKP.EQ. 3)L1NCR= : 1000 
DO 2900 ITTT=1.LINCR 

2900 CALL MOVXY2( 3,-1 ) 13 IS 2-AXIS, MOVE Z-AXIS UP 
IF (Cl I ECKP. EQ.3 )TH EN 

JXXX-1-1 !I JS Y-LOC; GO BACK TO PREVIOUS ROW 
W RITE( 1 5 ,REC— 1 1 5 ,FMT=52) JXXX 
1 15-1 15 * 1 
STOP 
ENDIF 
ENDIF 
900 RETURN 
end 
C 


2Q 
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SUBROUTINE MOVORG( SCANAXIS, NSTEPS. VOLT AGEL. VOLT AGEU, 

I LOOPTHRU_H) 

Run the thing back to the SCANAXIS origin 

integers SCANAXIS 
inleger*2 NSTEPS, LOOPTURUH 

DO 1000 K= 1 , LOOPTHRU II 
1000 CALL MOVXYZ( 3,-1) ' \ Move up 

DO 1200 K=l, NSTEPS- 1 

200 CALL MOVXYZC SCANAXIS, -l ) ! Move back 
DO 1400 K=I,LOOPTHRU_H 
1400 CALLMOVXYZ(3, 1 ) ! Move down 

CALL PRESSURE(VOLTAGEL.VOLTAGEU) hiothick 
CALL MOVXYZ( SCANAXIS, -1 ) ! Move back & forth 
CALL MOVXYZ( SCANAXTS, 1 ) 

RETURN 

end 

SUBROUTINE MOVORGXY( ZIGANS ; SCANAXIS,NSCANJNDXAXIS r NINDX) 

Run the thing back to the SCAN origin 

CHARACTER ZIGANS* 1 
integers SCANAX1S4NDXAX1S 
integer*2 NSCAN.NINDXJ 


IF (ZIGANS.EQ.'Y‘)THEN 
I=Nvndx 

IF (LEQ.LOR.LEQ.3.OR.IJBQ.5.0R.I.EQ.7.OR.I.EQ.9 
1 .OR.l.EQ.l 1.0R.I.EQ.I3.0R.I.EQ.15.0R.J.EQ.17.0R.I. 

I EQ. 1 9.0R.LEQ.2 1 .0R.I.EQ.23.0R.1.EQ.25.0R.I.EQ.27.0R. 

I I .EQ .29 OR.I.EQ .3 LO R.LEQ .3 3 .OR.1.EQ .3 5 .OR.I .EQ .37 . 

1 0RJ.EQ.39.0R.I.EQ.4 1 .ORLEQ.43.OR.LEQ.45.OR.LEQ. 

1 47.ORJ.EQ.49.ORXEQ.5I.OR.LEQ.53.OR.1JEQ.55.OR. 

1 LEQ.57.OR.LEQ.59.OR.LEQ.61.OR.LEQ.63.OR.LEQ.65 
1 .0R.1.EQ.67.0RXEQ.69.0RJ.EQ.71.0RXEQ.73.0RX 
1 EQ.75.ORXEQ.77.ORXEQ.79.OR.LEQ.8LOR.LEQ.83.OR 
1 XEQ.85.ORXEQ.87.ORXEQ.89.ORXEQ.91 .OR.I. 

1 LQ.93.0RJ.EQ.95.0R.IJBQ.97 0R.LEQ.99)TIIEN 


DO 1200 K=I,NSCAN-1 
1200 CALL MOVXYZ( SCANAXIS, -1 ) 1 
END IF 
END1F 

Move back X 

DO 1400 K=1,NINDX-1 
1400 CALL MOVXYZ( INDXAXIS, -1 ) ! 

Move back Y 


<21 
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RETURN 

end 


C 


Si- 
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SUBROUTINE ENTERP A R AM ( SCAN STEP, INDXSTE P,N S C AN ,N INDX,N A VE, VO LT AG EL 
1 ,VOLTAGEU J VOLTAGE_MAX.SHAPEANS,KOOPANSXOOPTHRU_H,LOOPTHRU _S ? 

1 ZIGANS.SCHEME,TSCHEME,UZEROO,FI LENAMEJ 1 5, SCANMODE, SLANT I ,SL ANT2, 
i SLANT3,SLANT4) 

integer*4 SCAN STEP, INDXSTEP 

real SCAN D 1ST, INDXD I ST.PU ,PL,S LANT,rs 1 ant I ,RSL ANT2 , SLANT I ,SL AN T2 

REAL SLANT3,SL ANT4,RSLANT3 .RSLANT4 

JNTEGER*2 LOOPTHRUJLLOOPTHRU S,ISLANTl,ISLANT2 

integer*2 NSCAN,NINDX, IDIST1,IDIST2, 16,115 

integers FREQ,NAVE,ISLANT3, ISLANT4 

character FILENAME*34, FTLEEXT*23, C HEADER* 32 

CHARACTER LOOP ANS*1, SCHEME* I 

CHARACTER SHAPEANSM RIGANS* I 

CHARACTER TSCHEME* 1 ,SCANMODE* 1 

DATA 16/1/ 


115=1 

10010 format(A) 
10012 format(A2) 
10020 format( I ) 
10030 foi*mat(F) 


Initialization and entry of parameters 


TYPE * ' ' 

TYPE * ' ***+*> DONSCAN.FOR AUTO SCAN’ 

TYPE*, 1 ‘ 

OPEN (UNIT=10,FILE= r NOTHlCK_ALLSH APE 1. DAT', STATUS-OLD’, 
t FORM-’FORMATTED') 

READ (10,93764) SCANMODE 

READ (10,93765) XSCANSTEP !(X 1ND) 

IF(SCANMODE.EQ.'P’)XSCANSTEP=2000. 

READ (10,93765) YINDXSTEP !(Y 1ND) 

IF (SCANMODE.EQ.’P')YINDXSTEP-2000. 
SCANSTEP=JNINT(XSCANSTEP) 

INDXSTEP-JNINT(YINDXSTEP) 

READ (10,93765) SCAND1ST !(X DIST) 

IF (SCANMODE.EQ.'P')SCANDIST=2000. 

READ (10,93765) INDXDIST !(Y DIST) 

IF (SCANMODE.EQ.'P^INDXDIST^OOO, 

IF (SCANMODE.EQ.'L/)THEN 
YINDXSTEP=20000 . 

1NDXSTEP=JN1NT(Y INDXSTEP) 

INDXDIST=20000. 

ENDIF 

READ (10,93764) ZIGANS 
READ (10,93765) PL 
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READ (10,93765) PU 
READ (1 0,93765) PM 
READ (10,93764) SHAPE ANS 
READ (10,93764) LOOPANS 
READ (10,93766) LOOPTHRUH 
READ (10,93766) LOOPTHRU_S 
READ (1 0,93764) TSCHEME 
READ (10,93765) UZEROO 
READ (10,93765) TH1CKN 

************** "SHAN'T” IS THE NUMBER OF NANOSECONDS / MICRONS 
THAT THE ECHOES MOVE DUE TO NONLEVELNESS OF 
EXPERIMENT OR THICKNESS VARIATION. IT MUST BE 
INCORPORATED INTO THE PROGRAM AS AN ADDITIVE OR 
SUBTRACTIVE VALUE TO ALL THE DELAY TIMES BY FIRST 
MULTIPLYING ’’SLANT" BY THE CUMULATIVE INDEXING DISTANCE 
IN THE X-D1RECTION (N*XSCANSTEP) AS THE SCAN IS 
BEING PERFORMED AND THEN ADDING / SUBTRACTING THIS 
VALUE TO EACH DELAY. "SLANT" MUST BE DETERMINED 
BEFORE THE SCAN IS RUN BY MOVING THE X AND Y STAGES 
OVER THE SCAN LENGTH AND NOTING THE SHIFT TIME OF 
AN ECHO, SAY B I . THIS SHIFT WILL BF. CONSISTENT FOR 
B1 AND B2 ECHOES IN MOST CASES BUT ANOTHER SLANT 
FACTOR WI LL HAVE TO BE DEFINED FOR THE ECHO OFF THE 
REFLECTOR PLATE. SO WE HAVE SLANT 1 AND SLANT 2. 

READ (10,93765) SLANT l 
SLANT1=SLANT1/10.**9. 

’CONVERT SLANT TO INTEGER FOR STORAGE IN 
.DATI2 FILE 


{CONVERT SLANT TO INTEGER FOR STORAGE IN 
.DATI2 FILE 


{CONVERT SLANT TO INTEGER FOR STORAGE IN 
.DATI2 FILE 


{CONVERT SI .ANT TO INTEGER FOR STORAGE IN 
.DATI2 FILE 


RSLANT 1 =SLANT 1 *10**13. 
C 

I S L ANT 1 =1 INT(R S LAN T I ) 

READ (10,93765) SLANT2 
SLANT2=SLANT2/I0.**9. 
RSLANT2=SLANT2*10.**13. 
C 

ISLANT2=UNT(RSLANT2) 

READ (10,93765) SLANT3 
SLANT3=SLANT3/1 0.**9. 
RSLANT3=SLANT3*10.**13. 
C 

1SLANT3~] INT(RSLANT3) 

READ (10,93765) SLANT4 
SLANT 4=S L AN T4/ 1 0 . * *9 . 
RSLANT4=SLANT4* 1 0.** 13. 
C 

I SLANT4H 1 NT(RSLAN T4) 


READ (10,93765) FREQUENCY 
READ (10,93764) FILEEXT 


5V 


126 


K=0 

K=K+1 
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IF( FILEEXT(K;K).NE.' ’ )GOTO 126 ! 

FILEEXT(K+1 :fC+4)=’.DAT' ! UNIT 6 - raw data 

READ (10,93764) CH FADER 

93764 FORMAT (A) 

93765 FORMAT (F) 

93766 FORMAT (1) 

C VOLTAGEL-PL/20. 

C VOLTAGEU-PU/20. 

C VOLTAGE_MAX=PM/20. 

C TYPE *,’RUN TEMPDIRSCANCOM - DTRSCAN SMALLTRANS_DATA.DAT 1 

TYPE */DATA FOR ULTRASONIC SCAN (ALL FLOATING PT. EXCEPT INDEXES)’ 

TYPE */SCAN DIMENSIONS/INDEXING DATA (MICRONS)’ 

TYPE */X INDEX - INTEGER EXPRESSION 1 
TYPE *,SCANSTEP 

TYPE */Y INDEX - INTEGER EXPRESSION’ 

TYPE MNDXSTEP 

TYPE *; X SCAN DISTANCE (MICRONS) - FLOATING POINT 1 
TYPE *,SCANDIST 

TYPE X SCAN DISTANCE (MICRONS) - FLOATING POINT’ 

TYPE *, INDXDIST 

C TYPE LOWER & UPPER PRESSURES FOR CONTACT SCAN’ 

C TYPE *,PL,PU 

C TYPE ^THICKNESS (MM)’ 

C TYPE *,THICKN 

TYPE *, ’SLANT 1 (X-LEVELNESS TIME CORRECTION FACTOR for bl,b2) (nsec / urn )’ 

TYPE *, SLANT! 

TYPE *,'SLANT2 (X-LEVELNESS TIME CORRECTION FACTOR for reflector echoes) (nsec / um )’ 
TYPE *,SLANT2 

TYPE VSLANT3 (Y-LEVELNESS TIME CORRECTION FACTOR for bl ,b2) (nsec / um )’ 

TYPE *,SLANT3 

TYPE VSLANT4 (Y-LEVELNESS TIME CORRECTION FACTOR for reflector echoes) (nsec / um )’ 
TYPE *,SLANT4 

TYPE ^TRANSDUCER CENTER FREQUENCY (FLOATING POINT)’ 

TYPE FREQUENCY 

TYPE VFILENAME FOR RAW DATA ’ 

TYPE *,FILEEXT 

TYPE VCHEADER INFO (UP TP 32 CHARACTERS)’ 

TYPE *,CHEADER 
TYPE ' 

1000 NSCAN = (SCANDTST/SCANSTEP) + 1 
NINDX = (INDXDIST/1NDXSTEP) - 1 

SCAN D 1ST = SCAN STEP* (NSC AN -1) ! SCAND1ST - negative 

INDXDIST = INDXSTEP*(NTNDX-1) 

IF (SCANMODE. EQ.’S')TUEN 

WRITE( 5,1 1 100 )SCANDIST, NSCAN, SCANSTEP, 

+ INDXDIST, NINDX, INDXSTEP 
ELSEIF (SCANMODE.EQ. , L')rHEN 
WRITE( 5,1 1200 )SCAND1ST, NSCAN, SCANSTEP 
ELSEIF (SCANMODE.EQ.'P’)THEN 
W'RITE( 5.11300 ) 

ENDIF 
11100 fbrmat( /' 

+ IP X scan to ',F7.0,\ ’,13,' steps of ’,15, 

+ IP Y index to \F7.0,', ',13,' steps of ’,15) 
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11200 l'ormat(/", 

- IP Line scan to \F7.0,\ \I3,’ steps of \I5) 

1 1 300 format( P \ IP Point Measure repeated 2 times’) 
SCIIEME-’B’ 


C 


11690 


THICK= IINT(THICKN * 1 000.) 

FREQ=FREQUENCY 

NAVE-64 

FILENAME = ’[R OTH.D A T A] 7/F 1 LEEXT 
WRTTE( 5,11690) 

fonnat( /'SNUMBER OF AVERAGES SET AT 64' ) 


DIRECT ACCESS STORAGE 


OPEN( unit-6, file-FILENAMEZ/'CIF, status-NEW', 

ACCESS= , DIRECT , ,RECORDTYPE= l FIXED’ ; RECL-32, 

-r form- FORMATTED\ORGANIZAT!ON- SEQUENTIAL' ) 

OPEN( unit- 15, file- FILENAME// 1 12', status-’NEW', 

+ ACCESS- DIRECT'.RECORDTYPE- FIXED’, RECL-2, 

+ form-'FORMATTED', ORGANIZATION- SEQUENTIAL’ ) 

OPEN( unit- 16, file=FILENAME//’R4\ status-TMEW\ 

+ ACCESS- D I RECT,RECORDTYPE- FIXED’, RECL-4, 

H form- FORMATTED\ORGANIZATION- SEQUENTIAL' ) 

OPEN( unit-17, fi Ic-FILEN AME// 1 W A V', slatus-’NEW', 

+ ACCESS- DIRECT 1 ,RECORDT YPE-T1XED’ ,RECL=2, 

+ form- FORMATTED’, ORG AN IZA'flON- SEQUENTIAL’ ) 

OPEN( unit-26, fi 1 c=F 1LEN AME// , RR4', status-’NEW\ 

I ACCESS- DIRECT^RECORDTYPE- FIXED’, RECL-4, 

+ form=’FORMA'ITED\ORGANIZATION-’SEQUENTIAU) 

OPEN( unit-27, fiIe=FlLENAME//’RWAV’, status-’NEW, 

+ AC CESS- DIRECTLRECOR DT YPn-TI XED', RECL-2 , 

+ form- FO R M ATTEDLORG AN IZATION- SEQU EN TI AL' ) 

OPEN( unit-36, file=FILENAME//’SR4', status-NEW, 

+ ACCESS - DIRECT*, RECORDT Y PE- F I X ED', REC L-4, 

1- form-’FORMATTED', ORGANIZATION-’SEQUENTIAL') 

OPBN( unit-37, file-FILENAME/7'SWAV’, stalus-’NHW', 

+ ACCESS-'DIRECT’.RECORDTY PE- FIXED', RECL-2, 

+ form -'FORMATTED’,ORGANIZATION- SEQUENTIAL’ ) 

OPEN( unit-46, file=FILENAME//TR4’, status-’NEW', 

+ ACC ES S- DIRECT, RECORDT Y PE- FI XED’, RECL-4, 

+ form- FORMATTED', ORGANIZATION- SEQUENTIAL’ ) 

OPEN( unit-47, file=FILENAME//TWAV’, status^NEW’, 

+. ACCESS- DIRECT 1 , RECORDTYPE _, FIXED', RECL-2, 

+ form=’FORMATTED’, ORGANIZATION- SEQUENTIAL’) 

51 FORMAT(A32) 

52 FORMAT(A2) 

53 FORMAT(A4) 

54 FORMAT(A2) 

C 

C Write header : X, Y dist, X, Y scan points, dummies 

C 

WRTTE(6,REC=16,FMT=5] )CHEADER 
16=16+1 

WRITE(6,REC=I6,FMT-5 1 )SI IAPEANS 


3C 



non 


69 


5,629,865 


70 


16-16+1 

WRITE(6,REC=I6,FMT=5 1 )ZIGANS 
16=16+1 

IDIST1 - - HF1X(SCANDIST/1000.) 

IDIS'17 = I1FLX(INDXDIST/1000.) 

IF (IDIST2.LT. 1 )1DIST2=1 1TRICK TO ALLOW DISPLAY OF IMAGE ON PSIDD 
WRITE(15,REO=ll5,FMT=52 )IDIST1 
115=115+1 

WRITE(1 5,REC=I 1 5,FMT=52)1DIST2 
115=115+1 

WRITE(1 5,REC=n 5,FMT=52)NSCAN 
115=115+1 

IF (SCANMODE.EQ. , L , .OR.$CANMODE.EQ.’P')NINDX=! !MOD FOR LINE/POINT SCAN 

WRJTE( 1 5,REC=I 1 5,FMT=52)NINDX 

115=115+1 

C W RITE( 1 5,REC=1 1 5.FMT=52)T111CK 

WR1TE( 1 5,REC=1 1 5,FMT=52)ISLANT 1 
115=115+1 

WRITE( 1 5 ,REC=1 1 5,FMT=52)i SLANT2 
115=115+1 

WRITE(15,REC~r 1 5,FMT=52)ISLANT3 
I15=ri5+I 

WR ITE( 1 5 ,REC~1 1 5 , FMT=52)ISL ANT4 
115=115+1 

WRJTE( 1 5,REC=I 1 5,FMT=52)DENS 
115=115+1 

WRITE( 1 5,REC=I 15 ,FMT=52)FREQ 
115=115+1 

WRITE( 15,REC=I 1 5,FMT=52)N A VE 

115=115+1 

RETURN 

end 

SUBROUTINE XPARAM( SCANSTEP,INDXSTEP,NSCAN,NINDX ) 
inleger*4 SCAN STEP, JNDXSTEP 
integer*2 NSCAN.NINDX 
SCANSTEP = 50 
INDXSTEP — 50 
NSCAN =50 
NINDX =50 
THICK =1 

OPEN( unit=6. file^XXX.DAT*, status=’OLD\ 

+ fonn- UNFORM ATTED 1 ) 

SCANSTEP = 50 
RETURN 


SUBROUTINE MOVKLING (AXIS,DIST) 

MOVE KLINGER A SPECIFIED DISTANCE IN SPECIFIED DIRECTION 

byte IXYZ 
character* 1 XYZ 
equivalence( XYZ, IXYZ ) 
common WRK0,BUFFER 


37 
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CALL STRTGP1B 
CALL INITINSTR(2) 

CALL INITINSTR(j) 

CALL INITINSTR(4) 

10010 format(A) 

10020 format(I) 

10100 format( /'SDi stance ?’ ) 

100 WR1TE( 5,10000 ) 

10000 format( /'$X, Y, or Z ?' ) 

READ( 5, 1 00 1 0 )XYZ 

1F( 1XYZ.GT.90. OR .IXYZ.LT.8S )GOTO 900 

IPORN=l 

WR1TE( 5,10100) 

READ( 5,10020 )IREQ 
1F( IREQ.LT.O )THEN 
IP0RN=-1 
1REQ=ABS(1REQ) 

END IF 

CALL SETXYZ( 1REQ,IREQ,IREQ ) 

I AXIS = IXYZ - 87 

CALL MOVXYZ( IAXIS,1P0RN ) 

GOTO 100 
900 STOP 

RETURN 

end 

INCLUDE ’[ROTRMENUJBASEOJSECJvl.FOR 1 
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C nothick cmnch.for — 4 SCAN METHOD 
C Crunch program for Nothickness velocity scan 
C Don Roth 20-scp- 1 994 


integer*2 SCANDIST,INDXDIST,NSCAN,N1NDX,JJJI,SCANKNT 

integers DENS JRFREQ,AVES,NPOrNT,I6,I15,I 1 6,1 1 8,126 

integer*2 ARRAYSIZE,KC0UNT,SELFREQ,ENDFREQ,IS5 ? 136,146 

INTEGERS islanti,islant2,]SLANT3,ISLANT4,ACTPOINT,JXXX 

intcger*2 rawwavebl (5 12),rauwaveb2(512),rawwavers(512),rawwaverns(512) 

INTEGERM n7,I7,REC 1 I27,17SUBJ37,I47J27SUB,I37SUBJ47SUB,I17SUB 

real *4 L_VEL,U_VEL,AV VEL 

realM ASPECB 1(1 0241, ASPECB2(1024) 

real *4 B1ASPEC(256), B2ASPEC(256) 

real *4 ASPECRS(1024) } ASPECRNS(1024) 

real *4 RSASPEC(256), RNSASPEC(256) 

REALM B1PJHASE(256),PHASEB 1(1024) 

REALM B2PHASE(256),PH ASEB2( 1 024) 

REALM RSPHASE(256),PHASERS( 1 024) 

REALM RNSPHASE(256),PH ASERN S( 1 024) 

REALM BB 1 (5 12),BB2(5 12),RRS(5 12),RRNS(5 12) 

REALM VEL(5000),TRDIAM,BUFLENGTH 

character FILENAMEM2, FILEEXTM6, CHEADER*32 ? N01SANS*1 

character FS1 EXT* 1 6,CALNAME*32,FLAG* 1 

CHARACTER D1FFANS* LPHASE* 1 , PHASE 1 * 1 

CHARACTER SELFREQANSM 

CHARACTER SH APEANS* l ,ZIG ANS* 1 ,DIR*5 

CHARACTER CRUNCH_C0DE*2,CF1LTER*1 

character actual selfreqans* 1 

byte FILEFS1IN 

equivalence ( B1ASPEC,ASPECB1 ) ! ! 

equivalence ( B2ASPEQASPECB2 ) ! +ASPEC 256 ! 

equivalence ( RSASPEC,ASPECRS ) ! ! 

equivalence ( RNSASPEC, ASPECRNS ) f *ASPEC 256 ! 

EQUIVALENCE (BlPHAS^PHASEBl) ! f 

EQUIVALENCE (B2PHASE,PHASEB2) 

EQUIVALENCE (RSPHASE,PHASERS) ! 1 

EQUIVALENCE (RNSPHASE,PHASERNS) 
equivalence ( FS1EXT, FILEFSHN ) 

DATA 16/1/1 1 5/1/1 16/1 /1 1 7/1 /I7/1/11 8/1/ 

DATA I26/1/I27/1 /I85/1/136/I/137/1/I46/ 1 /I47/1/ 

NTPI=0 


10010 foimat(A) 
10014 formatp \ A) 
1 0020 format(I) 
10030 format(F) 

C 


X) 
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OPEN (UNTT-14,FlLE“’NOTHICK_SELFREQl . DAT, 

1 STATUS- OLD', FORM-’FORMATTED') 

READ (14,98652) FILEEXT 
READ (14,98652) DRIVE 
READ (14,98652) PHASE 
READ (14,98652) PHASE1 
READ (14,98652) DIFFANS 
READ (14,10030) BUFLENGTH 
READ (14,10030) BUFVEL 
READ (14,10030) TRDIAM 
CLOSE (14) 

98652 FORMAT (A) 

C ********** OPEN FILTER FILE FOR UPPER & LOWER LIMITS ********** C 


OPEN (UNIT=24,FILE= 1 N0TH1CKJLIM1TSLDAT', 

1 ST ATUS- OLD\FORM- FORM ATTED’) 

READ (24,98652) CF1LTER 
READ (24,10030) L_VEL 
READ (24,10030) U_VEL 
LVEL=L_VEL 
U_VEL-U_VEL 
A V__V EL=(L_VEL+U_ V EL)/2 . 

£ ************* ** *************** ********************* ^ 


K=0 

126 K-K+L 

IF(FlLEEXT(K:K).NE. M )GOTO 126 ! 

FI LEEXT(K> 1 : K+4)- .DAT' ! UNIT 6 = raw data 

IF (DRIVE.EQ.'A’)DIR-DUC2:’ 

IF (DRrVE.EQ. , C , )DIR='DUC0: r 


FILENAME - ’[ROTH.D ATA] V/FI LEEXT 
£********** * * **************** ** ********************************* 
DIRECT ACCESS STORAGE 

OPEN( unit=76. ftle=DIR//FILENAME//'GB , ? status=’NEW\ 

+ ACCESS- DIRECT, RECORDTYPE-FIXED',RECL=4, 

+ form- FORMATTED', ORGANIZATION- SEQUENTIAL 1 ) 

OPEN( unit-66, file-FlLENAME/yCH*, status-'OLD', 

+ ACCESS- DIRECT,RECORDTYPE= = 'FIXED\RECL=32, 

+ form -TORMATHi ^ORGANIZATION- SEQUENTIAL 1 ) 

OPEN( unit- 15, fiIe=FILENAME// , 12', status-'OLD 1 , 

+ ACCESS=’DIRECT\RECORDTYPE-FIXED\RECL=2, 

+ form- FOR MATTED', ORGANIZATION- SEQUENTIAL' ) 

OPEN(unit=16, file-FI LEN AME// 1 R4 status-'OLD', 

- ACCESS- DIRECT',RECORDTYPE-FIXED’,RECL=4, 

+ form-'FORMATTED 1 , ORGANIZATION- SEQUENTIAL’ ) 

OPEN( unit-17, ftle=FlLENAME/"WAV', status='OLD'. 

+ ACCESS- DIRECT, RECORDTYPE- FIXED’, RECL=2, 
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+ form=’FORMATTED',ORG ANIZATION-'SEQUENTl AL’ ) 

OPEN( unit-26, fIle=FILENAME//’RR4\ status= T OLD’, 

+ ACCESS- DIRECT', RECORDTY PE- FIXED' ,,RECL^4. 

4 form- FORMATTED' } ORGANIZATlON- SEQUENT! A L ) 

OPEN( unit-27, ffle=FI LENAME/./’RWA V', status-'OLD', 

+ ACCESS- DIRECT',RECORDTYPE- FIXED 1 , RECL==2, 

+ form- FORMATTED',ORGANIZATl ON-SEQUENTIAL' ) 

OPEN( unit-36. fiIe=FILENAM£/./ , SR4 r , status-'OLD', 

+ ACCESS-DIRECT,RECORDTYPE-'FIXED',RECL-4, 

-I- form- FORMATTED’, ORGANIZATION- SEQUENTIAL' ) 

OPEN( unit=37, file=FILENAME//’SWAV\ status-'OLD’, 

+ ACCESS- 1 DIRECT,RECORDTYPE-FlXED\RECL-2, 

+ form- FORMATTED' .ORGANIZATION- SEQUENTIAL ) 

OPEN( unit=46, file=FILENAME//'TR4Vstatus=’OLD', 

+ ACCESS- DIRECT’, RECORDTYPE-FlXED\RECL=4, 

+ form- FORMATTED', ORGANIZATION- SEQUENTIAL' ) 

OPEN( unit-47, file-FILENAME//'TWAV\ status-'OLD'. 

+ ACCESS=’DIRECT , J RECORDTYPE= , FIXED\RECL-2, 

4 form- FORMATTED 1 , ORGANIZATION- SEQUENTIAL ) 

C 

51 FORMAT(A32) 

52 FORMAT(A2) 

53 FORMAT(A4) 

153 F0RMAT{IX,I4,1X.14JX,14,1X,E12.5,1X,E12.5,1X,E12.5,IX,E12.5,1X, 

+ E12.5,IX,E12.5) 

54 FORMAT(A2) 

K-0 

120 K-K4-1 

IF(FILEEXT(K:K).NE.Y)GOTO 120 ! 

FILEEXT(K-1 :R4'3)- SPC' ( UNIT 6 - raw data 

FILENAME = ‘{ROTH.DATAJV/FILEEXT ! 

C 

C***+*******™*njNn^ 9 S?C p ILE for SPECTRA, vel, FOR DIRECT ACCESS* *******0 

C 

OPEN( unit=7, file=Dl R7/FIL3N AME, slatus-NEW^ACCESS^DIRECT’, 

+ RECORDTYPE- FIXED', RECL=4,FORVl=TORMATTED\ 

4- ORG AN I ZATI ON - SEQUENTIAL) 

WRITE( 6,10200 )FILENAME \ UNIT S = calibration 

1 0200 format( ' Filename of analyzed data is \A ) 

FILEEXT(K+I :K4“3)- CAL' 

CALNAME - '[ROTH.DATA]7/FILEEXT 

OPEN( unit=8, file-CALNAME, status=‘NEW\ ACCESS='DIRECT, 

+ RECORDTYPE='FIXED\RECL=4,FORM~FORMATTED', 

4- ORGANIZATION- SEQUENTIAL' ) 

89651 DO 76065 IOP1- 1,1024 
PliASEBl(lOPIH). 

PH A SEB2(I OP I)=0 . 

PHASERS(IOP1)=0. 

PHASERNS(IOPI)=0. 

ASPECB I{IOP1)~0. 

ASPECB2(TOPl)~Q. 

ASPECRS(IOPI)-0. 

76065 AS PECRN S(IOPt)-0 . 
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K.COUNT=l 

#### Read preliminary scan info # # // § 

79193 READ(66,REC=I6,FMT=5 1)CHEADER 
16=16+1 

READ(66,REC=16,FMT=5 1)SHAPEANS 
16=16+1 

READ(66,REC=I6,FMT=5 l)ZIGANS 
16=16+1 

IF (SHAPEANS.NE.’Y')GOTO 89714 

OPEN(UNIT= 1 2 > FILE- fROTH.DAT A]'//FILEEXT( 1 :K)// I DATHS\ 

1 STATUS- OLD', FORM- UNFORMATTED') !INFO ON WHETHER WE ARE ON HOLDER/SAMPLE 

89714 READ(15,REC=I15,FMT=52)SCAND1ST 
115=115+1 

READ(15,REC=I15,FMT=52)INDXDIST 

115=115+1 

READ( 1 5,REC=I 1 5.FMT=52)NSCAN 
115=115+1 

REA D( 1 5,REC=f 1 5 ,FMT=52)NINDX 
115=115+1 

READ(15 1 REC=I15.FMT=52)islantl 

115=115+1 

READ( 1 5,REC=I 1 5,FMT=52)islant2 
115=115+1 

READ(1 5,REC=1 1 5,FMT=52)isfant3 
115=115+1 

READ( 1 5,REC=11 5,FMT=52)islant4 
115=115+1 

READ( 1 5,REC=I 1 5, FMT=52)DENS 
115=115+1 

RE AD( 1 5,REC=I 1 5,FMT=52)TRFREQ 
115=115+1 

READ( 1 5,REC=11 5,FMT=52)AVES 
115=115+1 

READ( 1 5,REC=l 1 5,FMT=52)JXXX 
115=115+1 

ARRAYSIZE=(NSCAN)*(JXXX) !IF SCAN FINISHED PROPERLY, JXXX=NINDX 

RE AD(1 6.REOI1 6,FMT=53)TIMES ET 
116=116+1 

ENDFREQ=TRFREQ*2.5 

CALLPOINT„FROM_FREQ(ENDFREQ,TIMESET,ACTPOlNT) 


DELTAF = l./(TIMESET*20.) 

NP01NT=NSCAN*JXXX !IF SCAN FfMlSHED PROPERLY, JXXX=NINDX 
N WAVES = 3*NPOfNT 


C ! U Loop thru scan data 
1TEM_COUNT=0 


3-T. 
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9965 1 DO I S90 NIP- 1 ,NPOINT ! INPOrNT IS TOTAL OF ALL POINTS 
SCANKNT-NIP 

ITEM_COUNT=ITEM_COUNT+ 1 

IF (SHAPEANS.EQ.'N'-AND.NIP.GE.2)ITEM_COlJ r NT=2 
TYPE *7 ’ 

TYPE V ITEM COUNT* MTEM COUNT 
TYPE CRUNCH CODE= ',CRUNCH_CODE 
TYPE * ’ * 

CCVACCUM-O. 

67519 CONTINUE 

N IPJ-((NIP- 1 )/N SCAN)+1 1Y LOCATION 

NIPt-NlPI+1 IX LOCATION 

IF (NIPI.EQ.(NSCAN+1))NIPI=1 ! RESET X LOCATION 
C ***** MODIFICATION FOR ZIGZAG SCAN ******************* ***********q 

if(zigans.eq;y’)then 

IF (NIPJ.EQ.2.0R.NIPJ.EQ A .OR.N IPJ .EQ.6.0R.NIP J .EQ.8 .OR.N I PJ.EQ. 1 0 
1 .OR.NIPJ.EQ. 1 2.0R.NIPJ.EQ. 1 4.0R.NIPJ.EQ. 16.OR.NIPJ.EQ.! 8.0R.NIPJ. 

I EQ.20,OR.NIP/.EQ.22.OR.NlPJ-EQ.24.OR 1 NlPJ.EQ.26.OR.NIPJ.EQ.28.OR. 

1 NIPJ.EQ.30.OR.NIPJ.EQ.32.OR.NIPJ.EQ.34.QR.NIPJ.EQ.36.OR.NIPJ.EQ.38. 

1 OR.NIPJ.EQ-40.OR.NIPJ.EQ.42.OR.NIPJ.EQ.44.OR.N1PJ.EQ.46.OR.N1PJ.EQ. 

1 48.OR.NIPJ.EQ.50.OR.NIPTEQ.52.OR.NrPJ.EQ.54.OR.NIPJ.EQ.56.OR. 

1 NIPJ.EQ.58.OR.NIPJ.EQ.60.OR.NIPJ.EQ.62.OR.NIPJ.EQ.64.OR.NIPJ.EQ.66 
I .OR.NIPJ.EQ.68.OR.NIPJ.EQ.70.OR.NTPJ.EQ.72.OR.NIPJ.EQ.74.OR.NIPJ. 

1 EQ.76.OR.NIPJ.EQ.78.OR.NIPTEQ.SO.OR.NIPJ.EQ.82.OR.NIPJ.EQ.84.OR 
1 .N1PJ.EQ.86.OR.NIPJ.EQ.88.OR.NIPTEQ.90.OR.NIPJ.EQ.92.OR.NIPJ. 

\ EQ.94. OR.NIPJ.EQ. 96. OR.NIPJ.EQ. 98. OR.NIPJ.EQ.iOOJTHEN 
NIPG=NSCAN-NIPIf t !FOR ZIGZAG SCAN, REVERSE BACK FOR DISPLAY 
ELSE 

NIPG-NIPI 

ENDIF 

elseifCzigans.eq.N^then Igeorge wood addition 
nipg-nipi ! George Wood addition 
ENDIF 

TYPE ’ 

TYPE * 

IF (SHAPEANS.EQ.'N’)GOTO 89715 

READ(12)JJ,II,FLAG,SCANKNT 'READ WHETHER ’H’ OR 'S' 

TYPE * ‘X^NIPG,’ Y=',N1PJ7 FLAG-', FLAG 


IF (FLAG.EQ.'S')THEN 

897 1 5 READ( 1 6,R EC=I 1 6,FMT=53)DELA YB 1 
116=116+1 

READ(1 6,REC=1 1 6.FMT-53) VOLTSETB 1 
116=116+1 

READ(26,REC=126,FMT=53)DELAYB2 

126=126+1 

READ(26,REC-126,FMT= s 53)VOLTSETB2 

126=126+1 
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RE AD(36 ,REC=I 3 6,FMT=53)DEL AY R S 
136=136+1 

READ(36,REC=I36.FMT=53)VOLTSETRS 

136 = 136+1 

READ(46,REC=I46,FMT=53)DELAYRNS 

146=146+1 

READ(46,REC=I46,FMT=53)VOLTSETRNS 

146=146+1 

DO 54322 IJI=l t 512 

54322 RE AD( 1 7,REC=((I 1 7- 1 )*5 12)+TJ1 ,FMT=54)R A WWA VEB 1 (III) 
117=117+1 

DO 54323 IJI=1,5I2 

54323 READ(27.REC=((I27-1)*5I2)+IJI.FMT=54)RAWWAVEB2(1JI) 
127=127+1 

DO 54324 IJ 1=1, 5 12 

54324 READ(37,REC= s ((I37-i)*512)+IJI,FMT-54)RAWWAVERS(lJI) 
137=137+1 

DO 54325 IJI=I»512 

54325 READ(47,REC= s ((I47-l)*512)+lJI,FMT-54)RAWWAVERNS(IJI) 
147=147+1 

FXSE1F (FLAO.EQ.'H')THEN 

GOTO 1890 

ENDIF 


11 = 11+1 
A = 0. 

DO MOO 1=1,512 

RAWWAVEB 1 (I)=RAWW AVEB l(l)/AVES 

1400 A=A + RAWWAVEB 1(1) 

ZERO = A/512. 

DO 1560 1=1,512 

1560 BB !(!) = ( REAL(R AW WAVEB1 (1» - ZERO )*( VOLTSETB 1 )*10./5 12. 

II=IJ+1 
A = 0. 

DO 1401 1=1,512 

RAWWAVEB2(T)=RAWWAVEB2(I)/AVES 

1 40 1 A=A + RAW WA VEB2(I) 

ZERO « A/512. 

DO 1570 1=1,512 

1570 BB2(I) = (REAL(RAWWAVEB2(1)) - ZERO )*( VOLTSETB2 )*10./512. 


11 = 11+1 
A — 0. 

DO 1402 1=1,512 

RAW WAVERS (I)=RA W WAVER S(I)/A VES 
1 402 A=A + RA W WAVERS(l) 

ZERO = A/5 12. 

DO 15801=1,512 

1580 RRS(I) =( REAL(RAWWAVERS(1»- ZERO )♦( VOLTSETRS )*10./512. 
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A = 0. 

DO 1403 1-1,512 

RAWWAVERNS(T)=RAWWAVERNS(I)/AVES 
1403 A-A + RAWWAVERNS(I) 

ZERO = A/512. 

DO 15901-1.512 

1590 RRNS(I) - ( REAL(RAWWAVERNS(1» - ZERO )*( VOLTSETRNS )*10./5I2. 

IF (PHASE.EQ.'N’)THEN 

CALL CORR(BB 1 ,BB2,TWOTAU_DELA Y) ! for time delay in sample 
ELSEIF (PHASE.EQ/Y')THEN iphase inversion FOR BETWEEN B1 &B2 
CALL MCORR(BB 1 ,BB2,TWOTAU_DELAY) !for time delay in sample 
END1F 

IF (PHASEl.£Q. t N‘)THEN 

CALL CORR(RRS,RRNS,DELTAT_DELAY) !for delay between reflector peaks 
C (sample vs. no sample) 

ELSEIF (PHASE I .EQ.'Y')THEN IPHASE INVERSION FORBETWEEN RN & RNS 
CALL MCORR{RRS,RRNS ; DELTAT_DHLAY) !for delay between reflector peaks 
ENDTF 


CALL OBT AIN_MAGN 1TUDE__SPECTRA(BB 1 ,ASPECB 1 ) 
CALL OBTAIN MAGNITUDE_SPECTRA(BB2,ASPECB2) 
CALL OBT A IN MAGN1TUDE_SPECTRA(RRS, A SPECRS) 

callobtainImagnitude_spectra(rrns,aspecrns) 

ADELAYS = DELAYB2 - DELAYB 1 
ADELAYR - DELAYRNS - DELAYRS 


TWOTAU=ADELAYS+(TIMESE I75L2)*TWOTAU_DELAY ’for time delay in sample 
DELTAT-ADELA YR+(TIMESET/5 1 .2)*DELTATJDELAY !for delay between 
C reflector peaKS (sample vs. no sample) 


CALL NOTH I CK_V EL_CALC(DELT AT,T WOTAU,VEL_NOTHICK) 
type *,'vetynothick= ',vel_nolhiek 


CALL FILTER(VEL NOTHICK, VELPREV.CFTLTER,SHAP£ANS,L VEL,U VEL, 
1 AV_VEL,NIPG ,N IP J , SC A N KN T . 1 1 8) 

VEL(SCANKNT)~VEL_NOTHICK 


VELPREV=VEL_NOTIIICK 

CALL STORE _SPECTRA(ASPECB 1 , ACTPOINT,i7,I7SUB) ’STORE SPECTRAS AT EACH SCAN 

POINT 

CALL STORE SPECTRA(ASPECB2,ACTPOINT,I7,l7SUB) 

CALL STORE_SPECTRA(ASPECRS,ACTPOTNT,I7,I7SUB) 

CALL STORE'" SPECTRA(ASPECRNS,ACTPOINT,I7,I7SUB) 
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1890 CONTINUE 

IF (SHAPEANS.EQ.N')SCANKNT=N POINT 

CALL STORE VEL(I7SUBJ85,SCANKNT,ACTP0INT,VEL) '.STORE VELOCITIES AT END OF 
.SPC FILE 

END 


SUBROUTINE STORE_VEL((7S UB,I85,N POINT, ACTPOINT,Y EL) 

C !! STORE VELOCITIES AND FIND MAX.MIN & STORE IN CAL FILE 

INTEGERS NPOINT,T85,ACTPOINT 
INTEGERS I7SUB 
REAL*4 VEL{5000) 

VUP=0.0 

VLO=lO.**8. 

VACCUM-0. 

I7SUB=T7SUB+1 
DO 94326 |JI“l,NPOINT 

VM AX=AMAXl (VUP,VEL(I JI» 

VUP-VMAX 

VMIN=AMINl(VLO,VEL(IJl» 

VLO=VMIN 

VACCUM=VACCUM+VEL(U I) 

WR1TE(7,REC=I7SUB,FMT-53)VEL(1JI) 

94326 I7SUB=I7SUB+I 
53 FORMAT(A4) 

V A VE= VACCUM/NPOIN T 

WRTTE(8,REC=185,FMT=53)VUP 

185-185+1 

WRJTE(8,REC=I85,FMT=53)VLO 

185=185+1 

W RITE(8 ,REC=I 85,FMT=5 3 ) V A VE 


RETURN 

END 


SUBROUTINE POTNT_FROM_FREQ(FREQJIMEPERDIV,ACTPOINT) 
DETERMINE SPECTRA POINT FROM FROM SPECIFIC FREQUENCY IN SPECTRA 
SEE SUBROUTINE CENTERFREQ FOR SIMILAR PROCESSING EXPLANATION 
INTEGERS DUMMYFREQ 
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INTEGERS FREQ,ACTPOINT 

REALM TIMEPERDIV,DFXFREQ,TEMPPOINT 


DUMMYFREQ-FREQ* IE+06 
DELFREQ=(L/(2.*(10.*TIMEPERDIV))) 
T£MPPOINT=DUMMYFREQ/DELFREQ 
ACTPOTN T- JN 1 NT (TEMP POINT) 

C 

RETURN 

END 

SUBROUTINE OBTAIN MAGNITUDE SP ECTR A(TW AVE,DUMMY_ASP EC) 
COMPLEX C 1(1024),CSPEC( 1 024) 

INTEGERS ISTAT 

real* 4 DUMMY_ASPEC(1 024),PHASE( 1 024),TWAVE(5 12) 


DO 500 1=1,512 

500 C 1 (I)=CMPLX(TWAV E(1 )) 

1580 CALL LSPSFFT _COMPLEX( Cl, CSPEC, 1024,<USTAT ) 

CALL LSP$PHASE_ANGLE( CSPEC, PHASE, DUMMY_ASPEC, 1 024 ) 

RETURN 

END 

SUBROUTINE ST0RE_SPECTRA(DUMMY„ASPFX,ACTP01NT,I7,I7SUB) 

INTEGER*2 ACTPOINT 

INTEGERS R£C,I7,I7SUB 

REALM DUMM Y_ASPEC( 1 024) 


DO 94326 1 J 1= 1 .ACTPOINT 
l7SUB=((I7-l)*ACTPOINT)+IJI 
94326 WRITE(7,REC-I7SUB,FMT=53)DUMMY_ASPEC(IJl) 

17“ 17"r 1 

53 FORMAT(A4) 

RETURN 

END 

SUBROUTINE FILTER(VEL_NOTHICK > VELPREV,CFILTER,SHAPEANS,L_VEL,U_VEL, 
1 AV_VEL,NIPG,NIPJ,SCANKNTJ 1 8) 

CHARACTER CFILTER* 1 ,CRUNCH_CODE*ZSHAP£ANS* I 
INTEGERS NIPG,NIPJ,I18,SCANKNT 
REALM L_VEL,U_VEL,AV_VEL 

C ******** TEST FOR VELOCITY OUTSIDE FILTER LIMITS ***** C 
IF (SCANKNT.EQ. 1)VELPREV=AV_VEL 

IF (CFILTER.EQ. , Y , .AND.(VEL_nothickXT.L_VEL,OR.VEL_nothick.GT.UVEL))THEN 

VEL NOTHICK=VELPREV 

CRUNCHCODE-BC 
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ICRUNCH CODE_COUNT= 1 
TYPE *,’RAD POINT AT \NIPG,NIPJ 
ENDTF 

IF (ICRUNCH m CODE_COUNT.EOO)CRUNCH_CODE~ , G' 
WRITE(76JREC=I 18,FMT=53)CRUNCH_CODB 
118=118+1 

WRITE(7 6,REC=T 1 8,FMT=53 )N 1 PG 
118=118+1 

WRITE(76,REC=ll8,FMT=53)NtPJ 

118=118+1 


53 FORMAT(A4) 

RETURN 

END 

SUBROUTINE NOTHICK_VEL_CALC(DELTATJWOTAU,VEL_NOTHICK) 

VELWATER=0.149 ICM/USEC 

VEL NOTHICK=VELWATER*((DJiLTAT/TWOTAU)+l) 

RETURN 

END 


SUBROUTINE CORR(Bl,B2,CXX) 

C CORRRELATE TWO WAVEFORMS 

COMPLEX’^ NCI 1(5 !2),NC12(5 12), NCI 3(5 1 2) 
COMPLEX* 8 NC 14(5 12),NC 1 5(5 12),NC22(5 1 2) 

REAL*4 NC16(5 12),Bl(5 12),B2{5 12), CXX 
INTEGERS STATUS 

DO 333 1=1,512 

NC 1 1(I)=CMPLX(R 1 (I» 

NC22(I)=CMPLX(B2(I» 

333 CONTINUE 

CALL LSP$FFT_COMPLEX(NC 1 1 ,NC 12,5 12,0, STATUS) 
CALL LSP$FFT_C0MPLEX(NC22,NC1 3,5 1 2,0, STATUS) 

DO 77710=1,512 

NCI4(IO)=NC12(IO)*CONJG(NC13(IO)) 

777 CONTINUE 

CALL LSPSFFT COMPLEX(NC14,NC15,5 12,0, STATUS) 

IF (.NOT. STATUS) CALL LIB$STGNAL(%VAL(STATUS» 

AMAX=0.0 

CXX=0.0 

DO 765 1=1,256 

NCI 6(i+256)=RE AL(NC 1 5(1)) 

NC 1 6(1)=REAL(NC 1 5(1+256)) 

765 CONTINUE 
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DO 881 1=1,512 
I F(NC 1 6(1) .GT. AMAX)THEN 
AMAX=NC 1 6(1) 

CXX=T 

ENDIF 

881 CONTINUE 

CXX=CXX-257. 

RETURN 
END 

SUBROUTINE MCORR(Bl,B2,CXX) 

CORRRELATE TWO WAVEFORMS !!1 modified by using absolute value 
of minimum of correlation function 

to take into account phase inversion 
of B2 w/ respect to B 1 such as what 
happens w / PMCs 

COMPLEX* 8 NC 11(51 2),NC 1 2(5 1 2),NC 13(5 12) 

COMPLEXES NC14(5 12),NC 1 5(5 12),NC22(512) 

REALM NC 1 6(5 12), B 1 (5 1 2),B2(5 12),CXX 
!NTEGER*4 STATUS 

DO 333 1=1,512 
NCI 1(I)=CMPLX(B 1 ( 1 )) 

NC22(I)=CMPLX(B2(I)> 

333 CONTINUE 

CALL LSP$FFT_COMPLEX(NCl 1, NCI 2, 5 12,0, STATUS) 

CALL LSP$FFT_COMPLEX(NC22,N C 1 3,5 12,0,STATUS) 

DO 777 10=1,512 

NC 14(IO)=NC 12(10)*C0NJG(NC 1 3(10)) 

777 CONTINUE 

CALL LSP$FFT_C0MPLEX(NC14,NC 15,5 1 2,G,STATUS) 

IF (.NOT. STATUS) CALL LIB$SIGNAL(%VAL(STATUS)) 

AMAX=0.0 

AMlN=L0E+6 

CXX=0.0 

DXX=0,0 

DO 765 1=1,256 

NC 1 6(I+256)=REAL(NC 1 5(I» 

NCI 6(l)=REAL(NC15(I+256» 

765 CONTINUE 

DO 881 1=1,512 


C IF(NC16(!).GT.AMAX)THEN 
C AMAX=NC16(I) 

C CXX=1 

C ENDIF 


C CORRRELATE TWO WAVEFORMS ! ! ! modified by using absolute value 
C of minimum of correlation function 


C 

C 

C 

c 

c 


31 
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to take into account phase inversion 
of B2 w/ respect to B 1 such as what 
happens w/ PMCs 


IF(NC1 6(I).LT.AMIN)THEN 
AMIN=NC16(I) 

DXX=1 
EM DIF 

881 CONTINUE 


CXX=DXX 


CXX=CXX*257. 

RETURN 

END 
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CC PROGRAM: NOTHICK JMAGBMAKER.FOR 

Read Data file written by CRUNCH and place data into 
"Grmnell-ready" files 


Header: X3,Y3 Dimension of scan ( SCAN DIRECTION, INDEX DIRECTION ) 
Nl, N2 Number of points ( SCAN, INDEX ) 


INTEGER+2 ACTPOINT 
BYTE BFILEEXT(23) 

CHARACTER FILENAME M4,FILEEXT*23, DRIVE* 1 
CHARACTER CALM AME*44,SELFREQANS* 1 .SHAPE AN S* 1 ,DTSKT*5 
CHARACTER HEAD*23 ,DIR* 1 6 
EQUIVALENCE ( FILEEXT, BFILEEXT ) 

JMARK=0 
185=1 

OPEN (UNIT=10, FILE~NOTHICK_DADQ1.DAT', STATUS-OLD’, 
i FORM-FORMATTED') 

READ (10,34572) FILEEXT 
READ (10,34572) DRIVE 
34572 FORMAT (A) 

4321 CONTINUE 

HEAD=FILEEXT 

IF {DR1VE.EQ. , A , )D1R-DUC2:[R0TH.IMAG]' 

IF (DRIVE.EQ. , A’)D1SKT-’DUC2: 1 
IF (DRIVE.EQ.'C^DIR-DUCO^ROTH.IMAG] 1 
IF (DRI VE.EQ.'C')DISKT~ DUCO: ' 

ICHA=0 

DO 1=1 ,LEN(F1LEEXT) 

IF (1CHA.NE.0) GOTO 98564 
IF (BF1LEEXT(I).EQ.32) THEN 
ICHA-i 

FILEEXT (ICHA: 1CH A+4)=' . SPC‘ 

ENDIF 

98564 END DO 

FILENAME = ’[ROTH.DAT AJ7/FILEEXT 
K=0 

20 K=K+1 

IF (FILEEXT(K:K),NE.V) GOTO 20 
FILEEXT(K+1 :K+3) - ’CAL’ 

CALNAME - ’[ROTH .DAT A]7/FILEEXT 

CALL DONDISP1S (ICHA, FILENAME, CALNAME, HE AD,DIR,DISKT, 
SHAPEANS.SELFREQAN^ACTPOINT) 


v/ 
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CLOSE(IO) 

END 


SUBROUTINE DONDISP1S (ICHA,FILENAME.CALN AME.HEAD, DIR,DI SKT, 
1 SHAPEANS,SELFREQANS,ACTPOINT) 


Read Data file written by CRUNCH and write values to Grinnell 
-UNCORRECTED & CORRECTED DATA IS DISPLAYED— 

Header: X3 , Y3 Dimension of scan ( SCAN , IN DEX ) 

Nl, N2 Number of points ( SCAN, INDEX ) 


INTEGERS REC,SCANPOS,XXX 
INTEGERS IXLENGTH_A,DCLENGTH,ACTPOINT 
INTEGERS X3,Y3,Y87,N1,N2,NNNNN2,NNNNNI 
INTEGERS DENS,NFREQ,AVES,X4,Y4J87 
INTEGERS COUNT, NCOU N T,N SC AN,N IN DX 
INTEGERS C(1024),D(460) 

INTEGER*2 SCANDISTJNDXDIST,TRFREQ,ENDFREQJ 1 6 
INTEGERS D_SHAPE(460),189,I6,I 1 5 

INTEGERS ISLANT1JSLANT2,ISLANT3,ISLANT4,JXXX,IABC,JABC,SCANKNT 

INTEGERS PSEUDOGRNL(460, 1 024),SELFREQ 

REAL FREQ, MEAN 

REAL RNUP(7),RNLO(7),CBUFFER(S 1 ) 

BYTE BSCALEMARK(4),CHFREQ(3) 

BYTE BSUFFREQ 

BYTE BFILEN AME(44).BN AMELAB(20),BD AT A_ST ATUS(40) 

CHARACTER FILENAME*44,FILEEXT*23,SELFREQANS* 1 ,FFLAG* 1 
CHARACTER CALNAME*44, CHEADER*32.SHAPEANS* 1 ,ZIG ANS * 1 

CHARACTER NAMELAB*20,DATA_STATUS*40,NEWHEADER*35,ACTUAL_SELFREQANS*1 
CHARACTER CORNAME*I5,CHFREQ_C*3,PHASE*I,PHASEl*l 
CHARACTER SCALEMARK*4,SUFFREQ*4,HEAD + 23,DIR*16,DISKT*5 

EQUIVALENCE ( CIIFREQ, CHFREQ C ) 

EQUIVALENCE ( BNAMELAB, NAMELAB ) 

EQUIVALENCE ( BDATA_STATUS, DATA_STATUS ) 

EQUIVALENCE ( BSCALEMARK, SCALEMARK ) 

EQUIVALENCE ( BSUFFREQ. SUFFREQ ) 

SCALEMARK-' 1 MM' 

SUFFREQ- MHZ' 

TAKE CHARACTER ’’FILENAME” AND DECODE TO BYTE ”B FILENAME" 


DECODE (44, 5 33, FILENAME) BFILENAME 
533 FORMAT(44Al) 


V2L 


C 
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C DETERMINE SHORTENED (-.SPC) FILENAME POSITION 
C 

16=1 

115=1 

116=1 

KMARK-0 

DO 4321 1PR=12,44 

IF (KMARK.EQ.l) GOTO 4321 

IF (BFILENAME(IPR).EQ.46) THEN 1PERIOD/DOT IS ASCTI 46 

1END=IPR 

KMARK=1 

ENDIF 

4321 CONTINUE 


34572 FORMAT (A) 

34573 FORMAT (I) 

51 FORMAT(A32) 

52 FORMAT(A2) 

53 FORMAT(A4) 


OPEN( unit=15, filc=FILENAME(1 :IEND-I )/Adatl2', status='OLD\ 

+ ACCESS=’DIRECT\RECORDTYPE^FIXED\RECL=2, 

+ form^FORMATTED', ORGANIZATION- SEQUENTIAL' ) 

OPEN( unit=l6, file=FILENAME(l ilEND-iy/'.ciatRT, $tatus='OLD\ 

+ ACCESS= , DIRECT,RECORDTYPE=TIXED , ,RECL=4 > 

+ form- FORMATTED', ORGANIZATION -SEQUENTIAL’ ) 

C OPEN( unit=26, fi le=F!LENAME( 1 :IEND- 1 )7/\datRR4', status='OLD', 
C + ACCESS='DIRECT , ,REC0RDTYPE= , F1XED I ) RECL=4, 

C + form- FORMATTED 1 , ORGANIZATION- SEQUENTIAL' ) 

OPEN( unit=66, file=FILENAME( 1 :1EN D- 1 )//'.datCI I 1 , status='OLD’, 
+ ACCESS- DIRECT, RECORDTYPE= :, FIXED\RECL=32, 

+ form- FORMATTED’, ORGANIZATION -SEQUENTIAL’ ) 


C**************m\T=7„S?C FILE FOR SPECTRA, vel, FOR DIRECT ACCESS********C 
C 

C OPEN( unit=7, file -FILENAME, status^'OLD', ACCESS-' DIRECT, 

C - RECORDTYPE- FIXED’, R£CL=4, FORM-FORMATTED', 

C 4 ORGANIZATION-SEQUENTIAL’) 

OPEN (UNIT- 1 4,FILE=’N0TH1CK_SELFREQ 1 .DAT, 

1 STATUS -’OLD'.FORM- FORMATTED’) 


79193 READ(66,REC=I6,FMT-5 1 )CHEADER 
16=16+1 

READ(66.REC=I6 ,FMT=5 1)SHAPEANS 
16=16+1 


V 3 



103 


5 , 629,865 


104 


READ(66,REOT6, FMT=5 1 )ZIGANS 
16=16+1 


READ (14.98652) F1LEEXT 
READ (14,98652) DRIVE 
READ (14,98652) PHASE 
READ (14,98652) PHASE! 

READ (14,98652) DIFFANS 
READ (14,10030) BUFLENGTH 
READ (14,10030) BUFVEL 
READ (14,10030) TRDIAM 
CLOSE (14) 

98652 FORMAT (A) 

10020 FORMAT (I) 

10030 FORMAT (F) 

89714 READ( 1 5,REC=T 1 5,FMT=52)SC ANDIST 

I15=I15H 

READ(15,REC=I1 5,FMT=52)1NDXDIST 
115=115+1 

RE AD( 1 5 ,REC=1 1 5 .FMT=52)NSCAN 
115*115+1 

READ( 1 5,REC=I 15,FMT=52)NINDX 
115=115+1 

READ( 1 5,REC=I 1 5,FMT=52)islant1 
115=115+1 

READ(1 5 f REC=I 1 5,FMT=52)isIant2 
115=115+1 

READ(15,REC=1 1 5,FMT=52)isIant3 
I15-I15+1 

RE AD(1 5.REC=I 15,FMT=52)islant4 
115=115+1 

READ(1 5,REC=115,FMT=52)DENS 
115=115+1 

READ(15,REC=I 1 5,FMT=52)TRFREQ 
115=115+1 

READ( 1 5,REC=1 1 5,FMT=52)AVES 
115=115+1 

READ(15,REC=I 15.FMT=52)JXXX 
115=115+1 

C ARRAYSIZE=(NSCAN)*(JXXX) !IF SCAN FINISHED FROPERLY,JXXX=NlNDX 


ENCODE (3, 19878, CHFREQ) TRFREQ .‘ANALYSIS FREQUENCY 
19878 FORMAT (13) 

READ(16,REC=I 16,FMT=53)TIMESET 
116=116+1 

ENDFREQ=2. 5* TRFREQ 

CALL POINT FROM_FREQ(ENDFREQ,TIMESET,ACTPOINT) 
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DO 9999 ICHOICE-1,1 
187=1 


ACCUM=0. 

NBAD=0 

IF (ICHOICE.EQ. 1) THEN 

N AMELAB - VELOCITY CM/US' 

CORN AME- _ V EL_C J//CHFREQ_C 
N EW HE ADER- CORRECTED DATA FOR VELOCITY’ 
NEWU=25 
NEWU 1=45 
NEWU2=<59 

ENDIF 


OPEN (UNIT=NE WU I ,FILE=DISKT//FILENAME( 1 :IEND- 1 )//CORN AME# 
+ , .CORHEADER‘jSTATUS=‘NEW' J FORM- UNFORMATTED') 

12345 IF (JMARK.EQ. I) GOTO 90 1FILES ARE ALREADY OPEN 

30 CONTINUE 

90 CONTINUE 

IF (JMARK.EQ.1) GOTO 89306 

10149 CONTINUE 


IF (SHAPE ANS.NE.'Y’)GOTO 893 1 7 
OPEN(liNIT=52,FILE=FILENAME(l:ICHA+10)//^DATHS , , 

1 STATUS=OLD',FORM -UNFORMATTED 1 ) !INFO ON WHETHER WE ARE ON HOLDER/SAMPLE 


89317 SCANDIST=ABS(SCANDIST) 

AD=NSCAN 

XREAL = RE AL(ABS(SCAN DI ST)) 

YREAL = REAL(INDXDIST) 

NPTS = NSCANUXXX !!t!lF SCAN COMPLETED PROPERLY, JXXX=NINDX 

IF (SHAPEANS.EQ.’Y')THEN 
DO 5697 1=1, NPTS 

5697 READ(UN1T=52, END=893 1 8)IABC,JABC,FFLAG,SCANKNT 
ENDIF 

893 18 PTS = REAL(NSCAN)*REAL(JXXX) 

X3=ABS(SCANDIST) 

Y3=INDXDIST 

NI-NSCAN 

N2=JXXX 

IF (I.EQ.l)THEN 

NNNNN1-N1 

NNNNN2=N2 

ENDIF 
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XREAL = REAL(ABS(X3)) 
YREAL = REAL(Y3) 

FREQ = REA L(N FREQ) 

PTS = REAL(N1)*REAL(N2) 


Proportion GRINNELL window similar to scan dimensions 


200 Y87 - ININT{ REAL(Y3)*8./7. ) 

IF( YS7.GE.X3 )THEN 

R1 = FLOAT(X3)/FLOAT(Y3) 
X4 = 400 

Y4 - 1NINT(400./Rl) 

IF (Y4.GE.480)THEN 
X4=300 

Y4 = ININT(300./R1) 

ENDIF 

IF (Y4.GE.480)THEN 
X4=200 

Y4 - ININT(200./R1) 

ENDIF 

IF (Y4.GE.480)THEN 
X4=I00 

Y4 = ININT(I00./R1) 

ENDIF 

IF (Y4,GE.4S0)THEN 
X4=50 

Y4 = ININT(50./R1) 

ENDIF 

IF (Y4.GE.480)THEN 
X4-25 

Y4 = ININT(25./R1) 

ENDIF 

IF (Y4.GE.480)THEN 
X4-10 

Y4 = ININT(10./R1) 

ENDIF 

ELS£IF( Y87.LT.X3 ) THEN 
Rl - FLOAT(Y3) f FLOAT(X3) 
Y4 = IFIX( 400. *RI ) 

X4 = 400 

IF (Y4.GE.480)THEN 
X4=300 

Y4 — 1N1NT(300.*R1) 

ENDIF 

IF (Y4.GE.480)THEN 
X4=200 

Y4 = ININT(200.*R1) 

ENDIF 

IF (Y4.GE.480)THEN 
X4=1 00 

Y4 = ININT(100.*RI) 

ENDIF 

IF (Y4,GE.480)THEN 
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X4=50 


Y4 = ININT(50.*R1) 
ENDIF 

IF (Y4.GE.480)THEN 
X4=25 


Y4-ININT(25.*R1) 

ENDIF 

IF (Y4.GE.480)THEN 
X4=1C 


Y4^ 1NINT(10.*R1) 
ENDIF 


ENDIF 


CC 

XD2G - REAL(N 1 )/RE AL(X4) 
YD2G - REAL(N2)/REAL(Y4) 

XD2X -REAL(X4)/REAUN1) 
YD2Y = REAL(Y4)/REAL(N2) 

I_INT=JINT (XD2X)- 1 
JJNT=JINT(YD2Y)+1 


89306 CONTINUE 

10300 format{ //*$ Enter choice :') 

54322 FORMAT (A) 

43 1 9 FORMAT C+y VALUH^F) 

TYPE * ‘ ' 

TYPE 'Working...' 

TYPE V ' 

43191 FORMAT (M '/BAD PTS/TOT. PTS:\I5,7\I4, 
+ ' CURRENT VALUE= \F1 1.2) 

7493 FORMAT (FI 1.2) 

7494 FORMAT(I) 


C 

C RETURN MAX/MTN OF DATA FOR GRINNELL DISPLAY & PRODUCE UNF. 

C OF GRINNELL DATA (CBUFFER) 

CALL RETMAXMINFILTER (FILENAME, CORNAME,CALNAME,NAMELAB,DISKJ\IEND.JXXX, 

1 

NSCAN,SCANKNT,SHAPEANS,ZIGANS : ICH0ICE,TIMESET,I89,C0UNT,ACTP01NT,NEWU,GRB0T, 

I GRTOP,MEAN) 

7487 COUNT-NCOUNT 
RMIN-GRBOT 
RMAX=GRTOP 
187—1 


^7 
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RNLO(ICHOICE)=RMIN 
RNUP(ICHO 1C E)=RM AX 

DO 500 KI-1 ,JXXX 
DO 49721 J 1=1, NSC AN 
NIPI=JI 
NIPJHCI 

SCANPOS=((NIPJ-l)*NSCAN)+NIPI 

I87=SCANPOS 

READ (NEWU,REC=IS7,FMT=53) CBUFFER(Jt) 

49721 187=187+1 

462 FORMAT (FI 1.2) 

CALL SCGREAL (CBUFFER, C, Nl,GRBOT,GRTOP,MEAN) IGREY SCALED 
DAT A STATUS- CORRECTED DATA FOR 7/NAMELAB 
C 

C INTERPOLATE DATA GRID (EX.8 1X8 1 ) FOR X4 x Y4 GRINNELL DISPLAY 

C 

CALL 1NTR ( C,D,NSCAN,X4 ) 

C IF (KI.EQ.2Q)TYPE *, , D=',D 


C ********** mod TO GET RID OF FALSE BORDER ***************** C 
IF (shapeans.eq:y')then 

IIXCOUNT=0 

DO 555 10P=1,X4 ’FORWARD X-DIRECTION 
IF (D (IOP) . EQ. 1 )GOT O 555 
IIXCOUNT=lIXCOUNT+ 1 
IF (1IXCOUNT.EQ. 1)THEN 
DO 655 IJQ=1,I_INT 
655 D(IOP+IJQ-1)=l 

ENDIF 

555 CONTINUE 

IIXCOUNT=0 

DO 755 IOP=X4,L-l IBACK.WARD X-DIRECTION 
IF (D(!OP).EQ.l)GOTO 755 
IIXCOUNT=IIXCOUNT+ 1 
IF (IIXCOUNT.EQ. 1)THEN 
DO 855 IJQ=1,1JNT 
855 D(IOP-IJQ+I)=l 

ENDIF 

755 CONTINUE 

C IF (K1.EQ.20)TYPE */D=\D 

ENDIF 


C ********** MOD TO GET RID OF FALSE BORDER ********** C 
DO 480 1,X4 

480 PSEUDOGRNL( J,N1PJ ) - D( J ) 

500 CONTINUE 
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cc 


OPEN (UNIT=NEWU2, FILE=DIR//HEAD//CORNAME//\PDK\STATUS- f NEW\ 
+ FORM- UNFORMATTED') 

WRITE(NEWU2) X4,Y4,N1,N2,ICH0ICE 
WRITE(NEWU2) SCANDIST.TRFREQ 
WRITE(NEWU2) RMIN,RMAX,MEAN 

TYPE * t X4,Y4,N 1 ,N2,ICHOICE 
TYPE *,SCANDIST ? TRFREQ 
TYPE *,RMIN,RMAX,MEAN 


IXLENGTH-0 
605 DO 800 IJ= 1 , X4 

IX = IJ-1 

DO 540 J=1,N2 

540 C( J ) = PSEUDOGRNL( IJ,J ) 
CALL 1NTR ( C AJXXX,Y4 ) 
C TF (IJ.EQ.200)TYPE =VD=\D 


C ********** MOD TO GET RID OF FALSE BORDER ***************** c 
IF (SHAPEANS.EQ.’Y'yrHEN 

IIXCOUNT=0 

DO 1555 IOP=l,Y4 ’FORWARD Y-DIRECTION . 

IF (D (IOP ) . EQ. 1 )GOTO 1555 
ITXCOUNT=HXCOUNT+l 
IF (1IXCOUNT.EQ. 1)THEN 
DO 1655 IJQ=i,J_INT 
1655 D(10P+IJQ-1)=1 

ENDIF 

1555 CONTINUE 
HXCOUNT=0 

DO 1 755 TOP=Y4 J l.-1 ! BACK WARD Y-DIRECTION 
IF (D(IOP).EQ.l)GOTO 1755 
IIXCOUNT=TIXCOUNT+ 1 
IF (IIXCOUNT.EQ. 1 )THEN 
DO 1855 1JQ=I,J_1NT 
1855 D(IOP-IJQ+ 1 )= 1 

ENDIF 

1755 CONTINUE 
C IF (IJ.EQ.200)TYPE *,D-,D 

ENDIF 


C ********** MOD TO GET RID OF FALSE BORDER ********** c 


Hi 



115 


5 , 629,865 


116 


c#J(t ************ M0DIFTCA . ri0N for COMPLEX SHAPE ***********c 


IF (SHAPEANS.NE.'Y')GOTO 675 'NOT A COMPLEX SI IAPE 
IF (SHAPEANS.EQ.'Y’)THEN 


XXX=1 
IXSTART=IJ 
IYLENGTH=0 
NKEEP-0 
DO 8970 J-1,Y4 

IF (D(J).LT.2)GOTO 8970 ! 1 IS CODED FOR LUCITE (IN SCGREAL) 

IF (D(J).GE.2)THEN ! SAMPLE DATA 
IF (NKEEP .EQ.0)l Y ST ART=J 
NKEEP= I 

C**** FURTHER MOD FOR GEORGE'S DISCONTINUITY ON YBCO HEX SAMPLE ****** 
IF (IYLENGTH.GT. L AND.D(J).GE.2.AND.D(J-1 ).LT.2)THEN 
8875 XXX=XXX+] 

I YLENGTH-IYLENGTH+ 1 
D_S HAPE(IYLENGTH)= 1 
IF (D(J-XXX).LT.2)GOTO 8875 
ENDIF 

C**** FURTHER MOD FOR GEORGE'S DISCONTINUITY ON YBCO HEX SAMPLE ****** 
IYLENGTH=IYLENGTH+1 
D_SHAPE(IYLENGTH)-D(J) 

ENDIF 

8970 CONTINUE 

C TYPE VIXSTART=\IXSTART/IYSTART= , ,IYSTART, , IY[,ENGTH=\ 

C 1 IYLENGTII 

IF (IYLENGTH.EQ.O)GOTO 800 I ALL LUCITE IN THIS COLUMN 
IXLENGTH=IXLENGTU+ 1 
WRITE(NEWU2)IXSTART,IYSTART,IYLENGTH 
DO 8979 IJK=UYLENGTH 
8979 WRITE(NEWU2)D_SHAPE(IJK) 

GOTO 800 
ENDIF 

£ ********* ********* 11^*4 ******** ********¥¥****9******* ****** Q 

675 WR1TE(NEWU2) D 

800 CONTINUE 


TYPE V • 

TYPE HIGHEST ',NAMELAB,’ IS ',RMAX 
TYPE V LOWEST \NAMELAB,’ IS \RM!N 
TYPE V ’ 

WRITE (NEWUl)NEWHEADER 
WRITE (NEWU1)RMIN,RMAX 
WRITE (NEWUI)MEAN 
WRITE (NEWU1) 1SCALEF 


SO 
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97364 FORMAT (A) 

CLOSE (NEWU) 

CLOSE (NEWU1) 

CLOSE (NEWU2) 

TYPE *.DONE CLOSING \ HEAD//CORNAME//'.PDK' 
JMARK=i 

IF (SHAPEANS.NE.TOGOTO 9999 
5001 REWIND(52) 

9999 CONTINUE 

IXLENGTH_A=IXLENGTH 

WRITE( 1 5,REC=I1 5,FMT=52)IXLENGTH A 

CLOSE(I5) 

CLOSE(52) 

87124 CONTINUE 


CLOSE (6) 

RETURN 

end 


C 

C GRNLBASE.FOR 

C - A set of subroutines useful for applications 

C with the GRINNELL image processor. 

C GRINIT must be called in main program before calling these 
C subroutines. 

C 

C David B.Stang latest update 19-APR-88 


SUBROUTINE SCGREAL( CC,DD, N, BOT,TORMEAN ) 

Convert values in CC to 1 thru 255 ( same as SCGR, with CC REAL ) 

CC( N ) Input Data 
DD( N ) Output Data 
EOT Desired Min for Gray Scale ( units of CC ) 

TOP Desired Max for Gray Scale " 

INTEGERS N 
1NTEGER+2 DD(N) 

REAL CC(N),MEAN 

RAN - TOP - BOT 
IF (RAN.EQ.0.)RAN=1 . 

DO 2000 J- I,N 
CO - CC{J) 

Q “ ( CO - BOT ) / RAN 
IF ( Q.LE.O.) Q = 0. 

IF ( Q.GT.l.) Q - I. 
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Q - Q*25I. 

C DD(J)=NINT(Q)+1 
C***** FOR MY MODIFICATION ***************c 
IF (CO.EQ.O.)THHN 
DD(J)= I 
GOTO 2000 

ELSEIF (CO.NE.O.)THEN 

DD(J)-NINT(Q)+2 

END1F 

**** ******* *********** ************* 

2000 CONTINUE 

RETURN 

END 


CC SUBROUTINE INTR 
CC 

CC Interpolate values in C(N2) to "fit" into D(NN) 
C 

C SUBROUTINE INTR ( CC,DD,N2,NN ) 

C 

C INTEGERS N2,NN 

C INTEGERS CC(N2), DD(NN) 

C 

C RN * FLOAT(N2) / FLOAT(NN) 

C 

C DO 3000 K = 1,NN 
C R = FLOAT( K )* RN 

C IR = IFIX( R ) 

C IR “ IINT( R ) 

C IRPI = IR+I 

C RR - R - FLOAT( IR ) 

C IF( 1R.LT.1 ) IR=1 

C IF( NN.GT.N2 .AND. IRP I .GT.NN ) THEN 

C IR=NN-1 

C 1RP1-NN 

C ENDIF 

C CO “ FLOAT( CC(IR) ) 

C Cl = FLOAT! CC(IRPI) ) 

C DO = ( l.-RR )*C0 + RR*C1 

C 3000 DD(K) = IFIX(DO) 

C 3000 DD(K) = IINT(DO) 

C 

C RETURN 

C END 


SUBROUTINE RETMAXMINJFILTER(FILENAME,CORNAME,CALNAME,NAMELAB 5 DISKT, 
1 lEND.JXXX,NSCAN,SCANKNT,SHAPEANS,ZIGANS,ICHOICE.TIMESET,I89, 

1 COUNT, ACTPOINT,N EWU,GRBOT,GRTOP,MEAN) 


sx 


INTEGERS REC,SCANPOS,I7 
INTEGER*2 NSCAN,ACTPOINT 
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INTEGER*! 185.187, SCANKNT 
INTEGER*! IABC,JABC,JXXX,I89 
INTEGER *2 NCOUNT 
INTEGER*2 ZIGCOUNT 
INTEGER*2 NBAD,COUNTXCOUNTER 
REAL DAT ATHIN G( 1 024,9) 

REAL MEAN,VEL(5000) 

REAL CBUFFER(81) 

CHARACTER FILENAME*44,CORNAME*I 5.NAMELAB*20,FFLAG* l 
CHARACTER SHAPEANS* 1,ZIGANS*1,DISKT*5,CALNAME*44 
INTEGER *4 SCANPOS1 


TYPE *,NAMELAB 
LCOUNTER-I 
ZIGCOUNT=l 
ACCUM-0. 

SCANPOS=0 

SCANPOSl-O 

NBAEHCOUNT 

187-1 

\ \ XMAX=0. 

XMIN= 1 0000000. 

PTS = REAL(NSCAN)*REAL(JXXX) 

NPTS==NSCAN*JXXX 

CAL FILE WITH M AX/MIN/M E AN 

OPEN( unit-8, file-CALNAME, status-'OLD’, ACCESS-’DIRECT’, 
+ RECORDTYPE- FIXED', RECL=4,FORM='FORM ATT FX>', 

+ ORGANIZATION-'SEQUENTIAL' ) 

185=1 


read(8,REC=I85,FMT=53)VUP 

185=185+1 

read(8,REC=I85,FMT=53)VLO 

185=185+1 

read{8, REC-1 8 5 , FMT=53) V A V E 


IF (SHAPEANS. EQ.'Y')THEN 

NPTS-SCANKNT 

REWIND(52) 

ENDIF 

I7=NPTS*4*ACTPOINT+l [START OF VELOCITY READS IN .SPC FILE 

SPC FILE (UNIT=7) CONTAINS ALL VELOCITY & SPECTRA DATA 

OPEN (UNIT=7,FILE=DISK.T//FILENAME.STATUS- OLD', ACCESS- DIRECT, 
+ RECORDTYPE- FIXED', RECL=4, FORM- FORMATTED’, 

+ ORGANIZATION-'SEQUENTIAL') 


OPEN (UNIT-NE WU,F1LE=D1 S KT//FILEN AME( 1 : IEND- 1 )//CORNAME//'.COR\ 
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+ STATUS-NEW'jACCESS^'DlRHCT, 

+ RECORDTYPE-FIXED\RECL=4,FORM='FORMATTED', 

+ ORGANIZATION-'SEQUENTIAL') 


DO 500 NTPJ=1 > JXXX 
410 DO 450 NIPI=I ,NSCAN 

C*** + * GET INFO FROM HOLDER. DAT (FFLAG-S' -> SAMPLE) (FFLAG-’H’ -> HOLDER) 
KCON=KCON+ 1 
J=NTP1 

SCANPOSI=((NIPJ-l)*NSCAN)+NIPl [TRUE SCAN POSITION 
IF (SHAPE AN S .NE. ’ Y’)GOTO 87123 
READ(52)IABC,JABC,FFLAG 3 SCANKNT 
C TYPE *;XH=MABC/ YH-,JABC,' FFLAG-,FFLAG 

IF (FFLAG.NE. , S , )GOTO 8702 [TRANSDUCER IS ON HOLDER/GOTO NEXT PT. 
87123 $CANPOS=SCANPOS+l 

CC ! ! ’SCANPO S' FOR FILE RETRIEVAL FOR CIRCLE DISK SCANS 
READ(7,REC-I7 3 FMT^53)VEL(SCANPOS) 

17-17+ 1 

DATATH1NG(J,5)=VEL(SCANP0S) 

N PTS=(NIP J- 1 )*N SC AN+N I P I 

IF(NPTS.EQ.NSCAN*JXXX)NPTS=NSCAN*JXXX 


C IF (ICHOICE.EQ.l)THEN 

IF(SCANPOS.EQT)DATATHINGPREVIOUS-VAVE 
C IF (DAT ATH IN G ( J, 5).LT ,( V LO+ VLO *0.001). OR . 

C I DATATHTNG(J,5>.GT.(VUP-VUP*O.OOI))DATATHING(J,5)-DATATHINGPREVIOUS 
C ENDIF 

IF (SHAPEANS.NE.'Y’)GOTO 87651 

8702 IF (FFLAG.NE. r S')THEN !SET TO 0 // TRANSDUCER IS ON HOLDER/GOTO NEXT PT. 

CBUFFER(J)-0,0 

GOTO 450 
ENDIF 

87651 CBUFF£R(J)=DATATHING(J,5) 

390 CONTINUE 

Cm-****kEEP TRACK OF MAX & MIN LIMITS OF CORRECTED DATA********C 
C **** FOR FFLAG.NE/S' ™ WE SKIP THIS STEP SO THAT 'O’ IS NOT CONSIDERED 
RMAX=AM AX 1 (CBUFF£R(J),XMAX) 

RMIN^AMIN 1(CBUFFER(J),XM1N) 

XMAX=RMAX 

XMJN=RMIN 


ry 
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DATATH1NGPREVT0US=CBUFFER(J) 


450 CONTINUE 

IF (ZIGANS.EQ.'N , )THEN !NOT A ZIGZAG SCAN - A NORMAL SCAN 
DO 49720 JI=1,NSCAN 

WRITE (NEWU,REC=IS7,FMT=53) CBUFFER(JI) 

ACCUM=ACCUM + CBUFFER(JI) 

49720 187=187+1 

ELSE1F (ZIGANS,EQ.'Y\AND.ZIGCOUNT.EQ.l)THEN !Z1GZAG/LEFT->RIGHT 
DO 49721 JI=LNSCAN 

WRITE (NEWU,REC=I87,FMT=53) CBtJFFER(JI) 

ACCUM-ACCUM + CBUFFER(JI) 

C TYPE *;CBUFFER(JI)-’,CBUFFER(JI) 

49721 187=187+1 

BLSE1F (ZIG ANS.EQ.'Y'. AND.ZIGCOUNT.EQ.O)TI IEN IZfGZ AG/RIGHT->LEFT 

DO 49722 JI=NSCAN.1,-1 

WRITE (NEWU,REC=I87 1 FTV1T=53) CBUFFER(JI) 

ACCUM=ACCUM + CBUFFER(JI) 

C TYPE *,'CBUFFER(J1)=',C BUFFER(JI) 

49722 187=187+1 
ENDIF 


C***** FOR ZIGZAG SCANS - TO ALTER ROW STORAGE ***********C 
IF (ZIGAN S.EQ.'Y*. AND. Z1 GCOUNT.EQ. 1 )THEN 
ZIGCOUNT=0 
GOTO 500 

ELSEIF (ZIGANS.EQ.'Y'.AND.ZIGCOUNT.EQ.O)THEN 

ZIGCOUNT=l 

ENDIF 

C500 TYPE VY=',JABC.'CBUFFER-,CBUFFER 
500 CONTINUE 
5001 GRBOT = RM1N 
GRTOP = RMAX 

TYPE V ' 

TYPE VRMIN-. GRBOT 
TYPE VRMAX=’,GRTOP 
MEAN = ACCUM/SCANPOS 
TYPE VMEAN=\MEAN 
IF ( SHAPE AN S .NE.’Y ')G OTO 88888 
REWIND(52) 

88888 CONTINUE 
53 FORMAT (A4) 

43191 FORMAT ('+7BAD PTS =\I5/ AT LOCATION (SCANPOS) =\I5) 

RETURN 

END 

SUBROUTINE P01NT_FROM FREQ(SELFREQ,TlMEPERDTV,ACTPOINT) 
C 


5:r 
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DETERMINE SPECTRA POINT FROM FROM SPECIFIC FREQUENCY IN SPECTRA 

SEE SUBROUTINE CENTERFREQ FOR SIMILAR PROCESSING EXPLANATION 

INTEGERS DUMMYFREQ 

INTEGERS SELFREQ,ACTPOINT 

REAL*4 TIMEPERDI V, DELFREQ,TEMPPOINT 


DUMMYFREQ-SELFREQ* 1 E+06 
D ELFR EQ-( L/(2. *( 1 0 *TIMEPERDIV)» 

TEMPPOINT^DUMMY FREQ/D ELFREQ 
ACTPOINT=JNINT(TEMPPOINT) 

RETURN 

END 

SUBROUTINE CENTERFREQCDUMMYP.TIMEPERDTV,ACTFREQ) 

DETERMINE FREQUENCY FROM POINT IN B 1(F) SPECTRA 

INTEGERS DUMMYP,ACTFREQ 

REALM TIMEPERDI V , DELFREQ,TEMPFREQ 

WE KNOW THAT DELFREQ=( 1/(2N*DELTI ME)) 

=( 1 /(2N *( 1 0*TIMEPERDI V/N))) 

-( I/(2*(10*TIMPERDIV)» 

WHERE N - 5 12 POINTS/ACQUISITION (SCREEN) 

TIMEPERDIV = TIME SETTING PER DIVISION (EX. 50 NSEC) 

2 = WE EXTENDED 512 ARRAY TO 1 024 ARRAY FOR PROCESSING 
10 = THE NUMBER OF DIVISIONS/SCREEN 

FOR EXAMPLE: 

WE CAN DETERMINE THE FREQUENCY FOR WHICH THE POINT "MAXB1P" 
CORRESPONDS TO BY JUST MULTIPLYING (DELFREQ) * (MAXBIP) 

DELFREQ=( 1 ./( 2 *( 1 0.*TIMEPERDIV))) 

TEMPFREQ=DELFREQ*DUMMYP 
TEMPFREQ=TEMPFREQ/I E+06 
ACTFREQ=JNINT(TEMPFREQ) 

TYPE*,’ ' 

TYPE *;POINT=',DUMMYP; FREQ~,ACTFREQ,’ MHZ' 

TYPE*,’ • 

RETURN 

END 


SUBROUTINE INTR (FI, F2, N, M) 
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IMPLICIT NONE 
C 

INTEGERS 1, M, N 
FNTEGER*2 F1(N), F2(M) 
REALM X, XA(512), Y, YA(512) 


DO I ~ 1, N 
XA(I) = FLOATI(I) 
END DO 


DO I - 1,N 

YA(I) = FLOATI(F I (I)) 
END DO 


DO I “ 1, M 

X = KLOATI(M) * FLOATI(N-l) / FLOATI(M-l) + 1,0 
CALL LINEARJNTERPOLATION (XA, YA, N, X, Y) 
F2(I) ~ ININT(Y) 

END DO 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


RETURN 

END 

LINEARJNTERPOLATION computes the value of the piecewise linear spline 

interpolant of the points in the arrays XA and YA at the point 

X. This routine is a modified version of the FORTRAN routine 

SPLINT found in the book "Numerical Recipes: The Art of 

Scientific Computing" by Press, Flannery, Teukolsky, and 

Vetterling published by the Cambridge University Press, 1986. 

Variables: 

NMAX : A constant which determines the array sizes 
for U,X, and Y. 

K : The midpoint between KH1 and KLO... used to find the two 
indices into the XA table such that 
XA(KLO) < X < XA(KHl). 

KHI ; The greater of the two indices into the XA array such 
that XA(KLO) < X < XA(KHI). 

KLO : The lesser of the two indices into the XA array such 
that XA(KLO) < X < XA(KHI). 

N : The index of the arrays and Y that contain interpolation 
points. 

A : Weight given to YA{KLO). 

B : Weight given to YA(KHI). 

H : The difference of XA(KHl) and XA(KLO). 

X : The value of X at which to compute the value of the 
interpolation function. 

XA : The table of X coordinates of the interpolation points. 

Y : The value of the interpolation function evaluated at X. 

YA : The table of Y coordinates of the interpolation points. 

SUBROUTINE LINEAR INTERPOLATION (XA, YA, N, X, Y) 


67 
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IMPLICIT NONE 

INTEGER NMAX 
PARAMETER (NMAX=512) 

INTEGERS K, KHL KLO, N 
REALM A, B, H 

REALM X, XA(NMAX), Y, YA(NMAX) 


KLO = 1 
KHI = N 

DO WHILE (KHI-KLO.GT.l) 

K - (KHI + KLO) / 2 
IF (XA(K).GT.X) THEN 
KHI = K 
ELSE 
KLO = K 
END IF 
END DO 

H = XA(KHI) - XA(KLO) 

IF (H,NE,0.0) THEN 

A = (XA(KHI) - X) / H 
B = (X - XA(KLO)) / H 
Y = A * YA(KLO) + B * YA(KIIl) 
END IF 


RETURN 

END 
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c ****** NOTHICK_DRAW.FOR *******************c 


INTEGER*2 TRFREQJXXXiXLENGTH 

INTEGER*2 SCANDIST,INDXDIST ; NSCAN,N1NDX,DENS 

!NTEGER*2 SELFREQ,AVES,]6J15,ISLANTUSLANT2,TSLANT3,ISLANT4 

CHARACTER RESP1 * 1 ,SELFREQANS* 1 ,RESP2*2 

CHARACTER NEWOLD*3,SHAPEANS* 1,Z1GANS*1 ,CHFREQ C*3,CHDRI*1 
BYTE B, CHFREQ(3), BSUFFREQ(4) 

BYTE B SCALEMARK(4) 

BYTE BDATAJSTATUS(40), BLABEL(30). BF1LEEXT(23) 

CHARACTER LABEL*3Q. FILEEXT*23,CHEADER*32,TEBAR*1 
CHARACTER L*48, DIR* 16, DRIVE* KNOTS AN S*LDIFFANS*1 
CHARACTER SUFFREQ*4,FTLENAME*34,SCALEMARK*4,DATA_STATUS*40 
CHARACTER ACTUAL_SELFREQANS* 1 ,S AS ANS* 1 .PHASE* 1 ,PH ASE 1 * 1 

EQUIVALENCE ( BFI LEEXT.FILEEXT ) 

EQUIVALENCE ( CHFREQ.CHFREQjC ) 

EQUIVALENCE (L,B) 

EQUIVALENCE ( BSCALEMARK, SCALEMARK ) 

EQUIVALENCE ( BSUFFREQ, SUFFREQ ) 

EQUIVALENCE ( BDATA STATUS, DATA STATUS ) 

EQUIVALENCE ( BLABEL, LABEL ) 

NE WOLD- OLD' 

OPEN(UNIT=10,FILE=’NOTHICK„DADQ1.DAT' ! STATUS=NEWOLD, 

I FORM- FORMATTED') 

READ (10,51 QFILEEXT 
READ (10,511) DRIVE 
511 FORMAT (A) 

51 FORMAT (A32) 

52 FORMAT (A2) 

53 FORMAT (A4) 

677 FORMAT (1) 


OPEN( unit=66, file- [ROTH. DATA] ’//FILEEXT/AdatC IT, status='OLD\ 
+ ACCESS- DIRECT',RECORDTYPE—FlXED , t RECL=32, 

+ form='FORMATTED\ORGANIZATION- SEQUENTIAL' ) 

OPEN (UNI T= 1 4 .FILE-'N OTHICK SELFREQ 1 .DAT, 

1 STATUS- OLD', FORM-FORMATTED') 

OPEN( unit- 15, file- [ROTH. DATA] Y/FILEEXT/Adatl 2', status='OLD\ 

I ACCESS - DIR ECE.RECORDTYPE- FIXED', RECL-2, 

+ form- FORM ATTED',ORG AN rZATION- SEQUENTIAL' ) 


16=1 

791 93 READ(66,REC=I6,FMT=5 1)CHEADER 
16=16+1 

RE AD(66,REC=I6,FMT=5 1 )SH APEAN S 
16=16+1 

READ(66,REC=I6,FMT=5 1 )ZIC,ANS 


5*7 
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16=16+1 


READ (14,98652) FILEEXT 
READ (14,98652) DRIVE 
READ (14, 98652) PHASE 
READ (1 4,98652) PHASE 1 
READ (14,98652) DIFFANS 
READ (14,10030) BUFLENGTH 
READ (14,10030) BUFVEL 
READ (14,10030) TRDIAM 
CLOSE (14) 

98652 FORMAT (A) 

10020 FORMAT (I) 

10030 FORMAT (F) 


115=1 

89714 READ( 1 5,REC=I 1 5,FMT=52)SCANDIST 
115=115+1 

READ( 1 5,REC-I1 5,FMT=52)INDXDIST 
115=115+1 

READ( 1 5,REC=I 1 5,FMT=52)N$CAN 
115=115+1 

READ( 1 5,REC=l 1 5,FMT=52)N1NDX 
115=115-1 

READ( 1 5 ,REC=1 1 5,FMT*=52)islant 1 
115=115+1 

READ( 1 5,REC=I 1 5,FMT=52)islant2 
U 5=1 15+1 

READ(15,REC=I 15,FMT=52)islant3 
115=115+1 

READ(1 5,REC=I1 5,FMT=52)islant4 
115=115+1 

READ( 1 5,REC=I 1 5 ,FMT=52)D£NS 
115=115+1 

READ( 1 5,REC=I 1 5,FMT=52)TRFREQ 
115=115+1 

READ( 15,REC=I 1 5,FMT=52)AVES 
115=115+1 

READ(I5,REC=I15,FMT=52)JXXX 

115=115+1 

READ( 1 5,REC=I 1 5,FMT=52)1XLENGTH 

C ARRAYSIZE={NSCAN)*(JXXX) 1IF SCAN FINISHED PROPERLY, JXXX=NINDX 


ENCODE (3,19878,CIIFREQ) TRFREQ ! ANALYSIS FREQUENCY 
19878 FORMAT (13) 


C DETERMrNE SHORTENED (-.DAT) FILEEXT POSITION 

C 

C KMARK=0 

C DO 4321 IPR=1, 34 

C IF (KMARK.EQ.l) GOTO 4321 


Go 
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C IF (RFILEEXT(IPR).EQ.46) THEN JPERIOD/DOT IS ASCII 46 

C JEND=IPR 

C KMARK=1 

C ENDIF 

C4321 CONTINUE 
C F1LEEXT=FILEEXT( 1 : IEND- 1 ) 

TYPE V CURRENT IMAGE FILE IS ',FILEEXT 
TYPE ’ 

TYPE V IF YOU WANT A DIFFERENT IMAGE FILE, YOU MUST* 
TYPE •; RECALL AN OLDER "NOTHICKDADQl .DAT" FILE' 

type v * 

IF (DRIVEJEQ.'A')DIR= , DUC2:[ROTH.IMAG]' 

IF (DRI VE.EQ .'C l )DIR~ DUCO : [ROTH . I MAG] ’ 

5321 TYPE VCURRENT DRIVE TO READ IMAGES IS ’ S DIR(1:5) 

IF (DIR(1 iS^EQ.'DUCO^THEN 
TYPE *,’DO YOU WANT TO CHANGE TO DUC2: (Y/N 
I )?’ 

READ (5,51)CHDRI 

IF (CHDRI.NE.'Y'AND.CHDRLNE. ! N , )GOTO 5321 
iF(CHDRI.EQ. , Y , )DlR= , DUC2:[ROTHJMAG]' 

ELSEIF (DIR(1 :5).EQ.'DUC2:')THEN 
TYPE *,’DO YOU WANT TO CHANGE TO DUCO: (Y/N 
I )?’ 

READ (5,5 1)CHDRI 

IF (CHDRLNE.Y\AND.CHDRI.NE. , N , )GOTO 5321 
IF (CHDRI.EQ.’Y^DIR-DUCO^ROTH.IMAG]' 

ENDIF 

2H TYPE*,'' 

TYPE V DO YOU WANT TO SEE TEXT AND COLORBAR (Y/N)?' 
READ (5,51)TEBAR 

IF (TEBAR.NE, r Y'AND.TEBAR.NE. 1 N')GOTO 21 1 
CONTINUE 

FILENAME- [ROTH.DATA]7/FILEEXT 

CALL GRINIT 
CALL GRSRST 


TYPE * 1 ' 

TYPE *, 'IMAGE FREQUENCY (MHz) = \TRFREQ 
ENCODE (3,9878.CHFREQ) TRFREQ 
9878 FORMAT (13) 

234 TYPE * ' ****** NOTHICKNESS VELOCITY IMAGE ****** ' 
TYPE */ ' 

WRITE (5,394) 

394 FORMAT 0$Y COLOR (C) OR B&W (B) IMAGE? ’) 

READ (5,2222) RESP1 




2222 

499 

500 

I 

444 

1 


212 

501 
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FORMAT(A) 

IF (RESPl.NE.'C'.AND.RESPl .NE.'B')GOTO 234 

[F (RESP1.EQ.'C)THEN 
TYPE V ' 

WRITE (5.499) 

FORMAT "COLORRED’ 1 (R) OR NORMAL COLOR (C) SCHEME? ’) 

READ (5,2222) RESP2 
ENDIF 

CALL IMAGE(NEWOLD,FlLEEXT,CHFR£Q_C,RESPl,RESP2,SHAPEANS,IXLENGTII, 
FILENAME, DIR, TEBAR) 


CONTINUE 

STOP 

END 


SUBROUTINE IMAGE(NEWOLD,FILEEXT,CHFREQ_C,RESPl,RESP2,SHAPEANS,lXLENGTH, 
FILENAME, DIR, TEBAR) 

INTEGER *2 D(460), CHANNEL, INPUT, IXLENGTH 

INTEGER*2 X4,Y4,N 1 ,N2,TRFREQ,SCANDIST 

CHARACTER 5 * 1 RESP 1 * 1 ,SHAPEANS* I ,RESP2* 1 .TEBAR 

CHARACTER CHFREQ_C*3 

REAL RMIN, RMAX, MEAN 

BYTE B, CHFREQ(3), BSUFFREQ(4),BZERO(2) 

BYTE BFILENAME(34), BSCALEMARK(4) 

BYTE BDATA_STATU S(40), BLABEL(30) 

BYTE CHMEAN(5),CHRMAX(5),CHRMIN(5) 

CHARACTER FILNAM*48, LABEL*30, FILEEXT*23 
CHARACTER L*48, DIR*16 : CZERO*2 

CHARACTER SUFFREQ*4, FILENAME*34. SCALEMARK*4, DATA_STATUS*40 

EQUIVALENCE (L,B) 

EQUIVALENCE ( BSCALEMARK, SCALEMARK ) 

EQUIVALENCE ( BSUFFREQ, SUFFREQ ) 

EQUIVALENCE ( BDATAjSTATOS, DATA STATUS ) 

EQUIVALENCE ( BLABEL, LABEL ) 

EQUIVALENCE ( CZERO, BZERO ) 

CZERO-O 1 

LLLL-0 

LIIMAGE=1 

CONTINUE 

IF (IIIMAGE.EQ. 1 )FILNAM= T _VEL_CJ//CHFREQ_C/^.PDK’ 

CONTINUE 

CHANNEL=0 

IF (CHANNEL.NE.O . AND . CH ANNEL.NE .2. AN D . CH ANNEL .NE.3 ) GOTO 501 
IF (CHANNEL.EQ.O) THEN 
INPUT=0 


C,x 



MASK=l 

ELSEIF (CHANNEL.EQ.2) THEN 

INPUT=1 

MASK-4 

ELSETF(CHANNEL.EQ.3) THEN 

INPUT=2 

MASK=8 

ENDIF 

CALL ERASE 


1 9 OPEN (l JNIT-6,FILE=DI R// FILEEXT//F I LN AM, STATUS- OLD’, 

1 FORM='UNFORM ATTED 1 ) 

READ(6) X4,Y4,N1,N2J VALUE 
READ(6) SCANDIST,TRFREQ 
READ(6) RMIN,RMAX,MEAN 

IF ( IVALUE.EQ.I ) LABEL- VELOCITY CM/US' 

ICOLORBAR=0 
SUFFREQ- MHZ' 

SCALEMARK- 1 MM* 

DAT A_STATUS=LA BEL 

CALL GRNrN(l,INPUT,rNPUT,INPUT) 

CALL GRNBY(1 ,1,1,1) 

CALL GRSBFD 


c ****************h.***** MO dIFICATION FOR COMPLEX SHAPES ******C 
IF (SHAPEANS.NE.’Y')IX49-X4 

IF (SHAPEANS.EQ.'Y')THEN 
1X49=1 XLENGTH 
ENDIF 

DO 100 1=1,1X49 
IF (SHAPEANS.NE.’Y')THEN 
IXSTART=I-I 
IYSTART=1 
IYLENGTH-Y4 
READ(6)D 
GOTO 789 

ELSEIF (SHAPEANS.EQ.Y')THEN 
IF (I.NE. 1 .AND.IYLENGTH.EQ.1 )GOTO 100 
READ(6,END= 1 00)DCSTART,I YSTARTJYLENGTH 
DO 7897 1JK=1,IYLENGTH 
7897 READ(6)D(UK) 

789 CONTINUE 
ENDIF 

CALL ORW A W(D,IXSTART,I YSTART, 1 JYLEN GTH, 0,0.0, 1,1) 

Q ****** ********************************************** **>♦****+ Q 

CALL GRSBFD 
100 CONTINUE 
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IF(RESP!.EQ:B')THEN 
IF (TEBAR.EQ.'N')GOTO 213 
CALL PLACEBAR(O) 

GOTO 8999 

ELSEIF (RESPl.NR'B'JTHEN 

IF (RESP2.EQ.'C’)THEN 

CALL COLORS (!VALUE,ICOLORBAR) 

ELSEIF (RESP2.EQ. , R')TIIEN 

CALL COLORRED (IV A LUEJCOLORB AR) 

ENDIF 

ENDIF 


IF (TEBAR.EQ.N^GOTO 213 


8999 DECODE (34,8989,FILEN AME) BFILENAME 
8989 FORMAT (34A1) 

ENCODE (3,909,CHFREQ) TRFREQ 
909 FORMAT (13) 


CALL GRFCD (1,255,0,0, BFILENAME,!, 44 1,6,0,34, 0,0,0) 
CALL GRFCD (1 ,255.0, 0,CHFREQ, 200, 44 1, 6,0,3, 0,0,0) 
CALL GRFCD ( 1 ,255,0,0,BSUFFREQ,2 1 5,44 1 ,6,0,4,0,0,0) 

CC 

CC PRINT LABEL TO GRINNELL 

CC 


NCHAR~LEN(DATA_STATUS) 

DO 765 1KY=1,NCHAR 
IKX=400-(IKY-1)* 13 

CALL GRFCD (1.255, 0,0, BDATA STATUS(IKY),505.1KX,0,0J, 0,0,0) 
CALL GRSBFD 
765 CONTINUE 

CC 

CC CALCULATE MAX, MIN, AND MEAN 
CC 


IF (IVALUE.EQ.I.0R.IVALUE.EQ.2.0R.IVALUE.EQ.6.0R.IVALUE.EQ.7) THEN 
IF (IVALUE.EQ.l) THEN 

IF (RMAXJEQ.1.0) RMAX- 9999 
ENDIF 

RRMAX=RMAX* 10.**4. 

RRMIN=RMIN*10.**4. 

RMEAN=MEAN*10.**4. 


ELSEIF (IVALUE.EQ.3.0R.IVALUE.EQ.4) THEN 
RRMAX-RMAX/100. 
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RRMFN=RMIN/ 1 00. 

RME AN=ME AN/ 1 00. 

ELSEIF (IVALUE.EQ.5) THEN 
IF (RMAX.GE.100.)RMAX=99.9 
IF (RM1N.GE.1 00.)RM1N=99.9 
IF (MEAN.GE.I00.)MEAN=99.9 
IF (RMAX.LT.IO.)TTON 
RRMAX=RMAX*10.**4. 

RRMIN=RMIN* 1 0.**4. 

RMEAN=MEAN* 10**4. 

ELSE 

RRMAX=RMAX*10.**3. 

RRMIN=RMIN*I0.**3. 

RMEAN=MEAN* 10**3. 

ENDIF 

ENDIF 

IF (IVALUE.EQ.2.0R.IVALUE.EQ.6.0R.1VALUEJEQ.7) THEN 
JKMN=465 

ELSEIF (IVALUE.EQ. 1 .OR .IV ALUE.EQ.3 .OR.I VALUE.EQ.4)THEN 
JKMN=463 

ELSEIF (IVALUE.EQ.5)THEN 
JKMN=470 
ENDIF 

CC DRAW DECIMAL POINT FOR MIN AND MAX TO GR1NNELL 

DO 9879 IUY=0, 350, 350 
CALL GRFAR (1,255,0,0,JKMN,IUY,1 ; 1) 

CALL GRSBFD 
9879 CONTINUE 


IMEAN=JN1NT(RMEAN) 

IRMAX=JNINT(RRMAX) 

IRMIN=JNINT (RRMIN) 

TYPE V ’ 

TYPE *, ’CORRECTED:' 

TYPE *,’1RMAX-,1RMAX 
TYPE *,'IRMIN-,IRM1N 
TYPE *,’IMEAN=’,IMEAN 
type*,’ 1 

IF (IMEAN.GE.1.0*10.**5.) GOTO 73920 
ENCODE (5,9878,CHMEAN) IMEAN 

73920 IF (IRMAX.GE.1.0* 10.**5.) GOTO 73921 
ENCODE (5,9878, CHRMAX) 1RMAX 

73921 IF (IRMIN.GE.1.0*10.**5.) GOTO 73922 
ENCODE (5,9878,CHRMIN) IRMFN 

73922 CONTINUE 
9878 FORMAT (15) 


65 ''' 
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IF (IMEAN.GE. 1 .0* 10.**6.)GOTO 93920 
ZZZZZZ=RRMAX-RRMIN 
IF (ZZZZZZ.EQ.0.)ZZZZZZ=1. 

SC ALEF=((RME AN -RRMIN)/ (ZZZZZZ)) 4 3 50 . !SC. MEAN TO CLR BAR 
ISCALEF=ININT(SCALEF) 

93920 CONTINUE 

IF (IMEAN.GE, 1 .0* 1 0.**6.)GOTO 83920 

CALL GRFCDS (1,255,CHMEAN,458,ISCALEF,6,0,4) 

83920 IF (IRMAX.GE. 1.0*1 0.**6.)GOTO 83921 
CALL GRFCDS (l,255,CHRMAX t 458, 350,6,0, 4) 

83921 IF (IRMIN.GE. LO* 1 0. **6.)GOTO 83922 
CALL GRFCDS ( 1 ,25 5,CHRMIN,45 8,0,6, 0,4) 

CALL GRSBFD 

83922 CONTINUE 
9877 CONTINUE 


IF (IVALUE.EQ.2) THEN 
JKMN=465 

ELSEIF (IVALUE.EQ.LOR.JVALUE.EQ.3.0R.IVALUE.EQAOR.TVALUE.EQ.6 
1 .OR.IVALUE.EQ.7)THEN 

JKMN=463 

ELSEIF (IV ALUE.EQ. 5 )THEN 
JKMN=470 
ENDIF 

CALL GRFAR (1,255,0,0,JKMN,ISCALEF,1,I) 

CALL GRSBFD 
95631 CONTINUE 

CC 

CC DRAW SCALE ONTO IMAGE 
CC 


RSCANDIST=FLOATI(SCANDIST) 

PIXELS=(X4*1 ,)/(RSCANDIST) 

IP1XELS=ININT(P1XELS) 

CALL GRFVCS ( 1 ,255,340,420,340+IPlXELS,420) 

CALL GRFCDS (1,255,BSCALEMARK, 345,410, 6,0,4) 

CALL GRSBFD 

C! FOR DRAWING INITIAL ZEROES FOR ATTENUATION COEFF. < 0.1 
IF (I VALUE.EQ.2. AND.RRMIN .IX 1 000.)THEN 
CALL GRFCDS (1 >255,BZERO,JKMN,0,6,0 S 2) 

CALL GRSBFD 
ENDIF 


213 CONTINUE 
201 CLOSE(6) 


1 FORM ATf$V WHICH IMAGE TO VIEW [*.PDK IN ROTH.IMAGJ?’) 

2 FORMAT(A) 

3 FORMATCS', 'Channel (0,2,3) ?■) 


u 
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4 FORMAT( 'SDescription ?? ' ) 

5 FORMAT('O') 

6 FORMAT( 1 ',A ) 

7 FORMAT( ' Range = ',A,’ , ',A ) 

8 FORMATC Mean = \A15/ \A) 

9 FORMAT(I) 


TYPE V ’ 


999 CALL GRSEND 

444 CONTINUE 
RETURN 
END 


SUBROUTINE COLORSUB (IVALUEJCOLORBAR) 

INTEGERS NR(256),NB(256),NG(256),IARRAY(1 8) 
DATA I ARRAY/0,2, 1,0, 1,2, 1 ,0.2, 1 ,2 ,0,2,0, 1 ,2, 1,0/ 

IF (IVALUE.GT. 10) GOTO 9876 
CALL GRINITOLD 
CALL GRSRST 
CALL GRZFC(0,0) 

CALL GRZCL(0,255,255) 

CALL GRZON(0,1) 

CALL GRSBFD 


B=4 

A=2 


L=256 

IL=L/B 

JL=A*IL 


DO 100 1=1, IL 
NB(l)-255 
NR(I)=0 
NG(I)=B*I 

IF(NG(I).GT.255)THEN 

NG(I)=255 

ENDTF 

100 CONTINUE 

DO 200 I=IL+1,JL 

BB=((B/( 1 -A))*I)+((A/(A-1))*L) 

NG(I)=255 

NR(I)=0 

NB(I)=BB 

IF (NB(1).GT.255)THEN 

NB(I)=255 

ENDIF 

200 CONTINUE 


G1 
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DO 300 1=JL+UL+IL 

RR=(B*I)-(A*L) 

NG(I)=255 

NB(I)=0 

NR(I)=RR 

IF (NR(I).GT.255)THEN 

NR(1)=255 

ENDIF 

300 CONTINUE 

DO 400 I=JL+IL+1,L 

GG=((-B)*I)+(B*L) 

NR(I)=255 

NB(I)=0 

NG(I)=GG 

IF(NG(I).GT.255)THEN 

NG(I)=255 

ENDIF 

400 CONTINUE 


800 KC"0 


NR(1)=0 

NG(1)=0 

NB(1)=0 

NR(256)=255 

NG(256)=255 

NB(256)=255 

889 CALL GRNWR(I,NB,1ARRAY(1+KC),0.256,0) 

CALL GRNWR(I,NGJARRAY{2+KC), 0,256,0) 

CALL GRN WR{ 1 ,NR, LARRAY (3+KC),0,25 6 ,0) 

CALL GRSBFD 
CALL GRNIN(0,0,0,0) 

IF (ICOLORBAR.EQ. 1) THEN 
CALL GRNBY (0,0, 0,0) 

ENDTF 

CALL GRSBFD 

C IF (IVALUE.EQ. I .OR.1VALUE.EQ.3.0RTVALUE.EQ.5.0R.JVALUE.EQ.7) THEN 

DO 1111 1=1,300 
DO 1111 J=l,900 
1 1 1 1 CONTINUE 
GOTO 950 
950 KC=0 
9876 RETURN 
END 


SUBROUTINE PLACEBAR (ICOLORBAR) 
CC Draw colorbar at right side of screen 




INTEGER* 2 NA(45) 
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IF (ICQLORBAR.EQ.l) GOTO 1010 

CONTINUE 

DO 1000 1-0,350 

DO 2000 J=l,19 

G=l*(254.0/350.0) 

NA(J)=G 

2000 CONTINUE 

CALL GRWLW(NA,41 0,1, 1 9,0,0, 1,1) 
CALL GRSBFD 
1000 CONTINUE 
1010 CONTINUE 
RETURN 
END 


SUBROUTINE SAVEBAR 
INTEGER*2 BAR(6) 

TYPE * ' ’ 

TYPE *, 'Creating and saving the color bar.,/ 
TYPE * ' ’ 


OPEN(UNTT=34.NAME=’BARFILE', STATUS-NEW 1 , 
1 FORM- UNF ORMATTED\BLOCKS I Z £- 1 024) 

CALL GRNIN(1, 0,0,0) 

CALL GRNBY( 1,1.1, 1) 

CALL GRSBFD 

DO 2060 JBBEL-0,459,3 
CALL GRRLW (4I0,IBBEL,6,1,BAR) 

WR1TE(34) BAR 
2060 CONTINUE 

CALL GRSBFD 
CLOSE(34) 

RETURN 

END 


SUBROUTINE ERASE 

M-4095 
N— 4095 

CALL GRFER (M,N,0) 
CALL GRSBFD 
END 


SUBROUTINE COLORS ( I V AL U E,ICOLORB A R) 

a 
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INTEGERS NR(256),NB(256),NG(256),IARRAY(1 8) 
DATA I ARRAY/0,2, 1 ,0, 1 ,2, 1 ,0,2. 1 ,2,0,2 ,0, 1,2. 1 ,0/ 
INTEGER*2 NA(45) 


IF (IVALUE.GT.10) GOTO 9876 
CALL GR1NITGLD 
CALL GRSRST 
CALL GRZFC(0,0) 

CALL GRZCL(0,255.255) 

CALL GRZON(0, 1 ) 

CALL GRSBFD 


B =4 

A=2 


L=256 

IL=L/B 

JL=A*IL 


DO 100 1=1, IL 
NB(I)=255 
NR(I)=0 
NG(I)=B*I 

1F(NG(I).GT.255)THEN 

NG(I)=255 

ENDIF 

100 CONTINUE 

DO 200 I=IL-t-l,JL 

BB=((B/(1-A))*I)+((A/(A- 1))*L) 

NG(I)=255 

NR(I)=0 

NB(I)=BB 

IF (NB(I).GT.255)THEN 

NB(1)=255 

ENDTF 


200 CONTINUE 


DO 300 I“JL+1,JL+IL 

RR=(B*I)-(A*L) 

NG(I)=255 

NB(I)=0 

NR(T)=RR 

IF (NR(I).GT.255)THEN 

NR(I)=255 

ENDIF 

300 CONTINUE 

DO 400 I=JL+IL+ 1 ,L 
G G=((-B) * I)+(B *L) 
NR(I)=255 


70 
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NB(I)=0 

NG(I)=GG 

IF(NG(I).GT\255)THEN 

NG(I)=255 

END1F 

400 CONTINUE 


800 KC=0 


NR(I)-0 

NG(1)=0 

NB(1)=0 

NR(256)=255 

NG(256)=255 

NB(256)-255 

889 CALL GRNWR(1 ,NB,IARRAY(M-KC),0, 256,0) 
CALL GRNWR( l ,NG ,1 ARRA Y(2+KC), 0,256,0) 
CALL GRNWR(LNR,IARRAY(3+KC), 0,256.0) 
CALL GRSBFD 
CALL GRMN(0,0,0,0) 

IF (ICOLORB AR.EQ. 1 ) THEN 
CALL GRNBY (0,0, 0,0) 

ENDIF 

CALL GRSBFD 
00 11111=1,300 
DO 1111 J- 1,900 

111! CONTINUE 
GOTO 950 

950 KC=0 


CC Draw colorbar at right side of screen-- 

9876 IF (ICOLORBAR.EQ.l) GOTO 1010 
CONTINUE 
DO 1000 1=0,350 
DO 2000 3=1,19 
G=I *(254. 0/3 50.0) 

NA(J)=G 

2000 CONTINUE 

CALL GRWLW(NA,4 1 0,1, 1 9,0,0, 1.1) 
CALL GRSBFD 
1000 CONTINUE 
1010 CONTINUE 
RETURN 
END 


Subroutine colorred(l VALUE, ICOLORBAR) 

INTEGERS NR(256),NB(256),NG(256),IARRAY( 1 8) 
DATA I ARRAY/0, 1,2, 0,2, L 1 ,0 ,2, 1 ,2, 0,2,0, 1 ,2, 1 ,0/ 
1NTEGER*2 NA(45) 


7 / 
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100JO format(A) 
10020 format(I) 

NM256 = 0 


IF (IVALUE.GT.10) GOTO 9876 

CALL GRINITOLD 
C CALL GRSRST 
CALL GRSBFD 
L=256 
IL=L/3 
JL=2*1L 
NM256=255 

DO 100 1=1, 1L 
NB(I)=0 
NR(I)=I*3 
NG(I)=0 

IF( NR(I).GT.253 )NR(1)=253 
100 CONTINUE 

DO 200 I=1L+1,JL 
NG(I)-3*(I-IL) 

NR(I)=253 

NB(I)=0 

IF( NG(1).GT.253 )NG(I)=253 
200 CONTINUE 

DO 300 I=JLH,255 
NG(I)=253 
NR(I)=253 
NB(I)=3*(I-JL) 

IF( NR(I).GT.253 )NR(I)=253 
300 CONTINUE 

800 KC=6 

NR{1)=0 
NG( 1 )=0 
NB(1)=0 

NR(256)=NM256 

NG(256)=NM256 

NB(256)~NM256 

840 CALL GRN WR( 1 ,NB,1ARRAY(1 +KC),0,256,0) 
CALL GRNWR( 1 ,NG,IARRAY (2+KC),0.256,0) 
CALL GRNWR(1 ,NRJ ARRAY(3-KC) 1 0,256,0) 
CALL GRSBFD 
CALL GRNIN(0,0, 0,0) 

CALL GRSBFD 
CALL GRNBY(0,0,0,0) 

CALL GRSBFD 
C WRITE( 5,10860) 

Cl 0860 format( '$ <CR> for next color scheme, Q to quit ) 
C READ( 5,10010 )ANS 


7a. 
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C 1F( ANS.EQ.’Q' .OR. ANS.EQ.'q' )GOTO 9876 
GOTO 9876 
KC = KC+3 
IF (KC.GE.18) KC-0 
GOTO 840 


CC Draw colorbar at right side of screen- 

9876 IF (ICOLORB AR.EQ. 1 ) GOTO 1010 

CONTINUE 
DO 1000 1=0,350 
DO 2000 J=l,19 
G=I *(254.0/3 5 0.0) 

NA(J)=G 

2000 CONTINUE 

CALL GRWLW(NA,41 0,1, 19,0,0, 1,1) 
CALL GRSBFD 
1000 CONTINUE 
1010 CONTINUE 


C 900 CALL GRSEND 
return 
END 


73 
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What is claimed is: 

1. A pulse-echo, immersion method for ultrasonic evalu- 
ation of a material, employing automatic scanning and 
digital imaging to obtain an image of a property of said 
material, wherein said material is held in a holding apparatus 
which is positioned in an immersion liquid over an acoustic 
reflector, said reflector having an acoustic impedance which 
is greater than that of said liquid, and wherein nonlevelness 
in said holding apparatus and material thickness axe 
accounted for and eliminated, said method comprising: 

(i) ultrasonically scanning said material at a plurality of 
scan points and receiving the first and second echoes, 
each of which is a complete waveform, reflected off the 
back surface of said material and the first echo reflected 
off the front surface of said reflector both with and 
without the presence of said material; 

(ii) adjusting the time delay for each said received echo 
from each said scan point during said scanning in (i) 
above, gating each said received echo so that it is 
centered within its respective time window; 

(iii) automatically scanning said material at said scan 
points to receive said first and second back surface 
echoes and said two front surface echoes using the 
information obtained in (ii) above, so that each echo 
received from each scan point during said automatic 
scanning is centered within its time window; 

(iv) digitizing each echo received during said automatic 
scanning and determining the time delay between said 
first two successive sample back surface echoes, 2t, 
and the time delay, At, between the two different 
reflector front surface echoes received at each scan 
point during said automatic scanning and calculating 
the wave velocity, using a cross correlation function, at 
each said scan point from 

v=c (-§- +1 ) 

where c is the speed of the ultrasonic wave transmitted 
in said liquid, and 

(v) scaling the velocity values obtained in (iv) to corre- 
sponding proportional color or grey scale values and 
displaying the resulting image. 

2. A method according to claim 1 wherein a single 
transducer is used. 

3. A method according to claim 2 wherein said transducer 
is a high frequency transducer which emits a frequency 
between 1-100 MHz. 

4. A pulse-echo, immersion method for ultrasonic evalu- 
ation of a material employing a single transducer, automatic 
scanning and digital imaging to obtain an image of a 
property of said material, wherein said material has a 
uniform thickness variation and is positioned in an immer- 
sion liquid between said transducer and an accoustic 
reflector, said method comprising: 

(I) accounting for and eliminating nonlevelness in the 
set-up and said material thickness variation by; 

(a) performing a preliminaiy scan along both the x-and 
y-directions of the material to provide slant correc- 
tion factors which are input into a computer to 
account for said nonlevelness and thickness variation 
during the subsequent automatic scanning for said 
material evaluation in (ii) below; 

(b) adjusting the time delay during said preliminary 
scan for any received echoes which are not centered 
in the scan time window, so that each received echo 
is centered in its time window; 
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(ii) automatically scanning said material at a plurality of 
scan points in both the x- and y-directions to receive the 
first and second back surface echoes and front surface 
echoes with and without the presence of said material 

5 between said transducer and reflector using the infor- 
mation obtained in (i) above, so that each echo received 
from each scan point during said automatic scanning is 
centered within its time window; 

(iii) digitizing each echo received during said automatic 

10 scanning and determining the time delay between said 

first two successive sample back surface echoes, 2t, 
and the time delay, At, between the two different 
reflector front surface echoes received at each scan 
point during said automatic scanning and calculating 

15 the wave velocity at each said scan point from 

v=c (-§r + 1 ) 

20 

where c is the speed of the ultrasonic wave transmitted 
in said liquid, and 

(iv) scaling the velocity values obtained in (iii) to corre- 
sponding proportional color or grey scale values and 

25 displaying the resulting image. 

5. A method according to claim 4 wherein said back and 
front surface echoes received from the first and last scan 
points in both the x- and y-directions during said preliminary 
scan determine said time base adjustments needed for each 

30 echo to be centered within the time frame for it. 

6. A method according to claim 5 wherein the location of 
the time window during said automatic scanning for said 
material evaluation is automatically adjusted via computer 
control by using the formula: 

wherein W^ r is the correct delay time window at a particular 
scan location, T I is the time delay at the the initial scan 
location, X sc and Y sc are the x- and y-direction slant 

40 correction factors, X SN and Y SN are the scan point numbers 
in the x- and y-directions, and X SI and Y s/ are the x- and 
y-direction scan increments. 

7. A method according to claim 6 wherein two scans are 
automatically made to obtain said first two back surface 

45 echoes and said two different reflector front surface echoes. 

8. A method according to claim 7 wherein said first two 
back surface echoes and said reflector echo with said mate- 
rial present are made in said first scan. 

9. A method according to claim 6 wherein three scans are 

50 automatically made to obtain said first two back surface 

echoes and said two different reflector front surface echoes. 

10. A method according to claim 9 wherein said first two 
back surface echoes are received in one scan, wherein said 
reflector echo with said material present is made in another 

55 scan, and wherein and said reflector echo without said 
material present is received in yet another scan. 

11. A method according to claim 9 wherein said first back 
surface echo is received in said first scan, wherein said 
second back surface echo is received in said second scan, 

60 wherein said reflector echo with said material present is 
made in said third scan, and wherein said reflector echo 
without said material present is received in said fourth scan. 

12. A method according to claim 6 wherein four scans are 
automatically made to obtain said first two back surface 

65 echoes and said two different reflector front surface echoes. 

13. An ultrasonic, pulse-echo, immersion method employ- 
ing automatic scanning and digital imaging to obtain an 
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image of a microstructural property of a material positioned 
in an immersion liquid between a transducer and an acoustic 
reflector, said method comprising: 

(i) automatically scanning said material at least three 
times at a plurality of scan points in both the x- and 
y-directions to receive the first and second back surface 
echoes and front surface echoes with and without the 
presence of said material between said transducer and 
reflector, each of said echoes received being a complete 
waveform and gated within a time window; 

(ii) digitizing each echo received during said automatic 
scanning and determining the time delay between said 
first two successive sample back surface echoes, 2t and 
the time delay, At, between the two different reflector 
front surface echoes received at each scan point during 
said automatic scanning and calculating the wave 
velocity at each said scan point from 

v = c (-§- +1 ) 
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where c is the speed of the ultrasonic wave transmitted in 
said liquid, and 

(iii) scaling the velocity values obtained in (ii) to corre- 

5 sponding proportional color or grey scale values and 
displaying the resulting image. 

14. A method according to claim 13 wherein four scans 
are automatically made to obtain said first two back surface 

10 echoes and said two different reflector front surface echoes. 

15. A method according to claim 14 wherein said first 
back surface echo is received in said first scan, wherein said 

15 second back surface echo is received in said second scan, 
wherein said reflector echo with said material present is 
made in said third scan, and wherein and said reflector echo 
without said material present is received in said fourth scan. 

20 


* * * * * 



